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Chapter 1 

Introduction 



One of the major achievements in theoretical physics in the past 40 years was to find out that 
not only the electrowcak interactions can be properly described through gauge theories, but 
that the same principle is obeyed by the strong interactions. The reason why QCD had for 
so long eluded theorists' efforts is mainly due to its most peculiar distinguishing feature, i.e., 
confinement. Spectroscopy experiments in the 1960's already made it clear that the zoo of new 
particles discovered at accelerators could not be elementary. The successful Eightfold Way of 
Gell-Mann and Ne'eman strongly relied on the existence of quarks as constituents, but even 
for Gell-Mann at that time their physical meaning was more than dubious. However, in the 
long run quarks were accepted to be elementary particles, which, endowed with a dynamical 
SU(3)c gauge group, build up the so-called QCD Lagrangian, which is purported to yield the 
proper description of the strong interactions. Confinement is then responsible for the fact that 
a quark or gluon cannot be detected in isolation, but only colourless combinations of them, i.e., 
hadrons. Nowadays we still lack a proof of confinement (we only have some persuasive hints 
from lattice simulations) , but the whole scheme has been thoroughly and successfully tested in 
several experiments. 

The fact that the QCD Lagrangian is expressed in terms of quarks and gluons instead of 
hadrons makes the determination of physical observables inside the Standard Model a highly 
non-trivial task. At high energies, asymptotic freedom renders the strong coupling constant 
small enough to allow the use of perturbation theory methods. However, below 1 GeV, confine- 
ment enters the game and binds quarks and gluons together, thus requiring the introduction 
of non-perturbative techniques. Common strategies are to resort to models (Nambu-Jona- 
Lasinio, Constituent Quark Model, ...), to rely on numerical simulations (lattice QCD) or to 
use approximations to QCD (large- iV c QCD). 

Among the interesting issues that need to be addressed with the above-mentioned tech- 
niques, an oustanding one is to measure the amount of CP violation provided by the Standard 
Model as compared with the one observed experimentally. Since its discovery more than 50 
years ago, the source of CP violation still remains an open field in particle physics. The pre- 
ferred scenario to test the Standard Model has been for years kaon physics. More recently, the 
same phenomena has been observed in B physics, opening an era of high precision measure- 
ments. 

The present work is devoted to determining some of the most relevant parameters that 
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account for CP violation in the Standard Model in a model independent way. Our framework 
will be restricted to kaon physics, even though our methodology can easily be extended to B 
physics. Our analysis relies on the use of an approximation to large- N c QCD (coined Mini- 
mal Hadronic Approximation), to deal with the non-perturbative aspects inside kaon matrix 
elements. Matching between short and long distances is an issue in kaon weak interactions, 
which our approach solves in a very natural way. A key ingredient therefore is to provide the 
right matching condition onto the OPE of QCD. This is precisely how scale dependencies are 
removed from physical calculations, and no spurious cut-off shows up. As far as we know, 
this is the only existing technique which allows this consistent scale and scheme-independent 
determination of hadronic parameters. 

The outine of the present work will be as follows: in chapter 2 we give a brief reminder of the 
Standard Model together with its low energy chiral realizations in the strong and clcctroweak 
sectors. 

Chapter 3 introduces the 1/Nc expansion and then deals with extensions of the chiral 
Lagrangian to higher energies through the introduction of explicit resonance fields, so-called 
Resonance Chiral Lagrangians. We show how quantum corrections to such a Lagrangian can be 
performed provided we use the power counting rule coming from large- Nc QCD. We compute, 
as an example, the low energy coupling Liq and test thereby the validity of the Lagrangian. 

Chapter 4 is devoted to a review of kaon phenomenology and the introduction of the relevant 
parameters in the study of CP violation. Determination of these phenomenological parameters 
is the subject of chapter 5 and the main core of the present work: we compute Bk and 527 
in the chiral limit and to subleading order in the 1/N C expansion. Furthermore, we reassess 
the determination of the kaon mass difference ArriK and the parameter ek in the large- N c and 
chiral limits when corrections in inverse powers of the charm mass are taken into account. 

In chapter 6 we move to another subject and deal with the issue of quark-hadron duality 
violations, which arise due to the fact that the OPE is not a valid expansion in the whole 
complex q 2 plane, but breaks down (at least) in the Minkowski half real axis. Determination of 
the OPE coefficients using spectral sum rules is then possible up to some inherent uncertainty: 
the duality violating piece, whose effect is systematically ignored in all determinations of OPE 
coefficients. There is no theory behind duality violations and one has to resort to models to 
study them. Armed with a toy model of large- N c inspiration, we assess the amount of duality 
violation in a particular QCD two-point Green function, the < VV — AA > correlator, whose 
OPE condensates are relevant for kaon decays, looking for the best strategy to extract the OPE 
coefficients reliably. Finally, we end with the conclusions. 



Chapter 2 

The Standard Model at Low 
Energies 

The Standard Model has been on the market for the last 40 years and many good reviews 
on the subject exist. We will here skip the many details and content ourselves with a brief 
overview designed to serve as a quick reference. For more comprehensive accounts we refer the 
reader to, e.g., PE]. 

The Standard Model is a Yang-Mills theory invariant under the gauge group 

SU(S) C x SU(2) L x U(1) Y (2.1) 

The first factor SU(3)c accounts for the strong interactions, whose discussion will be the 
subject of section (2.2). The remaining SU (2)l X U (l)y factor is responsible for the electroweak 
interactions. We will first discuss the inner structure of the electroweak interactions, to move 
afterwards to the strong sector. 



2.1 The Electroweak Sector 

The electroweak interactions merge the weak and electromagnetic interactions together under 
the unified SU(2)l X U(l)y gauge group. This merging of the two interactions is however 
highly non-trivial. The notion of spontaneous symmetry breaking, to be introduced in short, is 
the mechanism by which one can single out and recover both interactions. Before we proceed, 
it is convenient to expand the Dirac fields for the fcrmions in the heliciy basis using the hclicity 
projectors, 1 

Pl=(^) , Pr^(^) (2-2) 



1 Indeed, it is straightforward to check that they verify 



Pl + Pr = 1 , Pl ■ Pr = , P( L<R) = P (L , R) 
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and introduce the short-hand notation, to be used hereafter, 

1-75' 



tpL = Pl 4> 



1 + 75 



(2.3) 



The matter content of the Standard Model consists of three families of lcptons and three 
families of quarks, to be organized in three generations as follows 
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together with the gauge vector bosons 

W + W~ 
and the would-be intriguing Higgs scalar boson 



-4,, 



H 



(2.4) 



(2.5) 



(2.6) 



Equation l|2.4|) shows that left-handed fcrmions are isodoublcts of SU(2)l, whereas the right- 
handed ones are isosinglets. Defining 



L, 



R e = e^j , R v = v 



R ' 



Ru = UR , Rd = d-R 



we can write the free fermion sector of the Standard Model as follows 



£*Dirac ^ ^ ^ 



Lk^d^Lk + Rk-fdpRk 



(2.7) 



(2.8) 



where k is to be understood as running over all multiplets. Interaction terms can then be 
read off by invoking the gauge principle. The usual trick is to promote ordinary derivatives to 
covariant ones, 
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(2.9) 
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where the superscripts refer to the different chiralitics. t j are the three generators of the weak 
isospin group SU(2)l in the form of the Pauli matrices, listed in Appendix A, and Y is the 
weak hypercharge, the generator of the U(l)y gauge group. and arc auxiliary fields, 
the so-called gauge fields, whose kinematical term is 

C Kin = -\ W$ Wfi - \ B^ (2.10) 
where the field strength tensors are defined through 

= d„wi-d v w;- 9 ^ k w^w!; 

B^ = dp B v — d u Bft (2.11) 

Once expanded, the covariant derivatives in the Dirac Lagrangian ({2.8(1 will generate the in- 
teractions between quarks and leptons through the gauge bosons. With the use of Noether 
theorem, we can determine the (conserved) isospin and hypercharge currents, 



J1 = E L kl^ y L k 



- -E 



(2.12) 



Whereas r 3 is diagonal, t 1 and r 2 are antidiagonal and give rise to charged currents. The 
interaction Lagrangian for weak charged currents reads 

C cc = g(J?Wi + J$W£) 

= JL (j^W^ + JHW-^ (2.13) 

where in the second line diagonalisation has been performed so as to single out each charged 
current. The new current basis is 

j£=2(J?^iJ$) (2.14) 
with their associated gauge bosons given by 



W+\ , / 1 -i 



W~ J vz \ 1 +i 




(2.15) 



Likewise, neutral currents appear in the Standard Model in the following manner 



Cnc = gJ^Wl + ^-gJ^B,, 



j r A ^^4^ (2-16) 



Again, in the second line we have rotated to the physical basis, 

Zp\ I cos0w -sin0 w \ / W* 

An J \ suvOw cos 9 W J \ B^ 



(2.17) 
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T 3 


1/2 


-1/2 





1/2 


-1/2 
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-1 


-1 


-2 


1/3 


1/3 


4/3 


-2/3 


Q 





-1 


-1 


2/3 


-1/3 


2/3 


-1/3 



Tabic 2.1: Quantum numbers of the particle content of the Standard Model. For simplicity, 
only the first generation is shown. 



where 9w is the so-called electroweak mixing angle, defined as 

cos0 w = -f==f (2.18) 

V.9 2 + g 2 

This way, can now be interpreted as the photon, since it actually couples to the electro- 
magnetic current 



k 

= Jl + \ J l ( 2 - 19 ) 
The charge of the electron e can be expressed as the harmonic sum of the electroweak couplings 

- = - + 4 (2.20) 

egg 

and eQk in Ij2.19|l is the electric charge in units of the electron charge. The neutral current 
coupled to the reads 

J l = E $>> > (9v ~ 9a 75 ) ^- (2.21) 

A; 

5 ^ = T 3 fc - 2 g fe sin 2 Bw, g\ = T 3 fc (2.22) 

where values of the generators for different particles can be extracted from table (2.1). Note 
that equation (|2.19(l assigns the electrical charge to the known particles from knowledge of their 
weak isospin and weak hypercharge. In terms of Noether charges, it leads to the well-known 
GeU-Mann-Nishijima relation 

Q = T 3 + ^Y (2.23) 
Gathering it all, the Standard Model Lagrangian can be expressed in the concise form 

£<SM = C-Kin + £-D + -CjVC + £-CC (2.24) 

with terms given in (|2.1()(l . I|2.8|l . Ij2.16|l and (|2.13|l . respectively. We have not yet talked about 
mass in the Standard Model. As it stands in l|2.24(l , a Yang-Mills theory does not allow a mass 
term for the gauge bosons, otherwise the gauge symmetry would be broken. This is rather 
annoying, for we know that all gauge bosons have mass with the exception of the photon, 
which remains massless. A way out is provided by the Higgs-Kibble mechanism 2 . 

2 See, e.g., 0. 
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2.1.1 Spontaneous Symmetry Breaking and the Higgs-Kibble mech- 
anism 

Generation of masses in the Standard model is driven by the spontaneous breaking of the gauge 
group SU(2)l x U(1)y, following the pattern 

SU(2) L xU(l) Y ^5 U(l) em (2.25) 

Symmetries in quantum field theory can be realized in two ways, depending on whether the 
Noether charges associated to the symmetries annihilate the vacuum or not. The latter possi- 
bility is termed spontaneous symmetry breaking. Formally, the simplest way of parametrising 
our ignorance about the vacuum structure is to postulate the existence of a scalar particle, 
called the Higgs boson, arranged as a complex isodoublct 

* = ( ^ ) (2-26) 
and enlarge the Standard Model Lagrangian with an extra piece 

£ H = - {D^D^ - fi 2 + A ($ f $) 2 (2.27) 

The wrong sign for the mass term shifts the minimum away from $ = and makes the field 
acquire a vacuum expectation value. Actually, there is a multiplicity of degenerate minima. 
In order to apply perturbation theory consistently, we have to single out one of the infinitely 
many degenerate minima. Typically, 



<*>o=(JJ , V = ] /Zf (2.28) 

The field <fi has to be shifted accordingly for perturbation theory to be valid. This can be done 
as follows 



v + H(x) 
~^2— eXP 



■ ' Xk(x) 

i 

2 v 



(2.29) 



Due to gauge freedom, the \ fields can be removed from the theory. Their degrees of freedom 
are then converted to longitudinal modes of the gauge bosons, i.e., mass terms. Plugging the 
expression above to the Higgs Lagrangian H2.27J1 we get 



-i 2 



= - {d^d^+ y -{v + Hf[W+W»- + n _ 2a 

-fi 2 K " ' " ' - A v " ' " ' (2.30) 



2 4 v \ " 2cos 2 <V 

(v + H) 2 .(v + H) 4 

— A 

2 4 

from which the masses of the gauge bosons can be readily determined to be 

m 2 H = -2/i 2 = 2Xv 2 (2.31) 



-^-(SOGeV) 2 (2.32) 
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u 2 



-£ [g 2 +9 2 j ~ (90GcV) 2 (2.33) 

This means that the existence of a non-trivial vacuum, to which the , W~ and Z® gauge 
bosons are not transparent, makes them acquire a mass, whereas the photon remains unaffected 
by the vacuum and thus massless. The remaining terms in (|2.3()(l arc the couplings of the Higgs 
boson to itself and to the gauge bosons through triple and quartic vertices. 

Explicit fcrmion mass terms in the Standard Model Lagrangian are also forbidden by gauge 
symmetry. Fortunately, the Higgs mechanism allows fermion masses to be accomodated into 
the Lagrangian as the Yukawa couplings to the Higgs field. The picture, beautiful as it is from 
a theoretical point of view, is nonetheless rather disappointing on the phenomcnological side, 
because the only thing we know experimentally is the vacuum expectation value, 

i_ 

v = (V2 G F ) 2 ~ 246 GeV (2.34) 

from which the gauge boson masses can be determined. However, neither \x nor A arc known, 
which means that there is no -prediction for the Higgs mass. Hopefully, the LHC will soon 
dilucidate if vacuum effects can be cast in the form of a Higgs field or more complicated 
structures are needed. 



2.1.2 Quark Mixing and CP Violation 

When introducing the quark doublets in (|2.4|) we denoted the lower quarks with a tilde without 
further explanation. The reason for this is that quark fields arc not simultaneous cigenstates 
of the strong and electroweak interactions. This is because the electroweak interactions do not 
respect the global SU(3)f flavour symmetry, while strong interactions do, and some ambiguity 
arises as to the definition of quark fields. Therefore, they can in principle couple as mixtures 
of these strong eigenstates. It is customary to shift this mixing to the lower quarks in each 
multiplet, something which can always be done with a proper redefinition of the quark fields. 
Therefore, the electroweak quark fields are rotated with respect to the strong quark fields as 

(2.35) 

where Vckm is a unitary matrix known as the Cabbibo-Kobayashi-Maskawa matrix, which can 
be parameterised with three Euler angles 6±, 62, O3 and a phase 5 as follows 

Vckm = 





cos 6 


'l 


sin t 


'1 


— sin 




cos ( 


h 




















1 


1 





)( 








e i6 









(2.36) 
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Due to the existence of a relative phase between quark fields, the Vckm matrix picks an 
imaginary part, which is the source of CP violation in the Standard Model 3 . This is unlike 
what happens if there only existed two families of quarks. The Euler angles are then reduced 
to just one angle, the well-known Cabbibo angle 6c- This mixing angle is responsible for the 
suppression of the bothersome flavour-changing neutral currents (FCNC) through the Gell- 
Mann-Illiopoulos-Maiani (GIM) mechanism 4 . The most important aspect of the two-family 
model is the fact that the imaginary phase in (|2.36(1 becomes global and can be factored out 
to render a self-adjoint matrix, thus forbidding any violation of CP. Therefore, the presence of 
the third family has a pivotal role in the origin of CP violation in the Standard Model. The 
CKM matrix structure (|2.36|1 can be simplified recalling that 62 and #3 are small angles, and 
to a very good approximation^ 



IK. 



\V ub 



\Vct\, s 



(2.37) 



It is customary, however, to employ the so-called Wolfenstein parameterization, which is an 
expansion in the parameter A = | V us \ 



( 



V c 



KM 



1 



2 



A AX 3 (p-i7 1 ) \ 



-A 



1 - 
± 2 



\ A\ 3 (l-p-ir)) -AX 2 



AX 2 



1 



0{X 4 



(2.38) 



/ 



Unitarity of the Cabbibo-Kobayashi-Maskawa matrix results in a set of unitarity relations that 
bind the CKM matrix elements, e.g., 



v ud v: b + v cd v; b + v td v; b = o 



(2.39) 



The above relation can be visualized geometrically in terms of the so-called unitarity trian- 
gle 5 of figure (|2.1|l . The amount of CP violation is then proportional to the area of the triangle. 

3 There is also another source of CP violation, the strong CP term in the QCD Lagrangian, to be discussed 
in the next section. Its impact is nonetheless negligible, and quark mixing bears the bulk of CP violation in 
the Standard Model. 

4 At the time of its introduction, only the three light quarks had been postulated. The GIM mechanism 
spurred the introduction of the fourth (charm) quark, with which the orthogonality of the Cabbibo matrix 
rules out any flavour-changing neutral current transition. The addition of a third family does not spoil this 
statement. 

5 Out of the number of possible triangles, the one coming from 12,^91 is the only one where all sides are of 
the same order in A. This is why people usually refer to it as the unitarity triangle. 



2.2 The Strong Sector 



16 



The strategy in phenomenology of K and B physics is to overconstrain the triangle to check 
the unitarity of the CKM matrix inside the Standard Model. Should inconsistencies be found, 
this would signal at physics beyond the Standard Model. 

After this aside to discuss quark mixing in the Standard Model, we complete our description 
of the Standard Model with the strong interactions. 

2.2 The Strong Sector 

So far, we have focussed on the Electroweak Interactions. Strictly speaking, one cannot speak 
of a unification of the Weak Interactions and the Electromagnetic Interactions, in the sense 
that we do not provide a single coupling constant. However, the picture introduced previously 
showed a deep intertwining between the two through the Higgs-Kibble mechanism. Likewise, 
the Strong Interactions could be seen as another patch of the Standard Model, since it comes 
with its own coupling constant. However, there is a deep non-trivial consistency relation 
between the three interactions in the cancellation of axial anomalies. Anomalies show up as 
a manifestation of a classical symmetry which does not survive quantization. In particular, 
anomalous triangle diagrams involving electroweak gauge bosons could spoil renormalization 
of the Standard Model. This potentially dangerous contributions are cancelled due to a subtle 
combined action of quarks and leptons in each of the three generations, provided that quarks 
appear with three different colours, Nq = 3. It is this inner consistency and self-dependence 
that make sense of the word Standard Model. 

Strong interactions 6 bind quarks inside the atomic nuclei due to the mediation of gluons, 
the gauge bosons of the color group. The interaction is described by a Yang-Mills quantum field 
theory under the non-abelian gauge group SU(3)c of colour. The Lagrangian of the theory 
reads, 

£qcd = -\ G# G% + i £ - "*« <*j ) M + ^2 ^ (2.40) 

q— 1 

where the second term is the Dirac Lagrangian for the six quark fields, and is the covariant 
derivative for the group SU(3)c, defined as 

{D lt ))=S i j a iA + ig a A ll (2.41) 

where is the gluon field, 

A* = E^r4 a) ( 2 - 42 ) 

a 

and A a the Gell-Mann matrices (see appendix A). The first term in (|2.40() is the gluon kine- 
matical term written in terms of the gluon strength field tensor, to wit 

G$ = d M Ai a) ~ &4 0) ~ 9s fate 4 6) 4 C) (2-43) 

where g s is the strong coupling constant and f a bc the SU(3) structure constants, whose expres- 
sions are listed in appendix A. From a mathematical point of view, gluons belong to the adjoint 

6 There is a vast literature on the subject. We refer to [21 El f° r further details. 
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representation of SU(3), while quarks sit in the fundamental representation. From a physical 
viewpoint, this means that quarks come as triplets of colour, whereas gluons are combinations 
of colour and anticolour and their number is fixed by the group structure to be N% ~ 1. Finally, 
the third term in (|2.4U[) is a source of CP violation allowed by symmetry arguments, where the 
dual field strength tensor G^ is defined as 

G> = i e^ pX G pX , £ 0123 = +1 (2.44) 

The parameter 6 is related to the neutron anomalous magnetic moment and happens to be 
extremely small. Present bounds 7 give 9 < 1CP 10 . This smallness constitutes the so-called 
strong CP problem, which still remains an open issue in particle physics. 

Contrary to what happens in QED, the strong coupling constant is small at high energies. 
In other words, at high energies quarks are loosely bound and behave as free particles, thus 
allowing perturbation theory techniques. This crucial observation, coined asymptotic freedom, 
has been awarded last year's Nobel Prize in Physics, and it explains, in particular, the successes 
of the parton model. But at low energies, quarks become more and more tightly bound, i.e., 
the strong coupling constants increases its value and we enter a non-perturbative regime. This 
is known as infrared slavery and leads to the notion of confinement. Due to its non-abelian 
nature, gluons not only interact with the quark families, but they also interact with each other, 
leaving a very intrincate picture of the strong interactions, which is far from being understood 
in detail. Confinement makes the quark-gluon picture transform to the hadronic picture we 
observe in particle accelerators. This leads to a misfortunate situation: we have been able to 
write down the QCD Lagrangian in terms of quarks and gluons but we do not know how to 
evolve it to its dual description in terms of the asymptotic hadronic states. Fortunately, effective 
field theories and symmetry arguments provide a way to progress in that direction. This goes 
under the name of Chiral Perturbation Theory and will be the subject of the remaining parts 
of the chapter. 



2.3 Chiral Perturbation Theory 

2.3.1 Effective Field Theories of the Standard Model 

Consider a toy QFT in which, for simplicity, the matter content is restricted to one heavy field 
"J and one light field ijj. For concreteness, its Lagrangian would be 8 



(2.45) 



At energy scales below the heavy mass threshold, only -0 is dynamical, while is frozen. At 
the level of the Lagrangian, this can be implemented by integrating out the heavy degree of 
freedom in the path integral sense. 



Z 



Dip r>* cxp 



i dx 4 C fadrf; 



Dip exp 



i / dx 4 C eff (tp, d^ip) 



7 See, for instance, 0. 

8 For general reviews on the subject, we refer to 171- 1121 . 



(2.46) 
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The generating functional H2.46J1 is obviously left invariant, but this is not so for the Lagrangian. 
Formally, what one has is 



which can be thought of as a series expansion in inverse powers of the heavy mass. All the 
information of the heavy degrees of freedom is encapsulated in the new coupling constants 
Cfc(m*), whereas the operators are made up with the (dynamical) light field alone 9 . Whereas 
symmetry requirements suffice to constrain the form of the (infinitely many) operators, the low- 
energy couplings are to be determined through a matching condition to the original theory. 
This ensures that physics is the same in the original and the effective Lagrangian as long as 
one stays below the heavy mass threshold. The advantage of working with an effective field 
theory is that the integration of non-dynamical degrees of freedom cleans calculational efforts 
from unnecessary complications. 

This quite simple idea can be shown to be consistent once it is provided with a criteria 
to organize the (infinite) operators that come out in Q2.47[l . This is what is known as a 
power counting rule. Predictability is then ensured even though the theory is intrinsically 
non-renormalizablc: with a proper renormalization procedure, all ultraviolet divergences can 
be absorbed in the couplings of higher dimensional operators in a controlled manner. In an 
effective field theory, renormalizability is ensured order by order. 

In the case of QCD, the mass gap between the lowest- lying pion octet and the first reso- 
nances of the spectrum (p(770), ai(1260), ...) induced by chiral symmetry is of the outmost 
importance, since it allows an effective field theory treatment. There is, however, an additional 
sublety to the example shown above: at low energies the relevant degrees of freedom change 
due to confinement, and one is forced to change language and speak in terms of mesons instead 
of the fundamental quark and gluon fields. Knowing the relation between the two sets of vari- 
ables would be tantamount to solving QCD: we would know how quarks and gluons assemble 
to form hadrons, something which seems quite out of reach by now. Fortunately, symmetry 
requirements alone provide a great deal of information. 

2.3.2 Chiral Symmetry 

Out of the six quark flavours in nature, at energies below the charm mass threshold only 
three of them are dynamically active, namely u,d,s, while the remaining decouple 10 . This is 
motivated by the separation of scales in the quark mass spectrum. Therefore, at low energies, 
only the light quarks are the relevant degrees of freedom. Using the projector basis of Ij2.3|l we 
can split the Dirac term of the QCD Lagrangian as follows 




(2.47) 



k 



(2.48) 



and rewrite the QCD Lagrangian in the following fashion 



Cqcd = - 1 G$ G(a) + 1 £ 9L 7^ q L + q R -fD„ q R - m q (q R q L + q L q R ) (2.49) 



g=i L 



9 See, for instance, 1131 . For a more general discussion, we refer to 1141 . 
°For more comprehensive accounts of chiral symmetry, we refer to 1151 - 1181 . 
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The first terms can be shown to enjoy a global U(3) R x U(3) L flavour symmetry, since left- 
handedness and right-handedness do not mix. This symmetry is explicitly broken by the mass 
term, which connects right and left chiralitics. However, since the light quark masses are small, 
this breaking is very soft and one can treat the mass term as a perturbation. Therefore, to 
a very good approximation, the QCD Lagrangian for three flavours is chiral invariant. This 
U(3) R x C/(3){ can be decomposed into the subgroup factors U(l)v x SU(3) R x 5f7(3)£ times 
an axial coset U{\)a- 

According to Noether's theorem, there should exist 8 conserved currents for each hand- 
edness, Ql and Qr, together with Qy and Qa- Qv expresses nothing but baryon number 
conservation and Qa was long ago shown to be broken at the quantum level due to the chiral 
anomaly. We are thus left with the 16 generators of the so-called chiral group SU(3) f R xSU(3){. 
Continuous symmetries, however, can be realized in two ways, depending on its response to 
the presence of the vacuum. We already commented on that when discussing the electroweak 
sector. If the vacuum is symmetric, then the currents annihilate it 

Q,|0 > = (2.50) 

and, following Coleman's theorem, there should exist degenerate parity multiplets in the spec- 
trum. On the other hand, if the operator does not annihilate the vacuum 

Qi\0 > + (2.51) 

the symmetry is said to be spontaneously broken. The multiplets are no longer degenerate 
and, according to Goldstone's theorem, a set of massless, spinless particles, as many as broken 
symmetries, appear in the spectrum. Furthermore, their parity and internal quantum numbers 
are the same as those of the broken generators. 

We happen to live in an asymmetric vacuum, as pointed at by the experimental fact that 
there are no parity multiplets in the QCD spectrum: axial-vector multiplets have higher ener- 
gies than their vector partners. What remains to dilucidate is the precise pattern of symmetry 
breaking. Starting from the chiral group SU(3)r x SU(3)l, Vafa and Witten showed that the 
subgroup SU(3)v could not be spontaneously broken by the QCD vaccum |19| . The remain- 
ing coset, SU(3)a, should therefore be spontaneously broken, so that the accepted pattern of 
symmetry breaking in QCD is 

G = SU(3) R x SU(3) L ^5 SU(3) V = H (2.52) 

The pion 11 multiplet is interpreted as the 8 Goldstone bosons of the theory, one for every 
SU(3)a generator. This explains why pions are lighter than they should be according to their 
quark content: chiral symmetry protects them from acquiring mass. 
An immediate consequence of Ij2.51|l is that 

(0|[Q,O]|0) f (2.53) 

where O is an operator which does not commute with Q. This quantity, signalling at spon- 
taneous symmetry breaking, is commonly called an order parameter. Clearly, according to 

11 Under pion multiplet we are actually referring to the full multiplet, which consists of the it states but also 
of the K and r\ particles. 
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our previous discussion, O must be a pseudoscalar operator. The lowest dimension one is 
O a = (775 A a q, which yields 

(Q\[Q%O b ]\Q) = - 1 - (0 | A 6 } g | 0)^^(01^10) (2.54) 

which means that the quark condensate is the natural order parameter of chiral symmetry 
breaking 12 . Our purpose is to find out the effective theory of the strong interactions in terms 
of the dynamics of the pseudo-Goldstone bosons. From a group point of view, they belong to 
the coset G/H . Let us parametrize them as <j> a under the representation 

W)=(^W,«bW) (2-55) 

It can be readily shown that £(0) satisfies the following transformation rule under a chiral 
rotation 

aW-^ffi&W hHha) > ZrW 9R rf{4>,g) (2.56) 



9={9l,9r) e G , 



h(<j>,g) G H 



For simplicity it is convenient to work with the combination U(4>) = £r(4>) The trans- 

formation rule is then independent of h, to yield 



U(<f>) 5 « U{4>) g{ 



(2.57) 



Obviously, physical results will not depend on the chosen representation. However, linear 
representations are known to lead to the appearance of extra particles. Since we only want to 
describe Goldstone boson dynamics, the only requirement we impose is that U be a non- linear 
realization. The conventional choice is to take £r(</0 = £z,(<W = an d 



U{4>) = u(& 2 = exp \jr Yl ^ A ° 



(2.1 



where the Goldstone bosons arc collected in a SU (3) multiplct as follows 

K+ \ 

X a 



a T2 



\ 



K° 



(2.59) 



Armed with all these definitions, we are ready to write down the Chiral Lagrangian, understood 
as the most general chiral-symmetric Lagrangian in terms of U and its derivatives. In terms of 

For an alternative view, see 
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the effective action, we have to solve for the following 



Z = 



Vq Vq VGfj, exp 



ax '-QCD 



VUexp 



I J dx^Ceffi^D^U) 



(2.60) 



This was the method employed by Gasser and Leutwyler [211 12"2"] to unveil the form of the 
chiral lagrangian order by order in the mass and momentum expansion. Since we eventually 
would like to compute QCD Green functions, we have to make sure that chiral Ward identities 
are fulfilled. This can be accomplished most easily with the inclusion of external sources in the 
QCD Lagrangian and afterwards imposing local gauge invariancc. The QCD Lagrangian with 
the addition of these auxiliary fields now looks like 



C 



QCD 



'-QCD 



Jif 



= ^-qcd + <ll t "( v ^+l 5a ^) < l- Qi s - H 5 p)q 



(2.61) 



where v^, a M , s and p stand for vector, axial- vector, scalar and pseudoscalar external sources, 
respectively, defined to be the following SU (3) elements, 



y — «i a) 



2 a » 



s = s - 



P = Po~ 



(a) 



(2.62) 

The QCD Green functions are then computed as derivatives of the generating functional with 
respect to the external fields around the point 

v (l = ap=p = 0, s = diag(m u ,TO d ,m s ) (2.63) 

Imposing local chiral symmetry induces the following transformation rules for the external fields 



q(x 

fn (x) = (v^ + a M ) (x 
(x) = - a M ) (x 
(s + ip) (x 
(s — ip) (x 



= (dR Qr + 9l Ql) (x) 

= (9r ?> g R + igR d^ R ) (x) 

= {9Ll^9 ] L +i9Ld^ L )(x) 

= {9R[s + ip]g\){x) 

= {9L[s-ip]g R )(x) 



(2.64) 



In principle, external sources are formal fields whose importance is to enforce Ward identities. 
However, they can also be thought of as background fields when one switches on the electroweak 
interactions. Then, the external sources l[2.62[l can be identified with physical fields via 



L = -eQAu- 



e 



sin Ow cos t>w 

e 



ip 



m u 






















m s 




—Q sin 2 9w 



Q 



Q { l ] wj- 



(2.65) 
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where Q is the matrix of the electric charge of quarks 

Q = diag(|-i,-l) (2.66) 

and Ql contains the flavour changing Kobayashi-Maskawa matrix elements, to wit 

V ud V us \ / 

Q ( l ] = ( I , Q<-> = I V: d I (2.67) 


2.3.3 The Chiral Expansion 

The effective action <|2.60ll can be expanded in powers of momenta and masses, once power 
counting rules arc provided for all parameters, 




U O(p ) 



s , p 0( P 2 ) (2.68) 

At lowest order, (D(p 2 ), we obtain 

C 2 = ^- (D^D^U + U^ X + X f U) (2.69) 

where the covariant derivatives are defined, as usual, by 

D„U = d^U-ir^U + iUl^ 

D„rf = dpU + iUirp-ilrrf (2.70) 

and the matrix \ is usually parametrized as 

X = 2B Q (s + ip) (2.71) 

Thus, at leading order in the chiral expansion we are left with only two low-energy constants, 
-Fo and Bq, which can be related to QCD parameters through a matching procedure, to wit 

F = -i ( 1 ^ 1 7r+(p) ) = U ~ 93.3 MeV (2.72) 



S o = ^2(0|^|0) = _-y<0|gg|0) , ( | qq \ )(2 GcV) ~ -[( 280 ± 30) MeV] 3 

(2.73) 

Fo can therefore be identified with the pion decay constant, and the parameter Bq is propor- 
tional to the quark condensate, which takes into account the effect of non-vanishing quark 
masses. This very first term was already written down by Weinberg |23| . and to lowest order 
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it reproduces the current algebra results. At 0(p 4 ), one has ten low energy couplings Li plus 
two contact terms Hi, to wit [221 



U = Li (D^D^U) 2 + L 2 (DfirfD v U) (D^D V U) + 

+ L s {D^D' 1 UD v U^D v U)+L i {D ll U^D^U) (U^ x + X f U) + 

+ L 5 (D^D^U (tf X + X^U)) + L e ( U^ X + X^U) 2 + 

+ L 7 {rf x - X ] Uf + L 8 { X ] U X ] U + tf X U ] x) - 

- i L 9 (FgDJJDjfl + F£ V D^D V U) + L 10 ( U*Fg'UF Lllv ) + 

+ Hx(F Rtlv F^ +F LfW Fn+H 2 (xh) (2.74) 

where the chiral field strengths are defined as follows, 

Fg v = d»l v -d u l» -i[l»J u }; F» v = 0V - 0V - * [r^r"] (2.75) 

As already emphasised, effective field theories have to be provided with a power counting rule. 
At (D(p d ), the diagrams that contribute are dictated by 

d = 2 + Y^ N n (n-2) + 2 N L , n = 2,4,6,... (2.76) 

n 

where N n is the number of vertices coming from 0(p n ) operators, and Nl is the number of 
loops. However, the drastic increase in the number of operators as one goes to higher orders in 
the chiral expansion, together with the poorly-known associated low-energy couplings, makes 
one stop already at 0(p 4 ) and only in certain observables to go to 0(p 6 ). 

As stated in our discussion of effective field theories, the low energy couplings Li should 
absorb the divergences of the theory at 0(p 4 ) to guarantee the renormalisability of the theory 
up to that order. Using dimensional regularization and the Gasser-Leutwyler renormalization 
scheme |2*2*] . 



Li = L r l (n) + -^ I --log 4tt + 7-1 + log// 



Hi = HT(f,) + ^ ^__i og47r + 7 _i + i og/ /] 

and evolution under the renormalization group is then given by 

r i , ( Mi 
16n 2 \fi 2 



LM = L r M + T ^log[^) (2.78) 



where I\ take the values listed below 



ri 


. 3 
32 


p _ 3 


r 9 = 


r 2 


. 3 
16 


p _ 11 
L 6 — 144 


Tio = 


r\ 


= 


r 7 = o 


fi = - 



(2.79) 
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Li 


Experimental value (10 3 • Li{M p )) 


source 


Li 


0.7 ± 0.3 


K e 4 , 7T7T — ► 7T7T 


L 2 


1.4 ± 0.3 


K e 4 7T7T — > 7T7T 




-3.5 ± 1.1 


K e 4 7T7T — ► 7T7T 


L4 


-0.3 ±0.5 


Zweig rule 


L 5 


1.4 ± 0.5 




L e 


-0.2 ±0.3 


Zweig rule 


L 7 


-0.4 ±0.15 


Gcll-Mann-Okubo, L5, Ls 


L s 


0.9 ±0.3 


M K o — M K + , L 5 


L 9 


6.78 ±0.15 


(r'Tv 


Lio 


-5.13 ±0.19 


it — > e^7 



Table 2.2: Values of the low-energy coupling constants appearing in £4 \2. 74\ ) YM\ \2(A \27j 
and references therein. See also \29^ . 



Since we ignore the details of the underlying theory (QCD), the matching procedure to de- 
termine the low-energy couplings cannot be performed in an analytical way, and one has to 
resort to approximations or experimental data, when available. Their values, as determined 
from experiment, are listed in table (2.2). So far, we have stressed that one key ingredient for 
building up chiral perturbation theory is the existence of a mass gap. However, we have not yet 
identified the expansion parameter, i.e., the scale A x which makes the expansion in momenta 
and masses 

(2.80) 



TO," 

-7— < 1 



meaningful. According to our discussion of effective field theories, the natural order of A x is 
that of the last integrated-out degrees of freedom, i.e., the first hadronic multiplet: p(770), 
ai(1260), ... Therefore, one expects 

(2.81) 



A x ~ 1 GeV 



This means that, even though the light quark masses are rather different 



m d 



0.5 



TO, 

— - 20 

m d 



convergence, at least for two flavours, is good 



my 
A v 



m d 
A v 



< 



m s 



(2.82) 



(2.83) 



Above that scale A x , new fields can become dynamical and have to be included in the form 
of new operators. Such an attempt was done in |24| . which proposed a Lagrangian made of 
the first resonance multiplets in the vector, axial, scalar and pscudoscalar channels. We will 
discuss this approach in detail in the next chapter. 

The Chiral Lagrangian terms of (j2.74(l have an additional even parity symmetry which QCD 
does not possess. This is so because in going from the first line to the second in H2.60J) attention 
has to be paid to anomalies, which manifest themselves in the integration measure jacobian. 
Once this is taken into account, one sees that the chiral anomaly is an 0{p A ) effect, known as 
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the Wess-Zumino-Witten anomalous term Following the conventions of 7;, it can be 

cast in the following fashion 



S[U,l,r] 



iN, 



wzw 



c 



240tt 2 
iN c 
~48tt 2 



da ijklm 
d 4 X 



, n i n j n k n l 



W(U, I, rY vali - W(l, I, rY uaf} 



where 



tf = Utd M U, 



(2.84) 



(2.85) 



The first line in ([2.84JI is only local in five dimensions and contains information about odd- 
parity processes purely among Goldstonc bosons. The second term contains the couplings to 
external sources that guarantee chiral invariance. The tensor W^ap in the second line reads 

W(U, I, r)^ = ( UlpUatfrf, + \ Ul^r v Ul a U^r p + %UdJ v l a U*rf, + id„r v Ul a tfrp 



1 

it 



(2.86) 



where the left-right exchange in the last line amounts to do the following substitutions 



U 



(2.87) 



2.3.4 Chiral Electroweak Lagrangians 

Having discussed in the previous section the low energy effective lagrangian of the strong 
interactions as an expansion in powers of the light quark masses and external momenta, we 
can now apply the same procedure to account for the electroweak interactions at low energies. 
We will consider them as perturbations to be included in the strong chiral lagrangian. 

At leading order, the electromagnetic corrections can be summarized in the following oper- 
ator 

C 2 = e 2 F*d 2 (Q L Q B ) +0(e 4 ;eV) (2.88) 

where 82 is the resulting low energy coupling, to be determined from experiment, and the 
following notations are adopted 

Ql = Q r = uQu^, Q = diagQ,-i -i) (2.89) 

As for the electroweak chiral lagrangian. we will be concerned with non-leptonic AS = 1 and 
AS = 2 processes, which are the relevant ones in the study of kaon physics. The procedure 
again relies on the method of external sources. Define a scalar source A to transform like 



A = 9l \9l 



(2.90) 
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Chiral operators linear in A are related to 0(Gf) processes, while chiral operators quadratic 
in A account for 0(G 2 F ) processes. Fixing the proper flavour content of A to be 



A = 



1 



A 6 - i A 7 



= A 



2:1 



(2.91) 



the nonlcptonic AS =1,2 chiral Lagrangians can be determined. At leading order, AS* = 1 
transitions are described by the following Lagrangian 13 



C 



AS=1 



G, 



A u F 



98 i-8 + -927 ^-27 + e ^0 9ew i-g 



h.c + 0{p 4 ;e 2 p 2 ) 



(2.92) 



,(o) 



(2) 



where Gf and Fq are factored out of the coupling constants gs, gn and g ew . C s , C s 
collect the operators which transform as (8l,1r), (8l,8r) and (27^,1^), respectively, under 
the chiral group. Superscripts stand for their power counting in powers of momenta. Thus, 



4 2) and are 0(p 2 ), while C ( 8 V) is O(p ). Their expressions are 



.(2 



(0) 



C 



(2) 



D»U) 23 



C 



c 



(2) 
27 

(0) 



(ux 32 c/ f Q) 



■ ,23 ( D»U)u + - ( D, L U ) 21 ( £W) 13 



(2.93) 



where the subscripts under the traces stand for the flavour content. 
As for AS* = 2 transitions, the effective Lagrangian is 



i -2 



A 32 (D^)UX 32 (D^)U 



0(p 4 



;eV) (2.94) 



where, again, Gf, Fq and the loop factor are kept explicit. A| =2 is the low energy coupling 
encoding the physics of the heavy degrees of freedom which have been integrated out from 
the Standard Model Lagrangian. Dealing, as we do, with effective field theories, this coupling 
plays the pivotal role in the description of the K° — K° mixing. We postpone its determination, 
which will be addressed in chapter 5. 

Contributions beyond leading order have also been computed for AS = 1 transitions. How- 
ever, similarly to what happens in the strong sector, the number of low energy couplings 
increases dramatically, to the point that the number of couplings exceeds the experimental 
skill. 



See, for instance, 1341 and 1351 and references therein to the original literature. 



Chapter 3 



Chiral Lagrangians in the 
Resonance Region 

3.1 Motivation 

In the previous chapter we have discussed chiral symmetry and its implementation to study 
the low energy dynamics of the pion octet fields. We estimated the radius of convergence of 
the chiral expansion to be A x ~ 1 GeV, which lies close to the first resonance multiplet. A 
natural extension would then be to further exploit chiral symmetry by incorporating these 
higher energy excitation states in a (chiral) Lagrangian, thus pushing A x above the 1 GeV 
threshold. 

One such Lagrangian would then allow a prediction of the low energy parameters of the 
strong interactions (Li defined in chapter 2) in terms of the masses and decay constants of the 
included resonances. Out of the infinitely many resonances one can include, phenomenological 
evidence seems to favour the lowest lying multiplets as having more specific weight in the 
Li than the remaining hadronic multiplets. If one splits the contribution from the hadronic 
spectrum to the Li as 

Li (//) = Lf + L t (ix) (3.1) 

where the first term corresponds to the contribution from the lowest lying multiplets and the 
remaining piece gathers the remaining contribution from the hadronic spectrum, this lowest 
meson dominance suggests that the first term basically saturates the previous equation. This 
was the philosophy behind one such construction of a resonance chiral Lagrangian, already 
proposed in the late 1980's jSl] 1 , which succeeded in reproducing the above-mentioned reso- 
nance saturation and yielded expressions for the Li in terms of a finite amount of resonance 
parameters. 

Later developments showed that that Lagrangian could actually be interpreted as an ap- 
proximation to large- Nc QCD |2Sll2Sj- This has many interesting implications. For instance, 
the predictions of the Li, once supplied with QCD short distance information, are constraining 
enough to provide a set of relations between the different parameters Li. To the degree that the 

1 There have been other attempts to describe the dynamics of the lowest energy resonances. See for instance 
1511551 . 
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resonance Lagrangian is an approximation to large- Ac QCD, those relations can be envisaged 
as manifestations of a hypothetical underlying symmetry of large- Ac QCD itself. However, for 
our purposes we will be mainly interested in the fact that once embedded in a large- Nc frame- 
work, the resonance Lagrangian inherits a power counting (that of large- Nc) which enables a 
consistent computation of quantum effects. This turns out to be important, because it allows 
one to check whether the alleged resonance saturation still holds at the quantum level or, on 
the contrary, the resonance Lagrangian needs more structure to reach that goal. 

Bearing this point in mind, it is convenient, before we proceed, to introduce the basics of 
the l/Nc expansion. 

3.2 The Expansion in the Number of Colours 

3.2.1 Motivation 

Based on ideas of Condensed Matter Physics, 't Hooft proposed to use the number of 
colours in QCD as an expansion parameter around Nc — * oo. The 1/Ac expansion is therefore 
a perturbative approach to QCD, not in the strong coupling constant but in a topological 
parameter, namely the index of the colour gauge group (actually, its inverse). As such, since 
the expansion parameter does not depend on energy scales, the approximation is valid over the 
whole range of energies. 

Somewhat countcrintuitivcly, increasing the number of colours leads to a much simpler 
description of the strong interactions. However, this is only true in some sense. We do not 
know how to solve for Large- Nc QCD, and presumably this is not easier than actually solving 
QCD itself, in the sense that there arc more dynamical variables in Large- Nc QCD than in real 
QCD. All we know is that the increase in number of colours allows a much simpler treatment of 
the gross features of the theory. In some sense, having Nc arbitrarily large sheds some light on 
generic features of the theory shadowed by the poor perspective of having Nc = 3. The natural 
question to ask is whether Large-Ac QCD has any resemblance to the real Nc = 3 world. The 
answer to this question relies heavily on the skill of the theory to pass the phenomenological 
test. Interestingly, a great deal of phenomena find a natural explanation in the large-Ac 
framework (sometimes it is the only explanation available). Just to mention some of them: 

• It provides a rationale to understand the Zwcig rule. 

• It is compatible with Regge phenomenology and, at least in two dimensions, reproduces 
Regge trajectories. 

• It favours two-body meson decays, as it is experimentally observed. 

• Under quite general hypothesis, it can be proven that chiral symmetry breaking follows 
the same pattern as that observed in real QCD |37j . 

3.2.2 Large-N c QCD 

The Lagrangian of the theory stands as in Nc = 3 QCD (after all, it is still a Yang-Mills 
Lagrangian) except for the fact that we have changed the gauge group 2 . That means that 

2 For general accounts of the 1/Nq expansion, we refer to 
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Figure 3.1: Feynman diagrams with its counterparts in the double line notation introduced by 't Hooft. 

each quark field, since they sit in the fundamental representation, appears as an Afc-plet. 
Gluons, on the contrary, live in the adjoint representation and enlarge their number to N c — 1. 
In practice, it is commonplace to approximate this to N c ; in other words, we are skipping 
the tracelessness constraint (we are taking U(Nc) instead of SU(Nc))- The missing gluon is 
numerically unimportant at sufficiently large Nc- Besides, it can be shown that the abclian 
factor is indeed suppressed at large Nc- Notice that the theory we are looking at differs from 
Nc = 3 QCD in that there exist far more gluons than quarks (the former scale with N c while 
the latter only with Nc). 

Our purpose will be to show the topological ordering of diagrams induced by the l&rge-Nc 
power counting scheme. For clarity, it is convenient to use 't Hooft double-line notation. 

3.2.3 The Double-Line Notation and Planarity 

Unlike the perturbativc approaches we are used to, in which we expand in powers of a coupling 
constant, we have to change our strategy. In perturbativc QED or QCD there is only one 
coupling constant which shows up to couple fcrmions and antifermions. That is why Feynman 
diagrams are so useful to organize calculations in powers of the coupling constant: you only 
need to count the number of vertices. 

In the l/Nc expansion, we need to keep track not only of the vertices of the theory (we will 
show later on that the coupling constant at large-iVc is colour-dependent) but also of colour 
flow inside the diagram. We would like to have a pictorical approach to be able to determine in 
an easy way the scaling of physical amplitudes with the colour factor. The double-line notation 
introduced by 't Hooft makes it transparent to extract the colour scaling of a given amplitude 
altering only slightly the Feynman diagrams we are used to. Consider the following pictorical 
recipee: represent every quark line by a colour line, right-faced arrow meaning colour flow, 
left-faced arrow meaning anticolour flow. This does not change much the Feynman picture. 
When it comes to gluons, however, we have to interpret their colour indices as a system of 
colour and anticolour line. Figure (3.1) shows some examples of converting Feynman diagrams 
to double- line diagrams. 

We shall next prove that the expansion in inverse powers of the number of colours is 
equivalent to a topological expansion, where the planar topologies are the leading ones. Let 
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us characterize a diagram by its vertices (V), propagators (P) and colour loops (L). Then 
every diagram is holomorphic to a polyhcdrum of L faces, V vertexs and P faces. Once this 
identifications are made, it can be shown that the equivalence classes of diagrams in 1/Nc 
power counting are holomorphic to the equivalence topological classes of polyhedra. Therefore, 
diagrams scale as 

M ~ ,g s - 2xE (g 2 N c f (3.2) 

where xe is a topological invariant called the Euler characteristic. That leaves the only non- 
trivial scaling that renders the theory finite to be 

9s ~ (JV )-i (3.3) 

which means that indeed the 1/Nc expansion is purely topological. This leads to a set of 
selection rules, which can be summarized as follows: 

• The leading diagrams are the subset of planar diagrams which minimize the internal 
quark loops and whose arbitrary number of internal gluon lines do not cross except at 
interacting vertices. Quark colour-flow lines have to be external. 

• Most commonly, a Green function is made of a product of currents, i.e., of quark bilinears 
of the form tpiT^ipj, where T is a Dirac basis element. Then, planar diagrams are those 
that contain a single quark loop sitting in the border. Non-planar diagrams are suppressed 
by N^ 2 . The introduction of an internal quark loop supplies a Nq 1 suppression factor. 

In other words, planar diagrams are those in which the colour factors compensate for the vertex 
factors. This explains why one can supply extra gluons in a planar diagram without altering its 
Nq scaling, while extra quarks cannot compensate vertices with combinatoric colour factors. 

3.2.4 Mesons in Large— N c QCD: selection rules 

In the previous section we have found that QCD with a large number of colours provides a 
power-counting scheme, in which the so-called planar diagrams are the leading ones in powers 
of Nc- However, as it happens with Nc = 3 QCD, nature forces us to deal with mesons instead 
of quarks and gluons. Therefore, we eventually want to translate the above statements into 
the meson picture of QCD. Surprisingly enough, this will provide us with far-reaching results 3 . 

First and foremost, we have to assume that large-TVc QCD is confining. This has to be 
imposed, otherwise it would be sterile to try to make any contact between QCD and large-TVc 
QCD. Confinement therefore is a hypothesis we make at the very beginning, not an output of 
large-TVc QCD. 

We will consider QCD two-point correlators as our starting point. In the following chapters 
we will need to consider three-point correlators and even four-point correlators. We will discuss 
its detailed large— Nc behaviour in due course. Consider 

IV(g 2 ) = z J Ae 1 '' 1 (0|T{^(x)j,(0)}|0) (3.4) 

where j M , j v are quark bilinears. Lorentz symmetry tells us that 



IW<? 2 ) = (q 2 g^ - q^U^q 2 ) 
We will follow the excellent review by Witten 1401 . 



(3.5) 
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Figure 3.2: A representative Feynman diagram contributing to the two-point function defined in \3-4\l 
together with its counterpart in double line notation. The cut illustrates the fact that in planar diagrams 
only one-particle intermediate states are allowed. 



Using the selection rules of the previous section, we know that planar diagrams will be those 
with no internal quark lines and arbitrary gluons. There arc infinitely many such diagrams, and 
we do not know yet how to resum them in an intelligent manner. However, once the diagrams 
are written in double-line notation, we can look for intermediate states by cutting the diagram 
(see figure (3.2)). Then one realizes that there is only one colour contraction, no matter how 
many gluons plague the diagram, all making up a single qq meson state. Summing all the planar 
diagrams would be equivalent to determine the large-Afc qq meson wave function in terms of 
quarks and gluons. This we do not know, but the crucial point to stress is that intermediate 
states are one-particle states, whatever its precise quark and gluon content. Multiparticle cuts 
are suppressed, i.e., they show up in non-planar diagrams, because internal quark lines are 
needed. 

Therefore, in full generality, 



n(<z 2 ) = E^ 

* — ' q — n 



(3.6) 



So far we have only discussed the analytical structure of the Green function, but indeed much 
information can be inferred for mesons themselves. The expression above has to match onto 
the parton model logarithm, and this requires the sum to be infinite. On the other hand, the 
right-hand side of (|3.6|) picks no imaginary part, so that the mesons are stable: their decay 
width is 1/Nc suppressed. 

This is as far as we can get with two-point functions: two-point functions provide infor- 
mation on the propagator of the mesons. To get further, one should consider three-point 
and four-point functions. Their large-iVc structure is depicted in figure (3.3). One can then 
show that indeed the singularity structure is restricted to poles: single poles, double poles, 
triple poles,... but no cuts. Additionaly, a detailed, though quite simple, study of three point 
correlators allow for a determination of the Nc scaling of meson decay constants, to wit 



(3.7) 



This is so because a three point function has the vertex of a meson decaying into two mesons. 
In an analogous manner, four-point correlators have information on meson-meson scattering. 
A similar analysis thereupon shows that meson scattering is 1/Nc suppressed altogether. 

We have ended up with a very nice picture of large-iVc QCD. Assuming confinement and 
planarity, the theory is made of infinitely many stable non-interacting mesons. This picture 
has some flavour of a semiclassical theory of hadrons, and indeed some efforts have been done 
towards going further in this very suggesting direction. 
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{jjjj) = J2 




Figure 3. 3: Schematic decomposition of two-point, three-point and four-point Green functions in the 
large-Nc limit. For instance, the last line shows the four-meson vertex, three-meson decay and two 
meson scattering in the t and s channels. 

3.2.5 Two-dimensional Large— Nc QCD: the 't Hooft model 

When introducing the foundations of large-Afc QCD, we have stressed the fact that even 
though there are some simplifications, these are quite akin to many-body systems, where gross 
features can be dealt with while the dynamical details are in general much more complicated. 

There is, however, an exception, and that is QCD in two dimensions, in short QCD( 2( j). 't 
Hooft showed that an analytical solution existed to solve the spectrum of resonances in 
terms of quark degrees of freedom in an explicit way. Also, unlike large-TVc QCD( 4d ), in the 
't Hooft model confinement is no longer an assumption but an explicit feature of the theory. 
This is not that surprising if one notices that in 2 dimensions the gluon propagator behaves in 
position space as 

V(r) ~ r (3.8) 

which indeed is confining. We will not delve here into technicalities and instead give a brief 
account of the most relevant results. 

The reason for the solvability of the model in the planar limit is the fact that the reduction 
in dimensionality allows one to choose a gauge in which gluonic self-interacting vertices simply 
cancel. One is then left with two diagram topologies, rainbow diagrams and ladder diagrams, 
which can be resummed to yield an eigenvalue equation, the so-called 't Hooft equation |41j . 
which in the notation of 021 reads 4 : 

Ml 4> n {x) = f q cfnix) -ml { dy (3.9) 
x(l-x) J \x-yy 

where x is the momentum fraction carried by the quark q inside the hadron. This equation 
shows that the resulting spectrum is discreet, labeled by index n, with masses piled up in a 
tower with characteristic mass parameter 



ml = 

4 See also E5I . 



9 2 N C 

7T 



(3.10) 
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This tower is asymptotically linear, i.e., at large excitation number n, masses behave as 

~ TT 2 m^n (3-11) 

and hadronic wave functions take the rather simple expression 

4>n{x) ~ V2 sin (nirx) (3-12) 

Additionally, one can show that the operator product expansion of two-point functions is an 
asymptotic series, as signalled by the factorial scaling of the Wilson coefficients, which is of the 
form g3j 

c 4 „ ~ (-I)"" 1 (^) (2n-l)I (3.13) 

We do not know how this whole picture changes when one is considering QCDf^j. Certainly 
it may very well be that none of them survive in 4 dimensions, but at least some of them 
seem to be independent of dimensionality. For instance, Regge phenomenology hints at an 
asymptotic mass scaling of the form of |3.1ip . Also, several considerations tend to catalog the 
OPE in QCD as an asymptotic series, a feature that QCD( 2 d) shares. Therefore, at least some 
qualitative features seem to survive. 



3.2.6 Large- N c ^PT 

Being x?T an effective field theory of the strong interactions, it is also possible to implement 
there large- N c methods in what has been termed large-N c xPT |45l 146) . The mass difference 
between the octet and singlet elements of a multiplct is a 1 /Ac-suppressed effect, meaning that 
in the large- N c limit octet and singlet fields are degenerate and should be merged in a nonet 
[4T| . an element of the enlarged U(Nc)l x U(Nc)r chiral group. The procedure to follow in 
the chiral limit will be the introduction of an additional diagonal Gcll-Mann matrix 



A (o) 



(3.14) 



with the following algebra 



a 60 



0, 



mbO 



rahO 



The matrix collecting Goldstone bosons gets therefore modified 



(3.15) 



1 >° 



hi 



7T 

1 



V 



K- K° 

v d • ■ y ' 

(3.16) 

but the chiral Lagrangian stays the same. The power counting induced by the l/A c -expansion 
can be inferred once the A c dependence of the low energy parameters are specified, i.e., 



K+ 
K° 



(3.17) 
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since it is a decay constant and 

B 0(1) (3.18) 

since it is proportional to a mass term. A quite more involved analysis also leads to the 
following. 

L\, L/2, L$, L5, L%, Lg, Liq, Hi, H2 0(N C ) 

L 4 , L 6 , L 7 , 2Li-L 2 0(1) (3.19) 

a hierarchy which is in agreement with experimental data. It should be noted that the scaling 
above is specific of the U(Nc)l X U(Nc)r colour group. For large Nc, we have shown that 
this is basically equivalent to SU(Nc) l X SU (Nc) r, but in real Nc = 3 QCD, this is no longer 
an allowed approximation. The 771 particle decouples from the pion-kaon octet and becomes 
massive due to the chiral anomaly. L-j receives contributions from r/i, and this means that 
when comparing with experimental data, L-j has a special status, since its large- Nc version is 
ill-defined .48 ]. Leaving this subtlety aside, we can immediately verify that the scale of chiral 
symmetry breaking A x introduced in the previous chapter scales as 

A x ~0(N C ) (3.20) 

which means that chiral corrections are suppressed by 1/N C powers. 



3.3 The Minimal Hadronic Approximation 

We have seen in previous sections that large-Ac QCD is able to constrain the analytical form 
of QCD Green functions to be meromorphic functions. For a typical two-point correlator, 

U (Q 2 ) = Y, 2 " 2 (3-21) 

where the sum extends to infinity, in order to reproduce the parton model logarithm. If we 
knew the solution to \aigc-Nc QCD, then we could derive analytically the infinite number 
of poles and residues of the Green function. Since this is not at hand, equation l|3.21(l is of 
not much phcnomenological use. A much useful approximation, which has been coined the 
Minimal Hadronic Approximation [491 15(J| , consists in truncating the series of l|3.21[l to retain 
only a finite number of terms. The poles and residues are then to be determined by matching 
onto the OPE to get the right short distances and to xPT to reproduce the long distance 
behaviour. From a mathematical point of view, such interpolating functions come under the 
name of rationale approximants. 

For the subset of Green function which are order parameter of spontaneous chiral symme- 
try breaking (SxSB) such an approximation should work better than with conventional Green 
functions (see figure (3.4)). Order parameters of SxSB receive no contribution from pertur- 
bativc QCD and are therefore expected to converge more rapidly at high energies (sometimes 
they are termed superconvergent correlators). This suggests that a small number of resonances 
should be enough to reproduce the behaviour of the Green function. In other words, that we 
do not expect dramatic changes as more terms are considered in (|3.21|) . but a soft convergence. 
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0.5 1 1.5 2 2.5 3 3.5 




Figure 3.4: The MHA to the two-point < VV-AA> correlator as compared to experimental data from 
ALEPH and OPAL (left) and to a model to be defined in chapter 6 (solid blue curve). Comparison 
is done in the Minkowski (Req 2 > 0) and Euclidean (Req 2 < 0) half planes (left and right figures, 
respectively). Notice that, even though the spectral functions are drastically different, their analytic 
extensions to the Euclidean fit impressively. 



Another fact supporting this point is that the region of intermediate energies is rather narrow: 
%PT is valid at about 1 GeV and common lore sets the onset of the perturbative regime below 
3 GeV. Not much room is left in between, and it is not natural to think that a bump will 
show up in the Euclidean region. Classical vector meson dominance also relied on some of this 
points and proved to be phenomenologically successful. Therefore, the approximation consists 
in truncating Ij3.21|l to 

(3.22) 

where in the second equality we have rearranged the finite sum as a product of poles and 
zeros. Such an approach is called an approximation because one can in principle add more and 
more resonances in the equation above to improve. Each new resonance in (|3~2^|l has to be 
included in a way compatible with the short distance and long distance behaviour of the Green 
function. Therefore, there is some freedom in choosing the constraints coming from high and 
low energies. However, one has to ensure that at least the leading OPE constraint is fulfilled 5 . 
This imposes a lower bound for the number of resonances to be included, which can be cast in 
the form of a theorem [^01 

Z-V = -pope (3.23) 

where M is the number of zeros and V that of poles in the Green function, whereas pope is 
the leading fall-off power (l/(Q 2 ) POPE ) of the Operator Product Expansion. Since Af > 0, this 
leads to 

V > pope (3.24) 
and so fixes the minimal number of resonances to be considered. 



n(Q 2 ) = 



F 2 



/ j 9 9 



1 




q 2 + mf 


n 



5 Actually, this will turn out to be the crucial ingredient to obtain the right matching between short and long 
distances in nonlcptonic weak interactions (see chapter 5). 
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3.4 A chiral Lagrangian with Resonance fields 

We now turn to our initial goal of including resonances in a chiral Lagrangian. The build- 
ing of an effective Lagrangian out of Goldstone bosons and resonance excitation fields in a 
chiral-invariant way can be achieved once the resonance fields are embodied with a chiral rep- 
resentation. Consider R and R\ as the octet and singlet components of a resonance multiplet. 
Their transformation under the chiral group has to be of the form 

R->hRh\ Ri-tR.!, heSU(3) v (3.25) 

This suggests, by analogy with the Goldstone fields, to collect the multiplcts as 

V, v {x)=V^=\ P - -^P u + ^8 K*" ] (3.26) 

for the vector channel and similarly for the axial, scalar and pseudoscalar channels. The 
transformation rules 1)3. 25[) and chiral invariance fix the interactions between the Goldstone 
and Resonance states. 

However, in our discussion of effective field theories we already stressed the fact that they 
consist of an infinite number of operators. The rigorous way to proceed is to identify a power 
counting to systematically order them according to their relevance. Hence, in the chiral la- 
grangian we saw that dimensional power counting provides a way to truncate the chiral expan- 
sion consistently. Unfortunately, no such power counting argument exists for the present case. 
As a starting approximation, one could consider the following Lagrangian |24j 





P + 


K*+ 


P 




K*° 


K*~ 


K*° 





£ 2 L4(1 + +)] = ^-^{A^D 

C 2 [S(0 ++ )} = ( C <U {SiW") + C m, (SiX+) + c di S u + ^ S u (x+)} 

i 

C 2 [P(1- + )] = ]T [id mi {P iX -) +id mi Pu(x+)} (3.27) 

i 

where 

u ll = iu^D IJ ,Uu^ f£ v = uF^ v V ±u t F^ v u, x± = ^x^ ± u^u (3.28) 

and the sum runs up to infinity to include each hadronic state. The previous Lagrangian is 
indeed chiral invariant, but only includes linear operators in the resonance fields. Without any 
power counting behind, we have to consider them as an ansatz and test it phenomcnologically. 
Integration of the whole tower of hadronic states would result in a determination of the low- 
energy couplings Li of the original chiral Lagrangian in terms of hadronic parameters (masses 
and decay constants). Dimensional analysis leads to a straightforward estimation [51j 

L ^J2jk (3 - 29) 
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For instance, 

L ™=^{^kr^k} (3 - 30) 

showing that the higher the resonance states the fewer the impact on the low energy couplings. 
This is supported by the phenomcnologically favoured vector meson dominance. A truncation 
of the sums in (|3.27|) to include the lowest lying multiplet for each channel thus seems to 
be favoured experimentally. This so-called lowest meson dominance can then be viewed as 
a natural extension of the vector meson dominance 5 . A step forward was later on provided 
by |25l I26|. which embedded the previous Lagrangian in a 1/Nc framework. As it stands, 
the Lagrangian l|3.27(l indeed contains an infinite number of narrow-width resonance states, 
as the large- Nc limit demands 7 . 1/Nc counting then yields the following behaviour for the 
parameters in l|3.27|> 

F v , F A , G v , c d , c d ~ O(0Vc) 

M v ,M A> M s ,M Sl ~ 0{N°) (3.31) 

Truncation of the sums then amounts to be working in the already mentioned minimal hadronic 
approximation. The importance of endowing the resonance chiral Lagrangian with a large- Nc 
framework is thus to give a rationale, namely the MHA, for the otherwise purely phenomeno- 
logical truncation of the infinite sums. In the following, it will prove convenient to adopt the 
MHA point of view. This eventually would allow us to compute quantum corrections (i.e., 
1/Nc corrections) with the Lagrangian (|3.27|) . However, we want to emphasise that large- Nc 
is able to provide a power-counting rule for quantum corrections once the Lagrangian is given, 
but it does not yield a power-counting criteria to build the Lagrangian 13.27f) . The issue of 
whether l|3.27|) gets any close to the large- Nc QCD Lagrangian or, on the contrary, fails, re- 
mains unanswered. All we know is that, at least for certain two-point Green functions, H3.27J) 
agreement with QCD is provided by imposing QCD short distances constraints, as we will see 
later on, but there are indications that for certain three-point Green functions this agreement 
ceases to hold. 

Enforcement of local chiral symmetry on the truncated version of (|3.27|) requires the defi- 
nition of a covariant derivative acting upon the resonance fields 

V„R i = d li R i + {r^r i ] (3.32) 

where the connection L M is defined as 

T M = ^ jw f [<9 M - ir M ] u + u[d^- il^) u f | (3.33) 

6 See, e.g., R>5| . 

7 Obviously, in the large-A^ limit the multiplcts have nine components, as we have seen previously, which 
means that the distinction between singlet couplings and octet couplings in the scalar and pseudoscalar sector 
is somewhat artificial, since one expects them to be one and the same. However, taking into account that the 
1/Nc expansion works, in general, worse precisely for these channels, it is phenomcnologically advisable to 
split them apart. If the I/Nq expansion happens to be a good approximation, this will show up in a nearly 
degenerate couplings. We therefore prefer to check if the splitting was unnecessary a posteriori, based on 
phenomenological grounds. 
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leading to the kinetic terms 
£ km (V,A,S,P) = 



E 

R=V,A 

E 



R=V,A 

E 

R=S,P 



1 M2 

-d^RJx^R^ + -^{RJ^R^ 



\ {V^RV^R-MlR 2 ) + \ (d^Rtd^ - M 2 Ri R\) 



(3.34) 



The use of a tensorial antisymmetric representation for the resonance multiplcts is especially 
convenient when dealing with gauge fields. A comprehensive treatment of tensorial represen- 
tation of vectorial fields is given in |24l I53| . For instance, the propagator is 



AP 



ic 



k 2 



9pp 



k v k\ 



k k 



(3.35) 



Other commonly used representations include the Yang-Mills representation (see j^lj and ref- 
erences therein) or the hidden symmetry representation |55j . Obviously some ambiguities arise 
between them, but they were shown to disappear once the right QCD behaviour is imposed on 
certain two-point Green functions |56| . A general proof of the equivalence between them was 
given in |57| , showing explicitly that they are linked through different redefinitions of the same 
Lagrangian. 

We can now integrate the newly introduced resonance fields in Ij3.27|l and (|3.34(l , something 
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which leads to a determination of the low energy couplings, to wit |24| . 



r 2 r 2 r 2 

T ., — I L d I '-d I 



L 2 =5^+ + + + + 

L *= S% + + WF S + + + 

L 4 = + 



L 5 = + +^+ + + 

(3.36) 

2-2 

L 6 = + - eM + 2MT + + 



L 7 = + + + 



dm d 2 
6MT 2MT 



2Mt 



L *- + + 2§l + 

L 9= l^f + + + + + 

L 10 = -§r + i^r + + + + 

With the previous determination of the low energy couplings the above-mentioned argument 
in favour of lowest-lying resonances can be tested. Experimental values as compared with 
predictions are summarized in table (|3.1J) . In the light of the results, there seems to be an 
amazing agreement, signalling at the fact that the lowest lying multiplets arc enough to account 
for the low energy couplings of QCD, a phenomenon coined thereafter resonance saturation. A 
more careful approach would be to reduce the phenomenological input and at the same time 
make the theory resemble QCD to a higher extent. For instance, one could demand the right 
ultraviolet behaviour of certain two-point Green functions, e.g., the pion electromagnetic form 
factor Fy (q 2 ), the axial form factor Ga (q 2 ) in 7r — *■ e f e 7 decay, the VV — A A two-point 



3.4 A chiral Lagrangian with Resonance Gelds 



40 



Li 


Experimental value ( (i = m p ) 


V 


A 


s 


Si 




Prediction 




0.7 ±0.3 


0.6 





-0.2 


0.2 


o 


0.6 


L 2 


1.3 ± 0.7 


1.2 














1.2 




-4.4 ± 2.5 


-3.6 


o 


0.6 





o 


-3.0 


Li 


-0.3 ± 0.5 








-0.5 


0.5 





0.0 


L 5 


1.4 ±0.5 








1.4 








1.4 


L 6 


-0.2 ±0.3 








-0.3 


0.3 





0.0 


L 7 


-0.4 ±0.15 














-0.3 


-0.3 


L 8 


0.9 ±0.3 








0.9 








0.9 


L 9 


6.9 ±0.2 


6.9 














6.9 


Lio 


-5.2 ±0.3 


-10.0 


4.0 











-6.0 



Table 3.1: Predictions for the low energy couplings Li from the expressions in \3.H(S\) versus 
experimental values. Masses and decay constants are estimated phenomenologically with the 
values given in 



function U LR (q 2 ) and the SS - PP two-point function II 5P (q 2 ), 

■2 



Fy(q 2 ) 
G A (q 2 ) 
Tl LR (q 2 ) 

n SP (q 2 ) 
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Fy i Gv i 



F 2 
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2Fy i Gy i — 



5 + E 



F' 2 



16B 2 



h 



Ml - q 2 



-E 
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E 



Mp. - q 2 



8 q 2 



(3.37) 
(3.38) 
(3.39) 
(3.40) 



Truncating the previous expressions to the first multiplet and comparing with QCD results in 
the following set of matching conditions |56| 8 



Fy Gy 
2Fy Gy 



Fy 



My 



together with jSS] 



(3.41) 
(3.42) 



FyM 2 



Fy - F 
F 2 A M- 



F{ 




(3.43) 
(3.44) 



8 Rccall that this is actually a big jump. We found a convincing argument to truncate the hadronic tower 
at low energies, namely, that higher mass resonances yield a (in principle) small contribution to the low energy 
couplings Li (sec 13.291 '!. However, what it is implicitly assumed now in the matching equations 13.4113.471 is 
that they also saturate the high energies. 
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and with 

Fy My — F\ M\ = -^"Ft (3.45) 

= (3.46) 

clMl-dlMl = 3 -^F* (3.47) 

where the last two equations are the akin Weinberg sum rules for the scalar-pseudoscalar sector. 
Combining l|3.41(l - (|3.43(l results in the following relations for the couplings in terms of Fq 

F v = 2G V = V2F A = V2F Q (3.48) 

The above relation, together with (|3.44|) . leads to the following expression for the masses 

1/2 



My = = AnF ( ^} 



where in the second equality use has been made of l|3.45(l . Therefore, we are able to express 
couplings and masses in terms of Fq alone. Additionally, this means that there appears a 
parameter- free prediction of the low energy couplings, to wit 

din = 3L 2 = - 8 -L 3 = 4L 5 = SL 8 = \l 9 = -L w = ^ = (3.50) 

where (|3.46|l and 1)3. 47|) were used in the prediction of the L5 and Lg couplings 9 . The set 
of relations found in the previous equation among the low energy couplings of the strong 
interactions, far from being a simple curiosity, seem to point at a more profound structure of 
the large--/V c limit of the strong interactions, as we commented on earlier. To illustrate this 
point clearlier, it is worth making an aside and consider an analogy with the determination of 
the p parameter in the Standard Model. 

The p parameter measures the ratio between the neutral and charged currents in the Stan- 
dard Model. Its expression at tree level reads 10 

1 M w (351) 
^cos^Mf {3 - 51) 

whose analog in our study would be the parameter Liq at tree level 

L w = —At (3.52) 

AM\ A M v V ; 

We already discussed in chapter 2 that the dynamical symmetry of the Standard Model is 
an SU(2)l x U(1)y- However, there also exists a global SU(2)l x SU(2)r symmetry, coined 
custodial symmetry 11 , which yields the prediction 

p=l (3.53) 



9 The analysis leading to the relations for L5 and Lg in 13.501 actually requires more input than just 13.461 
and l3"Tfl . We refer to H3 and H3 for details. 
10 See, for instance, 1591 . 

llr To be precise, it is SU(2)n which is referred to as custodial symmetry. 
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and the parameter is symmetry-protected. Much in the same fashion, enforcement of this 
would-be symmetry of the large- Nq limit of QCD leads to 



1 / 15 



4 V32tt 2 ^ 



(3.54) 



For simplicity, we will give a detailed accound of the previous statements for the p parameter 
in the linear sigma model. Consider the Lagrangian 

~o 2 { Tv T UfD » U ) "?E (Tt^D^u) +SSB (3.55) 



where U is the SU(2) Goldstone matrix, 



o=l,2 



7TT„ 



17 = cxp ( i-^ ) (3.56) 



and we have gauged U(l)y out of SU(2) x C/(l), so that the covariant derivative reads 



12 



D^U = 8^17 + 1^^11 (3.57) 



SSB in the last line stands for the spontaneous symmetry breaking terms, which are responsible 
for the vacuum expectation value /. The previous Lagrangian is therefore invariant under 
SU (2) global x U(l)i oca i. We can now enforce the right high-energy behaviour of certain Green 
functions, much in the same fashion as we did in (|3.37|l - H3.40ll . We consider the matrix element 
for pion scattering [SJj 



.s 



M(tt+tt- -> n 3 ir 3 ) = — [ I (3.58) 

If we require unitarity to be preserved, then p = 1, and as a by-product this leads to an extra 
custodial SU(2) symmetry of H3.55JI . We can turn the argument round by saying that the 
custodial symmetry ensures the condition p = 1 to hold. Recall that this is akin to what 
happens with the set of constraints H3.41JI - H3.47J) . which eventually lead to IJ3.54JI . 

If we now integrate out the a particle in IJ3.55JI . in a similar fashion as what we did earlier on 
with the resonances in IJ3.27JI and IJ3.34JI . we find the effective theory for the Goldstone bosons 

C =-^) i( T ^fiE (T^D^U) 2 , (3.59) 

which yields a prediction for the p parameter to be 

3 a' 2 m 2 

"<"> = 1 -4i!^7# (3 ' 60) 



2 See, e.g., chapter 2 of | 



3.4 A chiral Lagrangian with Resonance Gelds 



43 



P(fj) g £ io(m) 

g' <-> 1/N C 

<r <-> {V, A, S, Si} 
Custodial Symmetry SU(2) <-> Symmetry of QCD(A^ C — » oo) 

Table 3.2: Analogies between the p parameter in the linear sigma model and the parameter L±o 
in the Resonance Chiral Lagrangian, understood as an approximation to large-Nc QCD. 



where the first term is the prediction one gets from imposing the custodial SU (2) symmetry, 
while the second term comes from the integration of the a resonance. By the same token, 
integration of resonances in 13.27f) should lead to a prediction for Lio of the form 

where, if resonance saturation holds in QCD, A x would be a function of the resonance masses 
and couplings close to the scale of the integrated resonances, i.e., A x ~ m p ~ 1 GcV. For the 
sake of clarity, table 1)3.2(1 summarizes the parallelism between the determination of Lio in the 
resonance chiral Lagrangian and the p parameter in the linear sigma model. Another analogy 
can also be drawn with the well-known SU (5) grand unification group, where the seemingly 
unrelated couplings of the electromagnetic, weak and strong interactions (the analogs to the 
Li Gasser-Leutwyler couplings) happen to converge at sufficiently high energies, i.e., 

9su(3) c ( m gut) ~ gsu(2) w (Mgut) ~ 9u(i) Y (Mgut) (3.62) 
which provides a prediction for the electrowcak mixing angle Qw, to wit 

. 2 a / \ 3 55 a m GUT 

sin 6 W (ji) = - - — - log (3.63) 

a z4 7r /i 

The analogy would then be best illustrated as follows 



SU(5) 


sin 2 9 W = | 


a su(z) (Mgut) = (Xsu(2)(Mgut) = olu(X)(Mgut) 



Hypothetical Symmetry of Large —N c 


T 15 1 


6L 1 -3L 2 7 L 3 4 L 9 - Lio 8 jjt 8 ^ lfa2 



As stressed above, the fact that the lagrangian can be endowed with a large- N c behaviour has 
many advantages. First and foremost, it provides a consistent power counting rule, therefore 
allowing a consistent computation of quantum corrections. Obviously, the minimal hadronic 
approximation we are using throughout is not large-A'c QCD, but it can be systematically 
improved towards that limit with the addition of more resonance states. In the following 
we will concentrate on a prediction for Lio beyond leading order, as an example of how the 
adoption of the power counting supplied by large- Nc QCD enables one to treat the Lagrangian 
(|3.27|) at the quantum level. 
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3.5 Prediction for Liq beyond leading order 

Consider the following two-point Green function 13 

n^ R (q)S ab = 2i [ d A xe^ x (0\T (L£(z)J2£(0)t) |0) (3.64) 



where R£ and are the QCD currents 

K = q{x)r ^= q{x) , L% = q{x)r ^ q{x) (3.65) 
Lorentz symmetry constraints the tensorial structure to consist of just one form factor 

nf fl (<z 2 ) = Cs"V - <f<n n LR { q 2 ) (3.66) 

where the low energy expansion of U LR (q 2 ) yields 



F 2 

n M (Q 2 ) = - 7 |+4L 1 o + 0(0 2 ) (3.67) 
Q 



In other words, it provides a definition for the low energy coupling Liq, to wit 



(3.68) 



9 2 =o 



The previous result is just an example of the very general fact that low energy couplings of 
the strong interactions can be defined as coefficients of the Taylor expansion of QCD Green 
functions 14 . Our strategy hereafter would be to compute the contributions to I\ LR (q 2 ) both 
in the chiral lagrangian and in the resonance chiral lagrangian. The chiral lagrangian, up to 
0(p i ) and in the chiral limit reads 

+ L 2 {DyJJ^DyU) (D"U^D U U) + L 3 (D^U 1 D^UD^U 1 D U U) - 

- iL 9 (Fg'D ll UD v Ui + F^D^D U U)+L W (U^F£ U UF L ^) (3.69) 



13 We will follow closely P . 

14 A completely different thing happens for low energy couplings in the electroweak sector, where they are 
expressible in terms of integrals of Green functions (c/. chapter 5). 
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L 



10 



L u = -R^WW AAAA L v + R/j 



Lie 



Figure 3. 5: Matching condition for \3. 71\ at leading order, the left-hand side diagram coming from the 
chiral Lagrangian and the right-hand side ones from the resonance chiral Lagrangian. 

where we omitted the WZW term, which plays no role in the determination of two-point 
functions. By contrast, the resonance chiral lagrangian in the chiral limit (x ~> 0) reads 



C kin (V,A,S,P) = 



R=V,A 

i (v»sv„s- mis 2 ) + i (a»s,a„s, - mJX 
{wi {v ^ ) + t % {v ^ v) } + 



+ 



Fa 



2V2 

c d (Su^u**) + c d Si (u^) \ + 



F 2 



L w (U^F^UF Lflu ) 



(3.70) 



where Li and Li are as defined in (|3.1|) . The previous lagrangian lacks the pseudoscalar channel, 
since it does not couple cither to L^ nor to currents, as it can be checked by inspection. 
The same reasoning justifies the absence of vector and axial vector singlet states. Computation 
of H3.68(l at leading order is straightforward and leads to the matching condition 



Lw 



JL 

AM\ 



F 2 
AM 2 



Li 



15 



4 V327r 2 ^ 



L 



in 



(3.71) 



where l|3.48[) and 1|3.49[) were used in the second equality. This matching condition is depicted 
in figure (3.5). Next to leading order corrections are also relatively easy to compute in the chiral 
lagrangian side. We shall use the standard dimensional regularisation with the renormalization 
scheme chosen in the foundational papers of Gasser and Leutwyler. This leads to 



d 



? 2 n^ T (g 2 ) 



4£io(a0 



32tt 2 V 3 



log- 



(3.72) 



where the second term corresponds to the pion loop depicted in figure (|3.6() . 

Its counterpart with the resonance chiral lagrangian requires a bit more effort. Diagrams 
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to be taken into account are listed in figure l3~^l Their contributions render 15 

J \2 ° g (j, 2 ) + 2(4tt) 2 1," 
[GvY f 17 , Ml\ 3 1 (F A 



_d_ 

dq~ 2 



(4tt) 2 
3 1 



2 (4^) 2 V U 

Gy X 



17 . Ml 
30 n £ 



log 



M 2 



2 (4^) 2 V A 



1 M 2 
2- l0g ^T 



2 (4tt) 2 
10 1 



1 , M\ 

-3- bg ^ 



9 (4^) 2 V/J V6 
2 n 



log^ 



1 1 



2 (4tt) 2 



log^ 



1 



9(4^) 2 
4£io(m) 



log + 25 + 2B 2 - (2S 3 + 3S 2 ) log ^ 



M 2 

Vi 



32tt 2 



5 , 



(3.73) 



where is the mass ratio 



D 



M 2 

Vl 



(Ml 



Ml) 



(3.74) 



The first important point to stress is the vectorial dominance in the /3-function of the Lio 
coupling. Indeed, using H3.48J1 the logarithmic dependence on the axial channel identically 
vanishes. 



3 1 (F A y (1 Ml 
2(4^) 2 \ fj \2 ° S /x 2 



1 



2 (4tt) 



1 Ml 
3- l0g # 



5 1 



4 (4tt) 2 



(3.75) 



This again is intriguing: the relations between the parameters H3.50|l suffices to ensure this 
vector dominance. Also the scalar sector shows some degree of cancellation in the interplay 
between the singlet and the octet. It seems that somehow the hypothetical symmetry behind 
(|3.50() tends to protect the vector meson dominance even at the one-loop level. 

It is also worth mentioning that the previous equation is very insensitive to the 771 mass. 
Recall that the 771 mass arises as a 1/Nc effect, therefore decoupling from the octet away from 
the strict large- Nc limit. Should large- Nc be a good approximation to the Nc = 3 real world, 
one necessary condition is precisely this mild dependence on the r\\ particle. Therefore, we find 
this a very appealing feature, pointing at a smooth transition between the large- Nq limit and 
N c = 3 QCD. 

15 To be consistent, we will all along stick to the Gasser-Leutwyler renormalization scheme introduced in 
chapter 2. 
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R» WVV AAAAy L v 

7T 

Figure 3. 6: Next to leading order contribution to 13. 6&I) stemming from the Chiral Lagrangian and the 
Resonance Chiral Lagrangian. The pion loop is reproduced by both Lagrangians and identically cancels 
in the matching condition (see main text). 



Equating lj3.72|l and (|3.73|1 yields the determination for Lio up to next to leading order, 
15 



4 L r w (ji) 



32tt 2 \/6 
3 1 



2 (4tt) 2 
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'(4^)2 

3 1 
"2 (4tt) 2 
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"2 (4tt) 2 

4 1 



Fa 
U 

X 

1 
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2 (4tt) 2 
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2 (4tt) 2 
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(47T) 



1 M J 

log^ 



log 



10 1 
~9~(4tt) 



log^ 



9 (4tt) 2 
+4 (2io) r (M) 



M 



- + log — f + 2£ + 2B J - (2B J 



3B 2 )log|f 

x 171 



(3.76) 



We are now in a position to assess whether the resonance saturation survives at the quantum 
level, i.e., whether the contribution from the higher mass multiplcts to Liq, which we noted as 
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Lio, can be dropped from 

L 10 (M)=£fo(M)+£io(/i) (3.77) 

It turns out that one can not get rid of this term (see for details) or, in other words, the inte- 
gration of the resonances does not predict the right evolution for L\q under the renormalization 
group equation, as dictated by chiral perturbation theory 

1 A/f 

L 10 (/i) = L w (M p ) - — log (3.78) 

647T Z /! 

Our conclusion is that the Lagrangian of |24j . despite yielding successful phcnomcnological 
results at tree level, fails at the quantum level. We are not the first ones to claim that that 
Lagrangian is incomplete 16 . The message we wanted to convey in |62| is that one can go 
beyond tree level with such resonance Lagrangians, because the large- Nc framework provides 
an expansion parameter with which to do quantum corrections in a consistent way, namely 
1/Nc- 

In the past three years there have been considerable efforts to improve on the Lagrangian of 
|24j . by adding more operators and looking for agreement with QCD short distances of certain 
two and three-point Green functions at the quantum level in the 1/Nq expansion 17 . However, 
as already pointed out, the absence of a power counting rule to tell which terms have the 
bigger impact makes it difficult to make further progress in that direction 18 : the new operators 
included ad hoc to improve on certain Green functions might as well spoil other Green functions 
not yet considered. 

From this point of view, we think that we can provide one of the simplest starting tests 
any resonance Lagrangian has to pass. The addition of more operators has to render, among 
other things, the right running for Lio(p). This is a necessary condition to eventually yield a 
finite prediction for U LR (q 2 ), one of the simplest QCD two-point functions. Should we have 
that Green function under control, we could then move to more complicated ones. However, 
to the best of our knowledge, no prediction exists yet even for Lio(/z). 



16 See, e.g., 1531 1511551 . 

17 See, e.g., 1661 1671 and references therein. 

18 Rccall that we are actually attempting to model the large- Nq Lagrangian of QCD. Therefore, large-TV^ 
power counting rules can only tell us how to go to the quantum realm once the Lagrangian is given. For this 
ambitious task we would need to know much than we nowadays do about the l/Nq expansion itself. 



Chapter 4 

Hadronic Matrix Elements of 
Kaons 

In chapter 2 we saw that the CKM quark mixing matrix had a non-factorizable phase signalling 
at CP violation in the Standard Model. In this chapter we will show how this phenomenon 
manifests itself at the meson level in neutral particle-antiparticle systems, the paradigmatic one 
being the K — K° system, to which we will devote our attention 1 . Two mechanisms combine 
as sources of CP violation: mixing of K° — K° and kaon decay The first part of the chapter 
is oriented to characterize these effects in terms of a few phenomenological parameters, to be 
determined later on. In chapter 2 we already saw that the Standard Model at low energies 
admitted an expansion in powers of momenta with the Goldstone octet as dynamical fields. 
This will be the appropriate tool to be employed all through our analysis. 

4.1 Phenomenology of kaon CP- Violation 

The neutral kaons that appear in the chiral lagrangian are the ones predicted by the Eightfold 
Way of Cell-Mann and Ne'eman, namely 

K° = sd, K° = ds (4.1) 

They are eigenstates of strangeness, but they have no well-defined CP parity. They transform 
in the following way 2 

CP\K°) = -\K°) (4.2) 

In the absence of electroweak interactions, neutral kaons are stable and they behave as distinct 
states, since they are protected by their different strangeness quantum number, i.e., no mixing 
nor decay is to be expected, and since they are antiparticles, CPT guarantees that they have 
degenerate masses. However, once the electroweak interaction is switched on, they are allowed 

x We will follow the treatments given in 1141 and 1681 . 

2 Actually, in full generality they are related by a phase £k , 

CP \K°) = £ K \K°) 

which we choose to be £x = — 1 f° r convenience. Physical quantities should be independent of the phase choice. 
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to decay weakly into two-pion and three-pion states. Strangeness is no longer a conserved 
quantum number and, owing to the well-defined CP parity of the decay modes, it is convenient 
to build the CP conserving states, 

k "U r "- r °)' k "U r " +r °) (4 - 3) 

whence, 

CP\K 1 ) = +\K 1 ), OP\K 2 } = -\K 2 ) (4.4) 

Hence, were CP conserved in the neutral kaon system, K\ should decay exclusively into two- 
pion states and K 2 into three-pion states. The observed difference in lifetimes 

r^-lO 3 ^ 1 (4.5) 

is what one would expect by naive phase space arguments. However, these states are not the 
physical ones, because even CP is violated by the electroweak interactions. Quite generally, we 
can write the physical kaon basis as 



1 



p\K°)+q\K ' 



\K S ) = 1 (p\K°)-q\K )) (4.6) 

Departure from CP conservation is signalled by the parameter e, defined as 

- P-Q 
e = 

p + q 

in terms of which, the physical kaons look like 



(4.7) 



4.1.1 K° — K° mixing 

If we define the wavefunction describing the kaon system evolution as 

|*(t)) SCl (t)|A-°)+ C2 (t)|^ ) (4.9) 
elementary quantum mechanics of two-state systems renders the following Schrodinger equation 

i±Mt))=MMt)) (4.10) 



where Ai stands for the Hamiltonian 

M = 



K°\H SM \K°) (K°\Hsm\K ) 
K°\H SM \K°) (K°\Hsm\K°) 
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whose entries can be parametrized as follows 



The imaginary part in the equation above accounts for the fact that kaons actually decay. 
Conventional quantum mechanics perturbation theory then yields 

n 

which provides a link between kaon masses and kaon decay widths in terms of the dynamics 
of the Standard Model. CPT invariancc, which is a consequence of any quantum field theory, 
further constrains the mass matrix elements |50j . 

Mn=M 22 , M 12 = {M 21 )* (4.13) 

The diagonal terms My are just the masses of the neutral kaons, whereas the off-diagonal 
entries account for kaon mixing. In terms of kaon masses and decay widths the parameter e 
can then be written as 



p— q _ i ImMi2 - § ImT 12 



(4.14) 



p + q 2 ReM 12 - | Rcr 12 

This in particular means that if kaons were stable nonetheless they would mix. 3 This suggests 
to study kaon decay and kaon mixing separately. 

For kaon mixing, our main concern will be the kaon mass difference, to be defined as 

ArriK = tul — ms (4-15) 

A close look at equation l|4.12|) reveals that the kaon mass difference receives contributions 
from the AS = 2 sector of the Standard Model, which connects both neutral kaons through 
the well-known box diagrams, depicted in figure ^4.1[l . This contribution leads to the AS* = 2 
Hamiltonian we mentioned at the end of chapter 2, whose leading order in inverse powers of 
the masses of the integrated out particles is 



r 2 m 2 
r AS=2 _ <^F m W 

L2 ~ 4tt 2 

where 



Ag-Fi + X 2 t F 2 + 2X c X t F 3 



Ci{n) Qi(v) (4.16) 



A c = V cd V* s , Xt = VtdK (4.17) 

and 

Q = SL^dLSL^dL (4.18) 

Fi,F 2 ,F 3 account for the electroweak and strong loop corrections, to be defined in the next 
chapter. Conventionally, one defines the bag parameter Bk in the following manner 

(K°\ Q(0)\K°) = - z Flm\B K {i?) (4.19) 



There is increasing evidence that neutrinos indeed follow this pattern of mixing without decaying. 
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such that it becomes the parameter which governs kaon mixing. The determination of the Bk 
parameter is the subject of the next section. As it is clear from (|4.12|) . though, Bk alone does 
not provide the whole kaon mass difference. There are also long-distance contributions in the 
second term in (|4.12|1 . with which we will deal in the next chapter. Furthermore, the last term 
in H4.12f) is a genuinely long-distance contribution. If we restrict the sum over intermediate 
states to be saturated by the tttt exchange, which is a good approximation taking into account 
that it is this the dominant kaon decay mode, the long-distance contributions can be cast in 

pieces, as it appears in (|4.12() . There is some rationale in 



the form of 



(AS= 1) x (AS = 1) 



neglecting these 



(AS* = 1) x (AS = 1) 



has been performed so far. 



contributions but no estimation of these terms 



4.1.2 Kaon non-leptonic weak decays 



So far we have only remarked the fact that the electroweak interactions do not conserve CP 
parity, but we have not assessed yet to what extent. Were CP a good quantum number, then 
the aforementioned CP basis K\, K% would be the physically observed basis, K$, Kl> The 
experimental fact that the long-lived kaon indeed decays into a two-pion state rules out this 
possibility and it is seen as evidence for CP violation. To evaluate the size of this effect, several 
observables can be constructed, for instance 



V+- 



mo 



{ 7T 7T 



\T\K L ) 



\T\K S ) 



(7r°7r°|T|A's 



(4.20) 



It is however convenient, in view of the phenomenological treatments to come, to express any 
observable in terms of well-defined isospin states . Bose symmetry forbids 1=1 transitions, 
so that we are left with 1 = and 1 = 2 states, 



^l(- - )o> + ^l(- -°) 2 > 



We define 



((7nr) I \T\K°)=iA I e: 5 ' , 
The natural ratios one can build thereof are 

A\K L 



((7T7T)j\T\K°) 
' (7T7r)j =2 ] 



-i A} e iSl 



A[K S -> (7r7r) 7= o] 



(4.21) 
(4.22) 

(4.23) 



_ A[K S -» (7T7r)j =2 ] 

uj — 

A[K S -> (7T7r)/ = o] 



4 Actually, there is a deeper reason: CPT invariance combined with Watson's theorem relate 
kaons to the same isospin final state. For a proof, we refer to 1501 . 



(4.24) 

both neutral 
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_ A[K L -> (7T7r) /= o] 

Nonetheless, it is commonplace to use the following basis instead 

s K (4.26) 



u (4.27) 



C-sk-l;) (4.28) 



which relate to the physically measurable quantities rj^ and ?7oo as follows 

1 

V+- = £k + % — — — 
1 + 7T 

r/oo = e K -2e' K 1 (4.29) 

1 — VZld 



Plugging l|4.22|l in (|4.25|l we get 



and 



e Re^o + HmA 
£jf = t; — ; t- (4.30 



s , = i ( _ g2 W 2 ReA -- W ReA 2 (tfa _, o) 

The previous expressions can be simplified using the following approximations, supported by 
phenomenology 

I 2 < 1 , elmAo < Rc4 (4.32) 

which yield 

ImAo 

Ek ~ e + i ■ 



ReA 

1 Re^ 2 flmA 2 lmA 



'- K 



e^^-^+f) (4.33) 



V2 Rc^4o \RcA 2 ReA , 

Our expression for can be further simplified if we make the additional approximations 

Am K T s AT T s 
ReMl2 ~ ~ , Rer 12 ~ -f 

T 12 - -(ReAo + i W ) 2 (4.34) 

where in the last one we consider the width to be dominated by the tttt decay channel. Imple- 
menting the previous conditions onto the definition for e (|4.14(l results in 

1 ^.ImM^Wj (4 35) 



1 + i V Am K ReA c 
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from where, plugging the previous equation in (|4.33|) . we end up with our final expression for 

1 / ImA/12 lmA ' 



e K = -=[ + — e** (4.36) 

From the previous equations we can readily see that ek probes both kaon mixing and kaon 
decay, whereas e' K is sensitive to direct CP violation in the form of the interference of the 
different isospin decay modes. 



4.2 Effective Lagrangians for kaon processes 

Having discussed the basic formalism for the phenomenology of kaon physics, we turn our 
attention to its description in the framework of the Standard Model. In chapter 2 we introduced 
the low energy description of the Standard Model. We will focus on strangeness changing 
processes, AS = 1 and AS = 2, which account for kaon decay and kaon mixing, respectively. 



4.2.1 AS = 1 transitions 

AS = 1 are mediated in the Standard Model through the exchange of virtual W particles. 
Nevertheless, we can make use of effective field theory techniques to simplify our analysis. 
Operators in the flavour changing charged sector of the Standard Model giving rise to AS* = 1 
transitions appear as the convolution of two currents with a W. We can therefore integrate the 
W particle, which results in a series expansion, the first term of which (the lowest dimension 
operator) is equivalent to shrinking the propagator to a point-like vertex. We are then left 
with a product of two currents, and we recover Fermi theory of the weak interactions. The 
remaining terms of the expansion are higher-dimensional local operators which conform the 
Operator Product Expansion, a topic to be discussed further in the last chapter. We end up 
in this very first stage with a set of local operators whose Wilson coefficients contain inverse 
powers of the W-mass. 

We then have to run the operators down to lower energies until we find the next particle 
to be integrated out. The evolution of the operators and Wilson coefficients in the energy gap 
between two particles is dictated by the renormalization group analysis. Heavy quarks can be 
integrated following these steps: at a particle mass threshold one integrates out the particle 
form the effective action, and evolves down the resulting operators with renormalization group 
techniques. Proceeding this way one can end up with an effective Lagrangian for strangeness 
changing processes of the form 

10 

v i— 1 

where 

A„ - V ud V* s (4.38) 

are the CKM matrix elements defined in chapter 2 and Qi are the leading operators {i.e., the 
lowest dimensional ones) made up with the left-over light quark fields. 

The method sketched above is not free of subleties: renormalization evolution equations take 
into account the behaviour of the operators under the strong and electromagnetic interactions. 
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This not only affects the Wilson coefficients, but quite often generates new operators, which mix 
between each other under the rcnormalization group equation flow. So, unless the operators 
are multiplicatively renormalizable, new operators will appear 5 . That is why we started from 
one single operator in the lagrangian of the Standard Model and end up with ten, to wit 

Qi = 4 {s L Yd L ){uLl^u L ) , Q 2 = 4 {s L ^u L ){u Llll d L ) (4.39) 

which are the so-called current- current operators. Furthermore, there are the so-called penguin 
operators 

Qs = 4 (s L ^d L ) ^ («l7h9l) Qi = 4 ^ {s L ^q L ){q Llll d L ) 
9 9 

Qs = 4 (s L ^d L ) (QRl^QR) Qe = -8 J] (s L q R ){q R d L ) (4.40) 

9 9 

and also, with inclusion of the electromagnetic corrections, there appear the following elec- 
troweak penguin operators 

Qi = 6 {sLl^d L ) Y e 9 (QbI^Qb) Q& = 4 ^ e q (si7 A 'to)(9i?.7p d L) 

9 9 

Q9 = ^{sLl^d L )^2e q {q Rlll q R ) Q w = -12 ^ e q {s L ^q L ){q Llll d L ) (4.41) 

9 9 

Before turning to our low energy chiral description we have to identify the symmetries of the 
lagrangian, which will be our guiding line. The procedure we have sketched for integrating 
out particles breaks down as soon as one goes down to energies were perturbative analyses are 
no longer valid. We have to change our description in terms of quarks and gluons for one in 
terms of hadrons and connect them through symmetry arguments. As stated previously when 
discussing Chiral Lagrangians in the electroweak sector, the previous basis Ij4.39|) . H4.40J) and 
(|4.41|) transforms under the chiral group as (8l,1_r), (8l,8r) and as (27^,1^). In order to 
establish contact with the phenomcnological analysis we carried out in the previous section, it 
is also conventional to split these contributions in the isospin basis. Therefore, 

Q?} 3 = ({Q2 - Qi,Q 3 , Qi, Qs, Qe}) (4.42) 

groups the operators which transform as a left octuplet with AI =1/2 transitions. Not all of 
them are independent, actually they are subject to the constraint 

Q 2 - Qi = Qi - Q 3 (4.43) 

On the other hand, the following combination yields the 27-plet component 

Q< 27 > = 2Q 2 + 3Qi - Q 3 = | (Q$J + 5 Q$?) (4.44) 

where, contrary to what happens for the octuplet, both AI = 1/2 and AI = 3/2 are allowed. 
The expressions for them read 

(27) ~ 

Qi/2 = SLl^dhULl^UL + s L 'f Pl u L UL'-f fl d L + 2 s L ^d L d L '-f fl d L - 3 s L "f Pl d L s L j IJ ,s L (4.45) 

5 For an exhaustive review, see [I]. 
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,(27) 



(4.46) 



The low energy realization of the previous short distance four-quark operators can be achieved 
with the use of chiral symmetry requirements. This yields, to lowest order, 



r AS=l 



Gf 
'V2 



A« + .927^27 + e 2 .9 et0 (C/A 3 2[/ t Qi?,)| 



where 



(4.47) 
(4.48) 



C 8 = F* (D^D»U) 23 
groups the operators that transform under SU(3)l X SU(3)r as (8l, 1r), while 

£27 = F± [i^D,U) 23 (tfD»U) u + 2 - (U^D,U) 21 (^D,U) 13 ^j (4.49) 

gathers the operators that fall into the representation (27l, 1^). The subscripts in the previous 
equations mean that the expressions inside the brackets are to be multiplied by the spurious 
matrices (Sij is the Kronecker delta) 



^23 — $i2 Sj3 



An — o~n Sji , A21 — Si2 Sji 



A 



13 



<5ji 6 



j3 



(4.50) 



which project out the right flavour content. 

The last term in (|4.47(1 is the low-energy realization of the electromagnetic penguin operators 
Qr - Qw- where 

Q« = diag(|4-i) (4.51) 

and thus it transforms as (8l,8r). 

Experimental numbers on K — > tttt decays give 



|32 7 |~0.16 , 
whereas naive factorization would predict 



.98 



98 + 7.927 
5 



.927 



5.1 



(4.52) 



(4.53) 



The previous discrepancy between theory and experiment is usually coined the AI =1/2 rule, 
of which no convincing account has yet been given. The bulk of the enhancement needed 
seems to point at the long distance contributions of the Q\ — Q2 operator. A penguin-based 
explanation in terms of Qq is ruled out, unless 1/Nc corrections turn out to be really sizeable 6 . 



4.2.2 AS = 2 transitions 

AS = 2 transitions arise in the Standard Model from double virtual W exchange, the so-called 
box diagrams, depicted in figure 1)4.1(1 . After integration of heavy particles, the four-quark 

6 See, e.g., the discussion in 1681 . 
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Figure J h l: Box diagrams for AS = 2 transitions 
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effective Lagrangian can be written schematically as 



C 



AS=2 



(4tt 2 ) 



A c F\ + A t F 2 + A c At i 7 ^ 



Notice that only one operator appears, namely 

Q (/i) = SLjfidLS^dL 



(4.54) 



(4.55) 



Its Wilson coefficient has been unfolded in the so-called charm-charm, top-top and charm-top 
contributions . We have adopted the following conventions 



Ai — VidV* s , 



i = c, t, u 



(4.57) 



as a short hand notation for the CKM matrix elements. The Fi are functions of the integrated 
out particles down to the strange quark mass. They are usually cast in the form 



Fi = r\i Si 



(4.58) 



where the r\i function collects the strong quantum corrections and the Si are the Inami and 
Lim functions, which account for the electroweak leading order corrections 0]. As already 
commented on earlier, AS = 2 transitions are governed by the bag parameter Bk, defined as 



\K°\Q(n)\Rt 



^F 2 K m 2 K B K (fi) 



(4.59) 



Sometimes it is nonetheless more useful to work in terms of the invariant, i.e., scale-independent, 
bag parameter Bk, 

( K° | c s=2 {y) \K°} = ^ F 2 K m K B K (4.60) 

where cg = 2(/i) is the term in square brackets in l|4.54|l . In the so-called vacuum saturation 
hypothesis, using that 

(0\s L ^d L \K°(p)) = -iV2F K p» 



(4.61) 

7 This does not mean that the up contribution is missing. To get to 14.541 we have set m u = 0, and we used 
whenever necessary the unitarity bound 

X u + A c + A t = (4.56) 
of the CKM matrix elements to express the result solely in terms of charm and top parameters. 



4.2 Effective Lagrangians for kaon processes 



58 



one immediately concludes that 



B K 



1 



VS 



Different approximations suggest the following bounds, 

0.3 < B K < 0.8 
For instance, the large-Nc prediction yields 



B 



K 



N c 



(4.62) 



(4.63) 



(4.64) 



As we will discuss in the next chapter, next-to- leading order corrections in the 1/Nq expansion 
seem to push that value down to Bk ~ 0.4, at least in the strict chiral limit. Taking into 
account that lattice estimates suggest a much bigger value, say Bk ~ 0.8, one concludes 
that chiral corrections are clearly not negligible and may account for the bulk of Bk- Recall 
that this parameter determines the prediction inside the Standard Model for the kaon mass 
difference and ek, through 



M 



K, 



Re(K°\n s e f f 2 (0)\K°] 



M Ks 

e K ~ lm(K o \n!jf(0)\K ) (4.65) 
At the same time, at least in the chiral limit, there is a relation connecting AS = 2 and AS = 1 



transitions. Indeed, 



(AS= 1) x (AS=1) 



comes as a AS = 2 long distance effect, as already 



pointed out. Since the AS — 2 operator transforms as a (27 l, 1r) under the chiral group, it 
follows that it is closely related to g27, namely |35l I69| 



.927 — tBk 
5 



(4.66) 



Next chapter is devoted to the determination of the AS — 2 parameters mentioned above in 
the chiral limit including 1/Nc corrections. 



Chapter 5 

Determination of kaon mixing 
parameters 

5.1 The Problem of Matching in Nonleptonic Weak In- 
teractions 

In preceding chapters we have seen that there is a dual description of Standard Model physics. 
On the one hand, one has the quark-gluon description, accepted as fundamental but of very 
little use beyond the perturbative regime. For most practical purposes, we are interested in 
computing matrix elements of hadronic states, and one is faced with an obvious problem of 
language mismatch. Consider, for instance, 

(K \s Llf ,d L s L rd L \K ) (5.1) 

Computation of the previous matrix element is far from being straightforward, for we do not 
know how quarks and gluons assemble to conform K states. Possible way-outs are numerical 
simulations (lattice QCD), quark models (Constituent Quark Model, ...) or effective field 
theories. We will adopt the latter strategy, since we want to be analytic until the very end and 
approximations have the advantage over models that they are systematically improvable. 

In previous chapters we have already shown the systematics of EFT's and showed that 
indeed there exists an alternative formulation of the strong interactions in terms of mesons 
under the form of chiral lagrangians. It is not known how to do the transition (sometimes people 
speak of change of variables) between the quark-gluon description and the meson description. 
However, we saw that symmetry alone was able to relate the quark-gluon operators to their dual 
counterparts in terms of mesons. This change of language is often referred to as bosonization 
of operators, and it is precisely what we need to compute matrix elements such as (|5.1H . Our 
remaining ignorance as to how these two worlds are connected is parametrised in terms of a set 
of low-energy couplings, one for each operator. Therefore, after bosonization, computation of 
hadronic matrix elements is tantamount to determining the low-energy constants, something 
that unavoidably requires a matching between the quark-gluon and meson descriptions. 

In full generality, low-energy couplings are related to QCD correlators. We can distinguish 
between two different behaviours, depending on whether they are strong low energy couplings 
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or electroweak low energy couplings. The parameters of the strong chiral lagrangian modulate 
the dynamics of strong processes in the presence of external currents. As already stated, they 
contain the information of the integrated degrees of freedom up to ~ 1 GcV. Gcnerically, they 
are expressible as the coefficients of a Taylor expansion of QCD correlators at low energy. For 
instance, consider the paradigmatic two-point functions 



with the QCD currents given by 



J£ = d(x) 7" u(x) , 



J I' 

J A 



d(x) 7 M 7 5 u(x) 



We define 



tLR, 



(5.2) 



(5.3) 



(5.4) 



The last line follows from Lorentz invariance and H LR (q 2 ) has a well-known low energy expan- 
sion 

Q 2 IWQ 2 ) = -F 2 + 4L 10 Q 2 + O (Q 4 ) (5.5) 



which yields 



and 



F 2 = 



Q 2 IIlr(Q 2 ) 



d 



Q 2 TIlr(Q 2 ) 



(5.6) 



(5.7) 



Q- 



Consequently, derivatives of Green functions which are order parameters of spontaneous chiral 
symmetry breaking lead to low energy couplings of the strong chiral Lagrangian. 

However, unlike what happens in the strong sector, the low energy couplings arising in the 
electroweak sector are much more involved. As we saw in chapter 2, they come from the virtual 
exchange of W and Z particles. Before reaching the non-perturbative regime were bosonisation 
is required, they have long been integrated out. As a result, the low energy couplings turn out 
to be integrals of QCD correlators. Following the paradigmatic 11^^ used above, one can see 
that the coupling c~2 introduced in chapter 2 as the lowest-order electromagnetic low energy 
coupling is expressible as 



('2 



32tt 2 



dQ 2 1 - 



Q 2 



Q 2 



Ml 



Q 2 ^lr{Q 2 



(5.8) 



Contrary to the strong case, where the operators were truly products of conserved currents, in 
the electroweak sector, even though they look like products of currents, they are convoluted 
currents (in the example above, due to the Z and photon exchange). Therefore, whereas in 
the strong sector the determination of the low energy couplings demands knowledge of the low 
energy properties of certain Green functions, in the electroweak sector the full Euclidean range 
is needed. In view of the fact that an exact solution of the problem is out of reach, for it 
requires full-fledged QCD, we shall have to resort to some approximations. However, whatever 
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approximations we use to the computation of matrix elements, they have to make sure that 
the scale dependence of the couplings and the operators entering the expression 

e 2 F^ 2 (p){Q L Q R ){^) (5.9) 

must cancel in order to provide a sound output for any physical observable. Thus, the challenge 
is not only to find a reliable approximation to the integral 15. 8|) above, but one that also 
guarantees the right matching between short and long distances. In the following sections we 
will show how to get the right matching. The method basically consists in approximating 
the integrals of Green functions such as that in 1)5.8(1 by integrals over interpolators based 
on the MHA we introduced in chapter 3. Recall that such interpolators were determined 
through matching to the known QCD high and low energies of the QCD Green function we are 
approximating. It is precisely the right OPE behaviour of our approximant which, eventually, 
automatically ensures the right matching between long and short distances in nonleptonic weak 
interactions 1 . 



5.2 Determination of kaon mixing parameters 
5.2.1 Lowest order Lagrangian 

Our purpose is to find the Effective Hamiltonian which describes K° — K° mixing. Our starting 
point will be the box diagrams of figure ((4.1(1 . Following the steps outlined in previous chapters, 
we have to integrate out the heavy degrees of freedom sequentially, beginning with the top 
quark and going down in energy. For the sake of simplicity we here reproduce the Gilman- 
Wise calculation [S3] and consider the W gauge boson heavier than the top quark, since this 
will not bring any qualitative difference. At the end we will modify our result accordingly to 
account for the top quark being actually heavier. This means that in our box diagrams every 
W boson propagator will be shrunk to a point, i.e., 

(q 2 « m 2 w ) (5.10) 



q 2 - mly 



since the momentum scale is indeed much lower than the particle mass. Box diagrams, after 
W and top integration, which again can be accomplished by shrinking its propagator, turn into 
the diagrams depicted in figure ((5.1(1 . and we are left with only one four-quark operator 

1 Fot different applications of the method, we refer to I7UI - I77I . 
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C 



eff 



Gp 

16^ 



c{n) \s(x)^(l- 75 )d(x] 



(5.11) 



whose Wilson coefficient at the top mass threshold /i = m t is given by conventional perturbation 
theory techniques, to yield 

c(m t ) = -m 2 X 2 (5.12) 

This is precisely the matching condition that guarantees that the Effective Lagrangian (|5.11() is 
equivalent to the full-fledged theory below the top quark mass threshold. From the top quark 
threshold until the charm quark threshold, the Wilson coefficient has to run as dictated by its 
rcnormalization group equation, 



which can be easily integrated to yield 

2 

c{p) = -m\ X 2 + 2m 2 c X c (A c + A„) log -§■ 



(5.13) 



(5.14) 



Recall that in the previous expression we have set m u = 0, which is a reasonably good approx- 
imation. The next step is the integration of the charm quark. Disappearance of the charm 
from the effective theory requires that the Wilson coefficient picks an extra piece, 



< 



= - m t \ 



+ 2m 2 c X c (A c + A„)log— | 



to? 



2 \2 

to„ A„ 



(5.15) 



Again, the matching condition has been easily determined by using perturbation theory. We 
are still above A x and we can run down the Wilson coefficient. However, its evolution below 
the charm mass threshold is proportional to m u . Thus, to a very good approximation, the 
Wilson coefficient gets frozen at its /i ~ m c value. The full answer for the dimension-six quark 
operators is then given by the well-known Gilman-Wise Lagrangian 



r ^f / 



2ml x cXt log ^ - ml A 2 



TO? 



(5.16) 



In the Effective Field Theory approach the non-logarithmic pieces above are to be interpreted 
as the matching conditions at the top quark and charm quark mass thresholds, whereas the 
logarithmic contribution stems from the running between the two scales. 

As mentioned before, the step we made of integrating out almost simultaneously the W 
gauge boson and the top quark, i.e., the assumption of considering the W slightly heavier than 
the top quark can be mended without changing any of the conclusions we reached with the 
Gilman-Wise Lagrangian (|5.16|1 . The full-fledged result for the (leading) AS" = 2 Effective 
Hamiltonian, once the particle hierarchy is respected and electroweak and strong corrections 
are included, reads 



H 



S=2 
eff 



G 2 



c(/i) < ?7i m 2 A 2 + r) 2 (mt) eff A 2 + 2n 3 m 2 A c A t log — | 



+ h.c. (5.17) 
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The appearance of an effective top quark mass, defined as |84| 



'eff 



2.39 M 



mt 



167 GeV 



1.52 



(5.18) 



is the ammendment to account for the top quark being heavier than the W gauge boson. The 
remaining r/j and c(/j) functions appear after perturbative corrections are taken into account. 
Their values are [EI] 



771 = (1.32 ±0.32) 



1.3 GeV 



m c (m c ) 



1.1 



'/2 



c(/x) = K(/x)) 



1 



a a (n) /1433 



1936 8 



0.57 ±0.01 



1 

— K 



0.47 ±0.05 



(5.19) 



where k is scheme dependent, its value being or —4 depending on whether one is working in 
the naive dimensional or 't Hooft-Veltman regularization schemes. 

5.2.2 B K with leading short distance OPE constraints 

The evaluation of the K° — K° matrix element, as already pointed out, is tantamount to a 
determination of the bag parameter 



(K°\c<jM)Q(0)\K°) = -F^m 2 K B K 



(5.20) 



The computation of the invariant Bk is however not as straightforward, since there is a mis- 
match in the above equation, as stressed at the beginning of the chapter. One should be able 
to express the four-quark operator in terms of meson operators, or, in other words, to know 
how the four-quark operator changes when strong interactions bind quarks and gluons together 
so close that they appear as hadrons. 

As already emphasised, symmetry alone is capable of constraining this set of meson op- 
erators. As with every Effective Field Theory, there are a number of low-energy coupling 
constants that have to be determined through a matching procedure. Fortunately, there is 
only one AS* = 2 operator in the leading chiral Lagrangian, to wit 



r 2 

X ~ 16^2 ^0 A S=2 11 



X 32 (D^)UX 32 (D^)U 



(5.21) 



where A| =2 is the a priori unknown low-energy coupling constant. For convenience let us 



factor out the short distance content of A| =2 , 



A 2 



9S=2 



r/i m\ \\ + r/ 2 ( m t) eff + 2 % m 2 c X c \ t log 



i 1 
,2 



(5.22) 



Plugging in the previous equation in the definition of the invariant Bk, one gets 



Bk — t 9s=2 



(5.23) 



So far, we have changed the problem such that our determination of the bag parameter is now 
tantamount to knowledge of the low-energy coupling constant gs=2 (or equivalently ^$=2) 
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governing the AS = 2 transitions. The change in vocabulary has been so abrupt in the 
bosonization process that the matching procedure is, unlike the ones we have encountered 
when deriving the Gilman-Wisc Hamiltonian, highly non-trivial. 

In chapter 2, when dealing with chiral Lagrangians, we emphasised the convenience of in- 
troducing external fields in the formalism in order to be able to compute QCD Green functions. 
Here we will make use of its useful consequences. A close inspection of l|5.21|) reveals that it 
contains a quadratic term in the right-handed external field. Indeed, recalling that 



(5.24) 



it is not difficult to get 



A3 2 ( J D^ t )(7A 32 (^C/ t )C/ = • • • - r^rf + 



(5.25) 



The right-handed current is by construction the same one showing up in the QCD La- 
grangian. Furthermore, since this quadratic term in the right-handed external fields is a gen- 
uinely AS" = 2 operator, it can solely come from the chiral Lagrangian Q5.21 Jl . The matching 
condition is therefore imposed by comparing a QCD Green function with the Green function 
coming from the effective chiral Lagrangian, namely 




C 



QCD 



r S=2 
L eff 




(5.26) 



Equating the right-hand side of l|5.26|l . coming from the electroweak chiral Lagrangian (|5.2I|) . 
with the left-hand side, computed with the effective Hamiltonian (| B . 1 T|> . one obtains jjl] 



ds=2 = c(m) 



32tt 2 F 2 



F(2 - |) 



dzz ^W(z) 



(5.27) 



where /i^ ad is a scale introduced to rewrite the integral in terms of a dimensionless variable 
z = Q 2 1 nj lad - W(z) is a short-hand notation for 



W(z) 



2 

9 ^ had 



1 



LRLR 



(5.28) 



and w2jiLji{zn\ ad ) is the solid angle integral of a four-point Green function with two left- 
handed currents and two right-handed (soft) currents 



9 



a/3 



w 



(1) 



LRLR 



(Q 2 ) 



Q 



(5.29) 



unfactorized 

where the four-point Green function in the left-hand side is defined as 



IT', 



LRLR 



(9,0 



lim i 3 



(fx <Ty (fz e lqx e 



iqx il-(y-z) 



(0\T{Ll d (x) Rt s (y) Ll d (0) Rf(z)}\0) (5.30) 



the currents being 



75 



d(x), R%{x)=d{x)^ 



•75 



s(x) 



(5.31) 
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Figure 5.2: Resonance diagrams contributing to W(z). 



This four-point Green function with two additional hard left-handed insertions is one of the 
QCD Green functions underlying AS = 2 processes, but certainly not the only one. However, 
it is relatively simple and has the nice feature that it is an order parameter of spontaneous 
chiral symmetry breaking (S^SB), which is more than welcome 2 . Furthermore, its physical 
content is rather transparent: from its very definition, our four-point function is a two-point 
correlator of left-handed currents, which guarantee the K° — K° transition, together with two 
soft insertions of right-handed currents required by the matching condition (|5.26|) . 

Equation (|5.27|) is an example of what we anticipated at the beginning of the chapter, i.e., 
low energy couplings of the electrowcak sector can be expressed as integrals over the Euclidean 
of certain Green functions. As already outlined, the procedure hereafter will be the following: 
we collect as much information about the Green function as possible, on the long distance 
regime by means of chiral perturbation theory and on the short distance regime by computing 
the OPE coefficients. We then bridge the energy gap in between by means of the large-Ac 
version of the Green function. This constrains the analytical form of the Green function to be 
a meromorphic function. The MHA is then used to truncate the infinite number of mesons to a 
finite number of resonances where the undetermined parameters (poles and residues) can then 
be fixed by matching onto the low and high energy behaviour of the original Green function. 

The analytical large- Nc behaviour of the four-point function W£^£ R is depicted in figure 
l|5.2[l . where all possible resonance exchanges are taken into account. This constrains our 
Green function to be a sum over an infinite number of resonances of at most triple poles if the 
resonance is a vector or an axial vector and simple poles if they are scalars or pseudoscalars. 
Therefore, in the MHA to N c -> oo QCD, 

W(z)=Y (-^- + -^— + -^—) + Y-^- (5.32) 

where pi = mf / ^f lad is a normalized mass, to; being the mass of the resonance i; o,-, bi, c, 
are the constants to be determined by matching onto the high and low energy regimes of the 
Nc = 3 QCD Green function. 



Long Distance Constraints 



2 Recall that order parameters are preferred if we eventually want to use the MHA (see discussion in section 
3.3). 
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Contributions to be taken into account come from the diagrams coming from the Chiral 
Lagrangian connecting two left-handed and two right-handed external sources. We quote the 
result from |71) 

W xPT {z) =6-24 ^ ^2Lx + 5L 2 + L 3 + L^j z + O (z 2 ) (5.33) 
The more involved part will be the calculation of the OPE coefficients. 



Short Distance Constraints 



Since our four-point Green function has the right-handed currents as soft momentum inser- 
tions, one needs to consider just the operator product expansion of the left-handed currents. 
Thus, this separation of scales allows to factorize the problem to the more manageable OPE 
two-point computation. 

The diagrams one has to consider are the ones shown in figure (|5.3|1 . where q 2 is a high 
momentum insertion which flows through the diagram as a gluon exchange in all possible 
ways. This way we are computing the leading corrections to K° — K° mixing in a s . Obviously 
there are also electroweak and electromagnetic corrections which enter the Operator Product 
Expansion, but they are clearly suppressed when confronted with the strong interactions. 

A convenient tool to deal with OPE calculations is the use of the so-called Schwinger 
Operator Formalism |851 1861 |87 j . Essentially, it is a background field method in coordinate 
space with the ingredient that covariance is preserved at each step of the computation. A 
rather detailed introduction into the subject is given in an appendix. 

We set off from the AS = 2 four quark operator. The four diagrams of figure l|5.3l) come 
out as an expansion of the product of currents 

s L {x)^d L {x) s L (0)7 M d L (0) (5.34) 

Following the steps outlined in Appendix B, the previous equation can be cast in the form 



.(f)) 



,(0), 
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d 4 x (x) ig s G^{x) 7 M S(x — x) 



d A x (x) ig s G^{x) 7 M S(x — x) 



l,df{x) 



d?V$\v)ig.G ll (,v)>fS(0-y) 



7m4 0) (0) + 



+ 0{a 2 s ) 



d 4 yS(0-y)ig s G,(y)rd ( °\y) 



(5.35) 



where the first piece is the unconnected one. We are interested in terms proportional to a s , 
i.e., the one gluon exchange. Using the expansion for the gluon and quark propagators given in 
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appendix B, we can systematically extract the different pieces of the operator product expansion 
order by order in powers of the soft momentum P and inverse powers of the hard momentum 
q. In particular, the P° terms will give rise to the dimension six pieces of the OPE, while the 
P terms will yield the dimension eight pieces, and so on. Thus, formally, we can determine 
the OPE to whatever order. However, we eventually would like to compute matrix elements 
of the resulting operators and, on the phenomenological side, regretfully, little is known about 
the matrix elements of such operators of increasing complexity. Nonetheless, one expects their 
impact to be progressively smaller, due to the suppressing 1/q 2 powers. 

In the MHA framework, consistency demands a minimum of terms coming from short 
distances. In the present case, the leading dimension six contribution to the OPE would suffice 
to constraint an interpolator with one resonance state, either vector or axial vector. This 
is precisely what was done in |74| . However, shortly afterwards there were claims that the 
neglected dimension eight operators could yield a potentially sizeable effect [IB] , thus affecting 
the previous determination of Bk- This motivated a reassessment of Bk with the inclusion of 
next to leading terms in the operator product expansion. 

The full result for dimension six contributions is the following, 

[ d*xe^ T{Ll d (x)L» d (0)}=J2 c f ] iQ 2 ) Of (5.36) 

•* i 

where there is only one operator, namely 

(6) = s L (Q)^d L (0) s L (0)7^(0) (5.37) 
with the Wilson coefficient given, to lowest order in a s = gl/An, by 3 

c^(Q 2 )=i^ (5.38) 

The presence of just one operator is precisely what one expects: clearly it is the only possibility 
for AS = 2 dimension six operators. We will briefly see that this simplicity no longer holds 
and combinatorics arc larger in dimension eight operators. Plugging in the previously found 
OPE for the left-handed currents into the full four-point Green function, one gets 

W^ LR (q,l) = Umi 3 J d 4 yd 4 ze il ^\0\T^^\Q 2 )Ol 6) Ri s (y)Rf(z)\ |0) (5.39) 

which can be evaluated first using vacuum saturation and then inserting a sum over complete 
states in all ways allowed by quantum numbers. In the large- Nc limit single particle states 
are dominant over multiparticle states, and, since by definition the momentum I is soft, this 
singles out the kaon. The result one gets is, 

W^ RLR (q,l) = 2c e \imi 3 (0\s L ^d L \K(l))±(K(l)\d RlaSR \0) 

(0 \s L7fl d L \K(l)) i (K(l)\d RlpSR \0) (5.40) 

3 See Appendix B. 
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which, together with the definition 

(0|«j/vfc|tf(0) = -iV^^oip 

immediately leads to [71j 



IT" 



OPE 



24 TOs f 2 l r 



1 



12 



(5- 



(5.41) 



(5.42) 



We are provided with two conditions coming from the low-energy side and one condition com- 
ing from the high-energy side of our Green function. Since we have already determined the 
analytical form of our four-point Green function, we can apply the Hadronic Approximation 
consistently if we include at least one particle, i.e., 



a \ ■ 



Cy 



(z + pv) (z + p v ) 2 (z + pv) 3 



(5.43) 



The authors of choose the particle to be close to the real p(770) state. After fixing the 
residues ay , by and cy they found a value for 



B k = 0.33 ±0.11 



(5.44) 



which is in agreement with other estimates [751 EDI an d with chiral extrapolations of lattice 
simulations 1811 1821. 



5.2.3 Inclusion of next-to-leading short distance constraints on W^f R 

Inclusion of next to leading terms in the OPE modifies the expression <|5.36ll in the following 
fashion 

f d^e^T{Lf(x)L^(0)}=Y^ cf (Q 2 )Of-^ cf >(Q 2 ) Of + ... (5.45) 

where the new set of dimension eight operators can be straightforwardly determined (see Ap- 
pendix B) and read jHHl 



0\ 



(8) 



sV^Vld sT v a d + sT a v V^dsT v a d 



sT^dsV^V^ T a y d- 



sT^dsTlV^d 



sT u a V^d sTlV^d + sV^dsV^T^d , 
s% T v a d sT a v V» d + sT a v V^d sV^T v a d, 
s% Y v a dWD v V^d + sT a v V^d sT» V v d , 
s%,T v a d sT^V u d + sT^V^d sTTT^d , 
9s G% {sT^dsVd-sT^dsTld} 



Tff is defined as 



2 2 
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(5.46) 
(5.47) 
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and the following conventions were adopted 



Da = da- igsAfj. , = 



Gfiu = dpAv-dvAp-igAAvA,,] , G uu = 1 e prTUU G pa , with e 0123 = +1 (5.48) 
The corresponding Wilson coefficients are 

C| 0) =i f^t ^00 (5.49) 



Q 6 

where the rji are c-numbers, whose values are 

(8) _ 5 ( 8 ) 22 (s) _ 8 

Vi — g j % — g i % ~~ g 

„(«) - 15 „(«) - 15 (8) _ J_ . . 

The previous basis of operators satisfy a CPS symmetry inherited from the box diagrams we 
started from. Indeed, consider the diagrams in position space. Their expressions are 



Ti = / d i xd 4 x'd 4 zd 4 z' 



T-2 = I d^x d 4 x' d A z d 4 z' 



SLW-fJPV) q { l\x')^d L {x') s L {x) lv q%\x) $ \z) lv d L {z) 

■D(x - x') D(z - z') 

s L {z')^qf{z') <tf\x'h»d L (x') s^xh^ix) q£\z)<y v d L (z) ■ 

■D(x-x') D(z- z') (5.51) 



where D(x—x') is the W propagator between x and x' . A CPS symmetry is a CP transformation 
followed by an (s <-*■ d ) interchange. Using the well-known charge conjugation transformation 
algebra 

C- 1 7" C = -{^) T , C- 1 7" 75 C = ( 7 ^ 75 ) T (5.52) 
and the one for the parity transformation 

P" 1 P = (y)t , P 1 7 ^ 75 P = ( 7 " 7s )t (5.53) 

one can easily arrive at the CPS transform of the Dirac bilinear 

s L (x) Y l q L (x) -q L (x) 7 M d L (x) (5.54) 

where x is the parity-transform of x. Integration measures should also be transformed accord- 
ingly. The conclusion one reaches is 

Ti <-> T 2 (5.55) 

Therefore, effective operators at each order have to preserve this original symmetry. The 
dimension-six operator trivially does. A bit more involved calculation shows that the dimension- 
eight operator basis lj5.46J) is indeed CPS invariant. This is a non-trivial check which, in 
particular, the basis found in does not respect. 
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The steps to follow hereafter are the same as with the dimension six contribution. How- 
ever, as we already pointed out, the situation now is much more involved, since we have more 
operators and their matrix elements are a priori unknown. Use of the large-iV^ limit here 

(8) 

is extremely advantageous. First, one notices that the contribution coming from Oq is sup- 
pressed due to the Wilson coefficient having a 1/Nc out front. Second, after colour-ficrzing 
the operators appearing in 0\ to 0\ , one ends up with matrix elements of the form 

(0\S L ^D v d L \K(l)) = F^l 2 ^ + F 2 {l 2 )l^l v (5.56) 

whose tensorial behaviour can be inferred from Lorentz invariance. Contracting with the metric 
tensor g^ v and using the equations of motion in the chiral limit, i.e., 

Bd L ~O(m d )~0 (5.57) 

one concludes that 

Fx(/ 2 ) ~ I 2 + 0(l 4 ) (5.58) 

which means that contributions of of\ 0< 8) , Of\ of ] are suppressed by I 2 with respect to 
the leading terms. The remaining 0\ operator can be evaluated by inserting the identity 

D 2 = B 2 + l -g s a^ v , <r M „ = (5.59) 
which after some algebra (see Appendix B) leads to 

(0\s L ^D 2 d L \K(l)) = (0\g s s L G a ^X a ^d L \K(l)) = (-*V5f l v ) 5% (5.60) 

It is worth emphasising that the previous equation defines the dimcnsionless 5 2 K parameter. 
The notation is due to where the square was chosen a posteriori to stress that 5 2 K > 0. By 
comparing with 1)5. 41J) . we see that 5^ measures the relative strength of the dimension-eight 
quark gluon operator with the dimension-six operator. It turns out, as we will show in the 
next section, that 8\ ~ 0.12 GeV 2 . 

The last equation, in particular, means that 0\ is dominant in the chiral expansion, i.e., it 
scales in the soft momentum / as the dimension six term does. Therefore, of the six operators 
that conform the basis of the dimension eight operators, only the first one survives in the 
large- Nc and chiral limits. We have thus reduced the complexity of the problem by reducing 
the number of operators to a single one, whose matrix element is governed by the parameter 



5.2.4 Determination of b\ 

We will derive a sum rule which will allow us to extract the value of the parameter b\. It 
will be a good strategy to consider the whole combination 5 2 k Fq as a single parameter. We 
have already mentioned that Fq is an order parameter which appears in the well-known two- 
point correlator 11^. It is a good idea to follow this analogy and define the two-point Green 
function 4 

n^(Q 2 ) = if d^x e^ x (0 \T {^-s L G a ^ a d L (x) d Rl „ SR (0)} |0) (5.61) 

4 This came as a suggestion of Eduardo do Rafael in a private communication. 



5.2 Determination of kaon mixing parameters 



72 





Figure 5.4: Diagrams determining the lowest order contribution to the OPE of fS.fil)) . 



in which, much the same as with 11^, Lorentz invariance guarantees that, in the chiral limit, 
it has the following tensorial structure 



The large-g 2 fall-off of the function l|5.62[) can be straightforwardly computed to be 



n Li V) 



2ir a s < i/iip > 2 



Q 



9 Q 4 

where factorization of the four-quark condensate in the large- N c limit, .i.e., [HO] 

1 



(o|#^#^|o) 



N 2 



Tr Tr (r,-) - Tr (T t T 3 



(5.62) 



(5.63) 



(5.64) 



has been used. The normalization factor N is defined as 

(Ofe|0) = ^ WVIO) (5.65) 

such that the above expression is an average over colour and spin degrees of freedom. Diagrams 
leading to the previous result are depicted in figure (5.4). Using an order parameter as a Green 
function is a good choice in the sense that it considerably simplifies the calculations, since 
there is no contribution from the continuum and perturbation theory cancels to all orders. 
The usual strategy in Finite Energy Sum Rules (FESR) is to make an hadronic ansatz for 
the spectral function and relate it through a dispersion relation to the values of the Operator 
Product Expansion 5 . Our hadronic ansatz will be 



hmU LR (t) = -^^S(t) 

7T 4 



JV°V xn 2 \ 
-^S(t-m v ) 



and the dispersion relation that Ulr obeys is an unsubstracted one, 

1 1 



dt 



5 See 1911 for a detailed review. 



t + Q 2 



Imn LH (t) 



(5.66) 



(5.67) 
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much the same as it happens with Ulr- Plugging (|5.66(1 in (|5.67(1 and comparing with 
we end up with the following Weinberg-like sum rules 

F 2 8 2 K m 



4 

4^ = f « s <^> 2 (5.68) 



from which the unknown S v and 8 2 K parameters are readily determined to be 

167T a s < 'iptp > 2 



S 



2y 9 f v m\ 



*« = T FfmJ (5 - 69) 



Using as input parameters 

F ~ 0.087 GcV , / v ~0.15, 0.77 GeV, 



a s {2 GeV) ~ 1/3, (^)(2 GeV) ~ -(280 ± 30 MeV) 3 (5.70) 

one obtains 6 

5 2 K = (0.12 ± 0.07) GeV 2 , S v = (0.06 ± 0.04) GeV 4 (5.71) 

As a cross-check of stability of Ij5.71|l . we could construct a Laplace sum rule 7 . Laplace Sum 
Rules have the advantage over Finite Energy Sum Rules that they exponentially suppress the 
contributions from high energy resonance states, such that in principle our truncated hadronic 
ansatz would be more accurate. Taking 

hmU LR (t) = S(t) + &*V6(t - m v ) + %&S(t - m\) (5.72) 

as our spectral function and applying the Laplace transform to (|5.72J) and to the OPE of H5.63J) 
we find the following Laplace sum rule 

p2r2 f2r2 f 2 r2 o_ 

_£0^ + JA°A e -mlT + W>V e - m » v T = _^L as < ^ >3 T (5.73) 

4 8 8 9 

Notice that in this Laplace sum rule analysis we have included an axial resonance in our spectral 
ansatz (|5.72(l . Contrary to what happens in Finite Energy Sum Rules, where one has a finite 
number of relations, Laplace sum rules appear as functions of the Borel parameter r. Since the 
parameters one wants to extract should not depend on the Borel parameter r, one searches for 
stability in the Borel parameter space. Taking as input l|5.70|l . together with 

f A ~ 0.08 , m A - 1.2 GcV (5.74) 

6 We have added generous error bars in the quark condensate to include the present spread of values in this 
quantity. This error is the dominant one. 

7 We refer to the review on spectral sum rules of 1911 . 
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Wc find the Laplace sum rule to hold in the rather wide window < r < 1 GeV 2 , as shown in 
figure (5.5), for values of 8\ and 5 V compatible with the ones in l|5.71|l . together with 8\ ~ 0.05 
GeV 4 . 

The result we have obtained for 8\ is consistent with two previous analysis jHOl IH2I ■ They 
both determined this 8\ parameter governing dimension-eight effects in kaon mixing with a 
slightly different Green function, which is not an order parameter. This makes the result 
depend on the onset of the continuum Sq, and some discrepancy arose between both analysis as 
to the right value for sq. In our opinion, the use of 11^^ as the Green function to perform the 
analysis has the advantage of not having to introduce this extra parameter and thus making 
the extraction of parameters much cleaner. 

5.2.5 Numerical results for Bk 

Armed with a precise knowledge of the matrix clement relevant for dimension eight contribu- 
tions, we can parallel the analysis we performed for the dimension six operators and determine 
the next to leading OPE contributions 

W^qM&B = * 3 / *V rf t *e*(»-)<0|r{ £cf> (Q 2 )oWRi°(y)Rf(z)}\0) 

(5.75) 

Repeating the steps we followed on page 68, we get 

W^ RLR (q,l) = Jim e2cf\Q\g s s L G^ lp d^K{l))^{K{l)\d RlaSR \Q) 

(0\s L7fi d L \K(l)) i (K(l)\d R y s R \O) (5.76) 
Putting all things together, we obtain 

_ i r e + 10^ n + o / n 

vLd z L 12 9 Kad z \ \ z3 J 
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for the OPE up to dimension eight contributions. As expected, dimension eight effects are 
modulated by the b\ parameter (compare this expression with l|5.42)l ). Inclusion of these 
higher dimensionality effects results in an additional condition onto our interpolator (|5.32(1 . 
The first possibility is to try to fix the pole of the vectorial resonance we had before. This leads 
to an imaginary value for the mass, which clearly calls for the introduction of a new resonance 
in the interpolator to accomodate the new constraint. The simplest possibility is to introduce 
a scalar particle close to the /o(980) — ao(980) doublet. 8 
Therefore, the interpolator we choose has the form 



W(z) {v . s} 



b v 



c v 



(5.78) 



(z + pv) (z + Pv) 2 {z + pvY (z + p s ) 

where the free parameters ay, by, cy, as can be determined by imposing matching onto the 
two constraints coming from long distances and the two more constraints coming from short 
distances. This results in the following 4 constraints: 



ay + as = 



ay 
Pv 



by - a v py - a s p s 

by Cy 



h^- + ^- + — = 6 



Pv 



as_ 



a V by Cy t a S 

o T i o T O , 
P V P V Pv 



Pi 



24to s F 2 

Phad 

24ira s F 2 

Phad 



1 



10 



K 



24^( : 



(5.79) 



where the second one is the one coming from next to leading order OPE operators. Before 
giving the numbers for Bk , let us note that the correction we are computing can be understood 
qualitatively by the following estimation 



5gs=2 



Phad 

327r 2 F n 2 



24-Ka s (p had )F§ 



Phad 



10 



Phad 



dz 

~2 



15 a s (phad) 

18 7T 



P 



ope 



(5.80) 



where p, ope is the scale above which the large-z expansion (i.e. the OPE) starts making sense. 
The previous expression emphasizes the fact that corrections due to the inclusion of dimension- 
eight operators in the OPE are governed by the ratio between b\ and a typical scale hope- 
Numerically, one obtains |<5<7s=2| S 0-03 when pL ope ~ 1 GeV and the value for 8\ in Eq. 
(I5.71|) are used. Plugging in our interpolator in 15.27|) and performing the integral one gets 
the determination for gs=2, which can be shown to be scale and scheme independent. 5s=2 is 

8 A similar analysis introducing an axial vector particle was also done as a consistency check. Since an axial 
particle introduces three new parameters to determine, we decided to constrain the interpolator by requiring it 
to be equal to the OPE at two arbitrary points sufficiently high in energy for the OPE to be trusted. Results 
as compared with the ones obtained using a scalar particle differ by less than a 5%. 
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Figure 5.6: Plot of the large- and small-z expansions of the function W(z) as given by the OPE and 
Chiral Theory, respectively (dashed curves). The solid curve corresponds to the interpolating function 
W(z)(v,s) obtained in Eg. \5. 751 ). 



trivially related to Bk , and the final result is the following 



B K 



Us{Hhad) 



f'had 



a s {u. had ) 1229 



1936 



O 



N c a 2 s (u, had ) 



32tt 2 F 2 



1 c v 

-ay log — as log ps H h - — 

py 2p^ 

which is indeed scale and scheme independent. 

This is our analytical result for Bk in terms of the residues ay, as, by and cy, to be 
determined using the matching equations of H5.79(l . Putting numbers. 



(5.81) 



a s (2GeV) = 0.33, F = 0.087 MeV, m v = 0.77 ± 0.03 GeV 



m s = 0.9 ± 0.4 GeV , 1 GeV < p had < 2 GeV (5.82) 

we get 

B K = 0.36 ±0.15 (5.83) 

for our determination of the invariant Bk- Comparison with 15.44fl shows that including next 
to leading order terms in the OPE does not affect much the determination of Bk- This is 
what one should expect, taking into account that Bk comes as the integral under the curve 
shown in figure (5.6). Changes in the high energy tip, for which the OPE is responsible, should 
therefore be dramatic to make the area below the curve change notably. Furthermore, we have 
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explicitly shown that the parameter controlling the next to leading OPE contributions is 
rather small. 

As already discussed at the end of the previous chapter, in the chiral limit Bk is related 
to one of the electroweak AS — 1 couplings, namely <?27- Therefore, as a by-product, we get a 
prediction for (727, 

327 = -=B K = 0.29 ±0.12 (5.84) 
5 

which is in excellent agreement with a recent determination from K e % decays |93) . 



5.2.6 Corrections to the AS = 2 Lagrangian in inverse powers of the 
charm mass 

As pointed out at the beginning of the chapter, there are two main contributions to kaon mixing: 
one coming from direct AS = 2 transitions, which is conventionally named under short distance 

contribution and a long distance contribution arising from non-local (AS* = 1) x (AS = 1] 

terms. We already commented on the fact that the latter is estimated to be small. Short 
distances are meant to be responsible for the bulk of the kaon mass difference, especially 
through the charm mass contribution, and indeed it is so. However, there are long and short 
distance contributions arising from the next to leading order operators in the 1 /mj? expansion, 
where convergence can be dubious and so the impact of these next to leading operators sizeable. 
Our objective is to assess their actual impact. 

We already stressed in chapter 2 that one of the requirements for the existence of an effective 
field theory is the existence of a separation of scales large enough to be able to separate heavy 
from light degrees of freedom. We already mentioned the potential hazard of considering 
m s < A x as to the convergence of the chiral expansion. When we derived the AS* = 2 Effective 
Lagrangian at leading order, we were implicitly assuming that m c 3> A x , which is not true in 
the real world. Thus, it seems quite unjustified to truncate the Effective Lagrangian at leading 
order, since next to leading terms in the form of corrections in the inverse charm mass might 
as well be sizeable. By dimensional analysis, l/m^ corrections in the effective Lagrangian will 
come as dimension eight operators. Our purpose will be to check if the expansion in the 

AS* = 2 Hamiltonian converges good enough. 

From the discussion of the previous section, we already know the basis of dimension eight 
AS = 2 operators, and, since we will be working in the large-Ac limit, we can anticipate that 
this contribution will be eventually proportional to the mixed quark gluon operator 

g s s L G^d L (5.85) 

However, this is only the short distance contribution, which we can anticipate that will be 
proportional to the parameter 8\. There is also an additional long distance contribution, which 
comes out from computing the up-up contribution to the box diagrams. This contribution is 
expected to be much bigger and proportional to a hadronic scale A x ~ 1 GeV. Their interplay 
in the determination of the kaon mass difference Amjf and Ek is what we will try to determine. 

As in the derivation of the dimension six effective hamiltonian, we start from the box 
diagrams after W integration. Dimensional analysis guarantees that no mass terms can appear 
in the Wilson coefficient. Above the top mass threshold quarks are effectively degenerate and 
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the GIM mechanism is at work, yielding 

c 8 (m t ) = (5.86) 

In the energy region m c < fx < m t the Wilson cocfBcicnt captures the ultraviolet behaviour 
of top-top, top-charm and top-up diagrams. Due to the absence of masses, this contribution 
comes in the form 

A 2 + 2 A 4 A c + 2A t A u = -A^ (5.87) 

where the unitarity relation X t + X c + X u = has been used. Therefore, there is a GIM- 
induccd suppression in the Wilson coefficients down to the scale of the charm quark. The first 
non-negligible contribution is then 

c 8 (m c ) = - I X 2 U - H A C A„ (5.88) 
6 6 

at the charm mass threshold. In the equation above we have used X 2 + 2A C A U ~ — A„, which 
is consistent with taking At as negligible. Below the charm mass, the Wilson coefficient gets 
contributions solely from the up-up diagram, with the result |94| 



c 8 ( M ) = - — X C X U + Xi [log ^ - - j (5.89) 

Obviously, the above equation ceases to be meaningful below A x . Perturbation theory applies 
no longer and our effective Lagrangian reads 

CP 1 ~ Gp 

C-eff = c s(/«) 9sSLG^^„d L s L 7 M d L - s L l^u L ulJ^cIl + h.c. (5.90) 

where the second term accounts for the left-over up-up diagram. As it stands, this diagram is a 
purely non-perturbative object, whose ultraviolet behaviour has to cancel the scale dependence 
of ca(fJ^), since after all box diagrams are finite. The process of our calculation will eventually 
show this scale cancellation explicitly Observe that the evaluation of this long-distance con- 
tribution is the only difficulty it remains. The short distance piece can be easily evaluated in 
the largc-A^c limit, where factorization applies, once we know the relevant matrix elements. 
Not suprisingly, it turns out that the matrix elements we need are precisely the ones involved 
in the determination of the next to leading terms in the OPE for Bk, i.e., 

(0\s Ll)1 d L \K(l)) = -iV2F l^ 
(0\g s s L G a ^X a ^d L \K(l)) = (-iV2F l„)gjc (5.91) 

Let us therefore focus on the non-perturbative part of our analysis. As it stands in figure (5.7), 
the up-up diagram is a misrepresentation of the true picture, in which the up-up exchange is a 
cloud of confined gluons and quarks showing themselves as hadrons. As we did previously, we 
have to construct an appropriate Green function, schematically shown in figure H5.7(l , with two 
left-handed insertions and two soft right-handed ones 9 . In the large- Nq limit, the four-point 
Green function depicted in figure l|5.7|l can be split into a convolution of simpler three-point 
Green functions. This factorization of the Green function is illustrated in figure (|5.8J) . The 

9 This should come as no surprise. We are eventually computing parameters of AS = 2 transitions beyond 
the definition of Bk ■ Recall that Bk is defined as the kaon matrix element of the dimension six operator in 
the AS = 2 Lagrangian. 
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Figure 5. 7: Pictorical representation of the up-up exchange contribution. 




Figure 5. 8: Factorization of the up-up contribution leading to 



Green function we have to tackle is therefore 



(q,p) ee / d 4 xd 4 y e-^ x+ vy\Q[T{Ri s {y)L s ;{x)Ll d {Q)} |0) 



(5.92) 



where Ra(y), L s "(x) and stand for the QCD left and right-handed currents 



L™(x) = S(x) 7m ( ^) u(x) , Lf(0) = u(0) 7 , ( ^ ) 
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(5.93) 



The subscript A explicitly points at the fact that the Green function is affected by the axial 
anomaly, a feature which deserves further discussion. 

Anomalies in three-point Green functions were discovered more than thirty years ago in 
the pioneering papers of Adler, Bell and Jackiw 10 . The observation was that certain three- 
point functions did not obey the classical Ward identities and, even more worrisome, they were 
ambiguous in the sense that Feynman diagrams were rcgularization-depcndcnt. A much closer 
analysis revealed that this two facts are actually manifestations of the same phenomenon. The 
triangle diagrams are linearly divergent and depend on the regularization chosen. The simplest 
way to make this explicit is to compute the triangle with a momentum prescription and then 

10 See, e.g., 1951 for an exhaustive discussion of anomalies in quantum field theory. 
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shift it (see, e.g., jHSI for details). Therefore, defining T^ va as the result of the triangle in 
a given regularization prescription, any T^ va is also a solution of the triangle if they can be 
related as 

^iivgl — ^iivol 4~ ^p,va (5.94) 

the last term being a polynomial in momenta. Thus, this is the reason why the classical Ward 
identities are not fulfilled. It is conventional to choose the axial Ward identity to be anomalous 
rather than the vector Ward identity, even though we have just argued that one could as well 
do the other way round. This means that every triangle with an odd number of axial currents 
will be anomalous. Our three-point Green function indeed contains such anomalous pieces. 
Decomposing the left and right-handed currents 1)5.2.6(1 in terms of vector and axial ones 

L^x) = q{x) 7 M f 1 2 75 ^ q(x) = V M (x) - A^x) 

R^x) = q(x) 1 J^-^-)q(x) = V^x)+A^x) (5.95) 



where 



Vn(x) = -q{x) 7 M q(x) , A^{x) = -q(x) 7^75 q(x) (5.96) 



it is straightforward to check that 

+ V/iVuAa - V»A v V a - AftV„V a - A^A u A a (5.97) 

where the first line collects the non-anomalous terms, whereas the second contains the anomaly. 
It can be readily seen that ((5.92(1 satisfies de following (anomalous) Ward identities 

q^ a {q,p) = -in^(p)-i^ £ Wa » 
where II^j is our well-known two-point function 



nj&fcO = J d 4 xe-^y(0\T{R2ML» d (0)} |0) = [g ua - ^f) H LR (P 2 ) (5.99) 

However, we can now turn the previously discussed ambiguity to our advantage and choose the 
regularization prescription which makes the naive Ward identities valid. This can be done with 
a modification of the chronological product T (see, e.g., |96('l. Thus, making use of the inherent 
ambiguity of the anomaly we can always remove it from our analysis with a proper redefinition 
of the time-order operator T. To this end, we define a non-anomalous Green function 

T^ a {q,p) = f d 4 xd 4 y e-^ + Py\0\f{R d a s (y)L™(x)L: d (0)} |0) (5.100) 
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where T is defined such that r M „ Q satisfy, de fiat, the naive Ward identities 

q^™{q,p) = -»H^(p) 
(q+p)^ a (q,p) = *II{£(p) 

p Q r^ Q ( g , P ) = o 



(5.101) 



It can be easily shown that our non-anomalous Green functions is related to the anomalous 
one through 

N n . N, 



T^ a (q,p) = T% a {q,p) + l ^ e^q x + l ^ e^p x (5.102) 

This is the expression akin to l|5.94|) of our previous discussion. Having defined our three-point 
Green function to be anomaly-free, we turn to its characterization. The most general tensorial 
structure of T ^ va consists of 22 terms. This number is considerably reduced in the chiral limit, 
i.e., if we take p — > 0. The constraints of the above Ward identities (|5. 101(1 relate them to yield 
only two independent form factors Ii(Q 2 ) and ^(Q 2 )- The final form for our Green function 
is 

T^ a (<l,P) = J™ [*xd*ye-«^H0\T{l£(y)L?(x)Ly l (0)}\0) 



= U LR (0)T^ a + h(Q 2 )T*Z + h{Q 2 )T* a + 0( P ) 
where Q 2 — —q 2 , and Ulr(0) is the pion pole 



(5.103) 



IW0) = -i^f 



(5.104) 



Tavai TnJa an d ^^va are Quite lenghty tensorial structures 



p 2 q 2 



v^q^aav + p 2 q^9afi - PaP^q^ - p a Pnqv + qp.q v 



p ■ q 



n-ST (2 \ ( P ' q \ 



IS^vXaq 



r PaP 



*\q 



(5.105) 



where the superscripts make it explicit the symmetric properties of the tensors with respect 
to the hard momentum q. Thus, T^ va is symmetric under a, <-> v exchange, T^J^ is symmetric 
and transverse while T^ ua is purely antisymmetric. The expression for r /JlyQ ,(g,p) in its full 
glory then reads 

F^ a (q,p) 



Fo 2 



2p 2 q 2 

i2n (2 



p 2 qp.9a V + p 2 q v g a n - PaPpq^ ~ PaP^q v + q^.q v 



p ■ q _ p_ 

2 Pa 2 Q a 



+ h{Q 2 ) {q 2 g^ - q^qu) ( ^-p a - q a 



h{Q 2 ) 



1£ fivXaQ 



rPaP 



*\q 



0{ P ) 



(5.106) 
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Let us go back a little to our starting point. The whole dimension eight contributions that we 
are computing eventually show up in the chiral Lagrangian l|5.21(l as a next to leading order 
correction to the low energy coupling governing the AS = 2 transitions A| =2 . Knowledge of 
this coupling A^ =2 is what we eventually need. As already emphasised, a low energy coupling 
constant appearing in an effective Lagrangian can be related to the parameters of another 
Lagrangian through a matching condition. This can be accomplished once an appropriate 
Green function is chosen. When computing the dimension six contributions, i.e. Bk, we 
realized that the AS" = 2 low energy operator 



X 32 (D^)U X 32 (D^)U 



(5.107) 



contained a kinetic term in the right-handed external sources. We then wrote the matching 
condition as 




-QCD 




(5.108) 



which is equivalent to compute the contributions from 1|5.21[) and (|5.90() to the following Green 
function 



->S=2 
*a0 



(P) 



d A x e tp - x (0\T{Ri s (x)Rf(0)} |0) 



(5.109) 



where R d t s (x) = dn"f a SR(x). Straightforward calculations with the chiral Effective Lagrangian 
(|5~2T|) yield 

r 2 

= _i_LF f 4 A|_ 2 



G s a f( P ) 



PaP/3 
p2 



9a/3 



o( P 2 ) 



(5.110) 



whereas with the Effective Lagrangian in terms of quarks and gluons (|5.90(l one finds 



&P) 



eff 



G p 

3^ 



c 8 (/i) / d^xd^z (0\T{R d a %x)O 8 (z)R d /(0)} |0) 



-i ^Xl(4n^JdQ\Q 



2\l-e/2 



dn q T^ a (q,p)T^(q,p)(5.ni) 



where the first line above collects the short distances coming from quark integration and the 
second line the hadronic contribution coming from the up-up diagram. Equating (|5.110J) and 
(|5.111|l results in the following matching condition 



PaPfl 



- g a p ) A^ =s 



= ^c 8 ( M ) J d'xd'z j»* (0\T{R d J(x)O 8 (z)Rf(0)} |0) 
-^(4V) e/2 /^ 



2\l-e/2 



rf^r Aira (< Z ,p)L^' 3 (( ? ,p) (5.112) 



which is an expression for the dimension eight contribution to A$=2- Therefore, our problem 
from now onwards reduces to the determination of this low energy coupling, which, recall, is 
expressed as a matrix element of the dimension eight operators and an integral over the Eu- 
clidean regime of the square of the Green function T^ a . This is precisely how short distances, 
i.e., heavy quarks but also hadronic resonances down to the kaon scale, are encoded in low 
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energy couplings. The matching condition l|5.112[) is shown schematically in figure (|5.9|l . The 
convolution of the three-point Green function T ^ va with itself can be done without much effort 
to give 

(5.113) 



J dn q r, ua ( q ,p)T^( q ,p) = 



PaP/3 



~g af 3) w(Q 2 ) + o( P 2 ) 

where the tensorial structure is dictated by Lorentz symmetry. W(Q 2 ) is given by 



W(Q 2 ) = - lf(Q 2 ) Q 6 -- I 2 2 {Q 2 ) Q 2 + - ^ (5.114) 

Recall that, as a consequence of the symmetry properties of the tensors 7^„ a , there appear 
no crossed terms of the form factors in l|5. 114(1 . Once we have our function W(Q 2 ), the 
usual procedure is to find out its low energy and high energy expressions by means of Chiral 
Perturbation Theory and the Operator Product Expansion, respectively. Afterwards we shall 
connect them through an interpolator of the analytical form dictated by large- Nc QCD. A 



relatively simple analysis of the three-point Green function r^„ Q 
consist of simple and double poles, 

Ef <H bi, 

Nc 



reveals that it can only 



N c 



+ M 2 (Q 2 + Mf)2 



(5.115) 



as shown in figure (|5.10|) . Therefore, its square W(Q 2 ) has as much as quadruple poles. Its 
expression in the large- Nc limit looks like 



fW(Q 2 ) Nc 



E 



A, 



Bi 



C, 



D, 



Q? + M 2 (g 2 + M2)2 (Q2 + M 2)3 (Q2 + M 2 )4 



(5.116) 
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Rds 



L L L 



Rds 



Figure 5.10: Pole structure of the function 



in 



Eq. ]f>.115\) . "L" stands for L u d or L su . 



N C 



The above expression, though analytically simple, still requires an infinite amount of informa- 
tion. The strategy to truncate the above series is provided by the Minimal Hadronic Approx- 
imation. The number of resonances to be included in our analysis depends strongly on our 
knowledge of the high energy QCD behaviour of W{Q 2 ). At large momentum, for the form 
factors, one has 

i2a _ , /o ( 1 \ t lr a\ _ F o - ( 1 



W) = ^ + °{qe) • W) = -W + °W) (5 ' 117) 

It is rather difficult to get further in the OPE expansion. Knowledge of matrix elements such 
as 

(0\a L 'f'f'fG vp dz,\0) (5.118) 

is needed and, to the best of our knowledge, such matrix elements have not been determined. 
However, what we eventually need for our calculation is the function W(Q 2 ), for which the 
OPE gets simplified. The final result is 

QW)0Pi!S3 :| + !J( 1 + !) +0 (-l) (5.H9) 



11 



Recalling that 

M-V = -vope (5.120) 

we conclude that one can do with just one resonance in the interpolator (|5.116|l . At low energies, 
the contribution comes from the diagrams listed in figures (|5.11(l . (|5.12|) and Ij5.13|l . At low 
energies both momenta q and p are soft. However, taking into account the definition of our 
Green function T^ ai p has been sent to zero accordingly. At leading order (figure ()5.11(l ). the 
surviving diagrams are the first three ones. At 0(p 4 ), only Lg and Lio yield some contribution. 
Nonetheless, all terms containing Liq vanish when p — > 0. Actually, terms containing Liq can 
only appear in T1lr(p) at 0(p 2 ) and therefore consistently cancel as p — > 0. Even-parity 
contributions feed Hlr(0) and I\(Q 2 ). As for the Wess-Zumino-Witten term, its contribution 
is absorbed entirely in the odd-parity form factor ^(Q 2 )- This contribution is non- anomalous, 
since we already subtracted the anomaly. The presence of non-zero WZW contributions in 
non-anomalous processes is sometimes called the Cheshire cat smile of the WZW Lagrangian. 
At low energies, therefore, 

h(Q 2 ) = 2^+O(Q ) , I 2 (Q 2 ) = -J^ + 0(Q 2 ) (5.121) 

which, after squaring, yields 



l see the discussion in chapter 3. 
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Figure 5.11: Diagrams contributing to the leading order in the chiral expansion o/T p 



F n 4 - 



3 / iVc V 
2 ^24tt 2 J 



3Li 



Q 4 + 0(Q e ) 



(5.122) 



Thus, equations (|5.119|l and H5.122|) are the ones onto which we will match our interpolating 
function, to wit 



2 W{Q 2 )ha = a + 



B 



C 



D 



2 + A/2 (Q2 + M 2 )2 (Q2 +M 2)3 (Q2 + M 2 J4 

This will result in the following set of matching equations, 

7 
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(5.123) 



(5.124) 
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Figure 5.12: Diagrams contributing to the next to leading order even parity processes in the chiral 
expansion of r M „ a . 
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Figure 5.13: Diagrams contributing to F MJ/a coming from the WZW term. 



out of which we can extract the unknown residues A, B, C and D in terms of Lg, Fq, My 
and 5 K . Figure (5.14) shows the interpolating function (|5.123() for a given choice of input 
parameters together with its chiral and OPE extrapolations. Notice that the parameter a 
above, which accounts for the pion contribution, does not contribute to the integral of (|5.112[l 
(in dimensional regularisation). This is why we have subtracted it away in figure (5.14). 



5.2.7 Numerical analysis for AMk and ek- 

Using our expression for the Green function (|5.123|l and plugging it into (|5.112() we get 



A 2 

A S=2 



K 



dimS 



l0g Mf "3 



B 



c 



2My 



D 



(5.125) 



as our analytical result for the low energy AS = 2 coupling, once corrections coming from the 
next to leading order terms in the charm mass expansion are included. Recall that our result 
is scale and scheme independent, as it should. Indeed, the scale dependence coming from the 
integral of the A parameter in H5.124|) is exactly the one coming from the Wilson coefficient 
cs(fi) (see l|5.89|) V Taking as input parameters 



L q = 7 x 10" 



Fn = 0.087 MeV, 



0.77 GeV 



5 2 K = 0.12 ±0.07 GeV 2 , m c (m c ) = 1.3 GeV (5.126) 

we get 



A 2 

A S=2 



26 I 

T A C A U 5% + Xl (0.95 GeV) 2 \ (1 ± 0.3) (5.127) 

which shows the contribution from the two scales involved in the problem: 5k ~ 0.35 GeV, 
weighting the short distances and an hadronic scale A x ~ 1 GeV, responsible for the long 
distances coming from the up-up exchange box diagrams. The previous equation already 
shows that, since only A c has an imaginary part, contrary to A u , ex will be unaffected by the 
up-up diagram and will only pick a small correction coming from the short distance scale S K . 
However, the kaon mass difference picks the real part of A?j =2 and so gets corrected mainly 
by the hadronic scale A x . The equations relating our previously-determined A$=2 low energy 
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coupling to Ek and Am^ are as follows 



t/4 



G%FqMk 



6^tt 2 (Mk l — Mk s ) 



(5.128) 



ryiAfm^ + ??2 A* 2 (m?) e// + 2^^* to 2 log _| 
for whereas for the kaon mass difference, one obtains 



7 ^S=2Uim8 



Mk, - M, 



- Re<#°|«f;f (0)|i? 



3?r 



(5.129) 



Rc 



JJiAfme + '72A t * 2 {m 2 t ) eff + 2^X* c \* t m 2 c log 



A* 2 



S=2ldim8 



The first term in the curly brackets accounts for the dimension six contribution to A| =2 i 
whereas the second term is the result found in (|5.127|l . We stress that the previous results 
are in the chiral and large- Nc limits. Some chiral and 1/Nc corrections are known for the 
dimension-six part, which amounts to the replacement 



9S=2 







R ~ K 



(5.130) 



with Bk given in f5~£5|l and F K = 0.114 GcV. Lattice calculat ions of Bk with inclusion of 
chiral corrections push its value to Bk = 0.86 ± 0.15 [HZ]. Making the previous replacements, 
together with the input 



m t = 175 GcV, 



M K = 0.498 GcV 



(5.131) 



and the CKM matrix elements as given in jM], results in a 0.5% correction to Ek- As for the 
kaon mass difference, our determination yields 



M Kl ~ M Ks = (3.1 ± 0.5) x 10~ 15 GcV 
to be compared with the experimental value |98| 



M 



K L 



M 



K s 



(3.490 ±0.006) x 10~ 15 GeV 



(5.132) 



(5.133) 



cxp 



We therefore find that contributions arising from direct AS* = 2 transitions account for the 



(90 ± 15)% of the observed kaon mass difference. The remaining 



AS = 1 x AS = 1 



long 



distance contributions might account for the rest. We already know that they cancel at tree 
level due to the Gell-Mann-Okubo mass relation. However, calculations at one loop are stopped 
by the fact that there is no determination of the low energy couplings involved (See, for instance, 
the discussion in |14|). 



5.2 Determination of kaon mixing parameters 



<s() 




Figure 5.14-: Profile of the function Q 2 W(Q 2 )ha together with the chiral and OPE behaviour at low 
and high Q 2 (in D — A), respectively. The pion contribution has been subtracted away (the a parameter 
in Eq. \5.116\) ). according to the discussion in the main text. 



Chapter 6 



Duality Violations 



6.1 Motivation 



The OPE is one of the most important tools with which one can extract information about 
the QCD parameters and, obviously, the condensates. These parameters are however not 
unrelated, and, for instance, the condensates of the OPE of the two-point correlator H-lr{o 2 ) 
are related to kaon physics parameters 1 The OPE is known to be a valid expansion in the deep 
Euclidean, but physical information is only available on the Minkowski plane. Both regions 
can nonetheless be related with the use of Cauchy's integral theorem 



over the contour depicted in figure (6.1). Since knowledge of the full Green function is not 
available, the common practice is to replace the correlator by its OPE. Therefore, strictly 
speaking, we can rewrite Cauchy's theorem in the following fashion 



The substitution of the full Green function by its OPE is what goes under the name of duality. 
Differences arising from the above substitution, collected under 'D^(sq), are therefore coined 
duality violations 2 Since the OPE is believed to be a good approximation at large sq, there 
is hope that duality violations are small enough to be neglected. However, extraction of the 
condensates for Ulr has recently showed discrepancies between the different approaches |1()5| - 
|112) . pointing at the importance of gaining insight into these duality violations. 

Unfortunately, there is no theory behind duality violations and one has to resort to models. 
Our aim will be therefore to study duality violations in the two-point correlator function 
II£r(<? 2 ) with a model based on large- Nc QCD. We will start with a model in the strict large- 
Nc limit to later on move to a more realistic scenario by adding widths through leading 1 /Nc 
corrections. Our hope is to be able to single out the features which are model-independent and 
therefore can be extrapolated to the real QCD case. Our discussion essentially follows |104) . 

1 See, for instance, ll00l - IT02l . 
2 For a review, see 11031 . 
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Figure 6.1: Contour integral leading to 16', The crosses on the real axis stand for the resonance 
poles in large-Nc or the physical cut when resonance widths are taken into account. 



6.2 A toy model for duality violations 

Wc define the axial and vector vacuum polarization Green functions as 

n^(g) = * J Ae^(J^(x) (J^) f (0)) (6.3) 

with the QCD currents 

J£ = d(x) 7^ u(x) , J% = d(x) 7 V u(x) (6.4) 

We will hereafter work in the strict chiral limit. Lorentz invariance then fixes the tensorial 
structure of (|6.3J) to be 

n^(g) - (q„q v - 9,»q 2 ) H^ A (g 2 ) (6.5) 

Moreover, both II y and IX4 can be shown to obey a dispersion relation (up to one subtraction) 
of the form 

r°° dt 1 

KvAq 2 ) = / 5 -lmTl v , A {t) (6.6) 

Jo t~q z - it 7T 

Wc choose as our spectral functions the following ansdtze |71j 

hmll v (t) = 2F 2 p 5{t-M 2 p ) + 2Y j F*5{t-A4{n)) 

- 00 

-ImlU(i) = 2F 2 5(t) + 2S^ Fl6(t- Ml{nj) (6.7) 

7T * — ' 

n=0 

which describe axial and vectorial towers of narrow resonance states, whose decay constants 
Fv,A ( n ) and masses My, a (n) show an n-dependence to be specified. We take the p and 
7r analogs out of the tower to make our model more realistic. F p and m p are therefore the 
electromagnetic decay constant and mass pole of the p and Fq is the pion decay constant in 
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the chiral limit. Notice that (|6.7(l can be interpreted as an ansatz for the spectrum of large- N c 
QCD. The n-dependence of masses and couplings is then constrained and a possible choice is 



Fy A (n) = const. 



M v An) 



™v,a 



+ n A 2 



V.A 



(6.8) 



which means that our spectrum now displays a Regge-like behaviour of equally-spaced reso- 
nance states, much the same as what happens in QCD2 (N c — > oo) 3 . Furthermore, we choose 
F v = F A = F 2 , which will prove to be a convenient choice when we later on want to provide 
the resonances with a non-zero width. Matching our model to QCD short distances (see next 
page) will additionally lead to Ay = A^ = A. We will somewhat anticipate things and use 
hereafter a single decay constant and a single resonance spacing for both channels. 
We consider the following two-point Green function 



U v (q 2 )-U A (q 2 ) 



whose explicit expression, obtained from (|6.6|) and (|6.7(1 . reads 

^ F 2 



IW) = f 



F 1 



n=0 H 



M 2 (n) 



(6.9) 



(6.10) 



The infinite sums in l|6.10|l are by itself divergent and a regulator has to be defined to render 
them finite 4 . We introduce a cutoff A as our regulator. In order for chiral invariance to be 
preserved, the cutoff has to satisfy 



A 



N v K 2 V ~N A A 2 A 



V,A 



OO 



(6.11) 



As it was already emphasised in |11H| . this choice yields a well-defined answer, with infinities 
in each channel cancelling each other. Using the identity 



lim 



N- 



* — ' z + n * — ' n z 

(n=l n=l ) 



IE 



(6.12) 



where je is the Euler-Mascheroni constant and tp{ z ) the Digamma function, i.e., 



t l-e-t 

we can perform the sums in l|6.10|l explicitly to yield 

«2 



dt 



dz 



F 2 
zJL 

q 2 



-r 



M 2 p 



El 

A 2 



A 2 



iogr(z) 



A 2 



(6.13) 



(6.14) 



In order to ensure the right short distance behaviour of our model, we match l|6.14ll to the first 
OPE terms of QCD. This induces the following constraints: 



JP2 

A 2 



N c 
24tt 2 



(6.15) 



3 See the discussion in chapter 3. 
4 See discussion in 11131 . 
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to guarantee that the parton model logarithm is recovered; 



F* = F 2 



1 



A 2 



F 2 = F 2 



A 2 



and 



-IF 2 M 2 + F A 



m v m v J- i 2 2 

A 4 A 2 " 6 



mi m; 



A 4 A 2 6 



(6.16) 



(6.17) 



which are the constraints arising from the Weinberg-like sum rules, i.e., the absence of Q~ 2 
and (J~ 4 terms in the OPE. 

We are now in a position to evaluate the OPE as a l/Q 2 expansion for our model, defined 
as (Q 2 = -q 2 ) 



ngr (-Q 2 ) 



E C 2k 
CP 

fc=l,2,3.— ^ 



Using the derivative of l|6.13|l and the well-known expansion 

oo 1 

e z - 1 ^ w n! 



(6.18) 



(6.19) 



where B n (x) are the Bernoulli polynomials defined in Appendix C, the analytical expressions 
for the OPE coefficients then read, 



C2A: — — F 5k,i+ 

the first ones of which arc, explicitly 
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(6.20) 
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A naive 1 / Q 2 expansion of the spectral function in the dispersion relation 



where p{t) is a common short-hand for 



dt 



t + Q 2 -ie 



p(t) 



p(t) = - lmU LR (t) 

7T 



(6.21) 



(6.22) 



(6.23) 



6.2 A toy model for duality violations 



94 



would yield to the so-called spectral moments 



M n (s ) 



dt t n p{t) 



(6.24) 



which are not the coefficients of the OPE. To establish the actual expression for the OPE 
coefficients in terms of the spectral function, consider the Laplace transform of 



/>oo />C 

n(-<2 2 ) = dr e- TQ2 p(r) with p(r) = 

Jo Jo 

Plugging in the OPE expansion one gets 

C 2 k 



dt 



P(t) 



E 

k=l 



fc-1 



(A -1)1 



dt e' tr p{t) 



and finally, 



dt t k ~ l e- tT p(t) 



(6.25) 



(6.26) 



(6.27) 



showing that coincidence between spectral moments and OPE coefficients is accomplished if 
and only if the spectral function fall-off is of an exponential type. This is neither the case in 
QCD nor in our model. Actual computation yields 



M (s ) 

Mi (so) 
M 2 (s ) 

Ms (s ) 



C 2 -F 2 

-C 4 - F 2 
Ca-F 2 



1 



FA 



-C$ - F 
—F 2 A i 



Bi(x v ) - B 1 (x A ) 
Bi(x v ) - B^xa 
Bi(x v ) - Bi(x A ) 
B^ixv) - B 5 (x A , 
B\{x v ) - Bx(x A ) 



so + ^ 2 A 2 



B 2 (x v ) - B 2 {x A ) 



si + F 2 A 2 



B 2 {x v ) - B 2 (x A ) 



so 



■F 2 A U 



B 2 (x v ) - B 2 (x A ) 



B 3 (x v ) - B 3 (x A ) 



«0 + J™ 



B 4 (x v ) - B4(x A ) 



(6.28) 



where < xy, A < 1 is the fractional part of (so — rriy A )/ A 2 . Spectral moments show therefore 
an oscillatory behaviour around the OPE coefficient. This oscillatory behaviour is provided 
by the Bernoulli polynomial terms in l|6.28|) , in the form of a multistep function, which pushes 
upwards at vectorial mass thresholds (so = m^+nA 2 ) and downwards at axial mass thresholds 
(sq = m A + jiA 2 ). The expressions above are valid for so > my, i.e., above the threshold for 
axial and vectorial towers. The result is depicted in figure (6.2), with our set of free parameters 
chosen to be 

F a = 85.8 MeV , F p = 133.884 MeV , F = 143.758 MeV , 
M p = 0.767 GeV, m A = 1.182 GeV, m v = 1.49371 GeV , A = 1.2774 GeV (6.29) 



As it can be readily seen, duality points Sq, which would correspond to cuts in the sq axis 
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Figure 6.2: From left to right: Mo, Mi (first row), M2 and M3 (second row) as defined in the main 
text. 



of figure (6.2). arc absent 5 . Step-like discontinuities happen at mass pole thresholds for all 
the moments, but without any predictability: the OPE value is bounded by the maximum 
and minimum of the oscillations, but otherwise completely undetermined. Moreover, as the 
moment increases the oscillation gets amplified, because moments scale as M n ~ Sq. Even if 
no prediction is available, the qualitative information we get from this zero-width model hints 
at some general features to bear in mind. First, it provides a rationale for FESR, in the sense 
that duality points are common to all moments and exactly at mass poles. We will later on 
see that with the introduction of widths, this picture gets blurred: slopes become milder and 
therefore predictions are possible, but duality points get dissaligned. Second, the amplification 
of oscillations seems to select low duality points as the best candidates to make predictions for 
the OPE values. Again, the passage to adding widths will show that this picture gets distorted. 

The previous analysis concentrated on the comparison of moments and OPE coefficients. 
We already mentioned that duality violations are much bigger than the actual OPE values. A 
sound strategy would then be to reduce duality violations to a minimum. An approach in this 
direction is based on a slightly more sophisticated use of the Cauchy theorem 



/ dtw(t)-Imn LR {t) =—^-.<t dq 2 w(q 2 )Tl LR (q 2 ) 
Jo T 2m J\q*\=s 



(6.30) 



where w(t) is a polynomial chosen so as to suppress the duality violating oscillations where they 
are more harmful, i.e., on the Euclidean real axis. Such a strategy is usually coined pinched- 
weight sum rules (pFESR) I114j . Due to the recursiveness we observe in the oscillations 

5 Strictly speaking, duality points happen when Mo(sJ) = C2, Mi(sq) = C4 and so on. However, it happens 
that the OPE coefficients are much smaller than the typical size of the oscillations in figure (6.2), so that the 
duality points actually lie very close to the sq axis. 
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in Ij6.28(l . appropriate combinations of moments would be, by inspection, 

Mo-— M x ~ F 2 A 2 [B 2 (x v )-B 2 (x A )] 
so 

Mo --M 1 + \m 2 ~ F 2 A 4 [B 3 (x v )~B 3 (x A )] 
so s^ 

M - — Mi + \m 2 - 4m 3 ~ F 2 A 6 [B 4 {x v )-B 4 (x A )} (6.31) 
s s s 

which would correspond to pinched-wcights of the form 

w(t) = fl-1 



«;(*) 



so 
f 

so 



w(t) = (l-J-J ( 6 - 32 ) 

As one can see, such combinations indeed kill the dominant terms of the duality oscillations, 
leaving a remnants shown in (|6.31|) . Quite generally, a pinched- weight of the form 

w(t) = (l - (6.33) 

would involve the first M n +i moments, killing duality violations up to 0(F 2 A n B n+ \), Other 
pinched weights proposed in the literature are |114| 

t\ ( t^ 2 



Wl (t) = 1-3— 1 

so J V s 



w 2 (t) = - (l--) (6.34) 
so V s J 

which are not optimal in our model, as figure (6.3) shows, but still will be useful for later 
numerical comparison. The usual procedure in the pFESR approach is to make a fit over an 
So window between the pinched- weight, which is a linear combination of spectral moments, 
and the same combination of OPE coefficients to determine them all. The latter combination 
is obviously unable to reproduce any oscillatory behaviour (see figure (6.3)). Therefore, even 
though the oscillations have been atenuated, one expects the fit to be unreliable and very 
sensitive to the window chosen. However, notice that in figure (6.3) the OPE curve seems to 
interpolate between the pinched- weight oscillation. This is a consequence of the fact that B 3 (x) 
and B±(x) display a sinusoidal- like behaviour. This suggests the possibility of extracting the 
OPE coefficients by fitting over a full oscillation. Later on, with the introduction of widths, 
we will check whether this strategy is still effective. 

So far we have discussed two of the most employed methods in the literature. However, 
having a model of XIlb (<7 2 ) we can determine the exact form of duality violations in our model. 
Recalling 

f S ° dt tn P® dq2 ^ U LR E (q 2 ) +2? W (so) (6.35) 

J0 lm J\q*\=s 
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Figure 6.3: Integrals of the pinched weights, Wip,z{t) defined in the text (solid lines from left to right 
and top to bottom) and the corresponding OPE prediction (dashed lines). As it can be seen, ws(t) 
shows the smallest oscillations, i.e., the smallest duality violations. 



where T>^(sq) measures the duality violations. Remember that the OPE expansion is not 
convergent, at least on the Minkowski axis. Determination of the analytical form for the duality 
violations will be done with the help of the reflection property of the Digamma function 



1p{z) = Tp( — z) — 7T COt(7Tz) 



(6.36) 



On the negative half plane of figure (6.1) we can use the asymptotic expansion of the Digamma 
function at large \z\, to wit 



ip(z) ~ log 



(6.37) 



On the positive axis we can use (|6.36() such that we find an expression for Ulr (q 2 ) in terms 
of the OPE, to wit 



-27.VI/A' 



iWg 2 ) 



(u° L p R E {q 2 ) + o 

\n^(g 2 )+A^(g 2 )+0 



, Re(q 2 ) < 



(6.38) 



where the exponential terms appear as a consequence of the OPE being an asymptotic expan- 
sion. The duality violating function A 00 (q 2 ) appears only in the Minkowski region and takes 
the form 



A 2 



^2 1 



cot 



A 2 



cot 



A 2 



(6.39) 
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Figure 6.4: Contour over the complex q 2 plane used to split the duality violating term in the two pieces 
of FTJfy and tfTJUf) . 



This contribution does not vanish nor die out as q 2 increases, as it can be readily seen in figure 
(6.2). 

In view of 1|6.38|) and Ij6.35|l . the family of functions 2>M(so) can be expressed as a counter- 
clockwise integral over the semicircle \q 2 \ = sq (Re q 2 > 0) as follows 

2? W (so) = -^r J dq 2 q 2n A 00 (g 2 ) + 0(e- a ™"/ Aa ) (6.40) 

\q 2 \ = so 
Re{q 2 ) > 

It is specially convenient to split the previous integral into a spectral integral and a countertcrm, 
as shown pictorically in figure (6.4). Applying Cauchy's residue theorem, one sees that the 
duality violation can be expressed as two integrals, one over Re (q 2 ) and the other over Im {q 2 ), 

^"'(so) = 2?ii.(*o) +Z>£L.(«o) (6-41) 

with 

r l r s ° 1 

V oLu( s o) = / dtt n -ImA^t + ie) (6.42) 
Jo 

= ~^i\j_ is + j it \ dq2q2n A °» iq2) (6 - 43) 

where in (|6.42|) use has been made of Schwartz reflection theorem. T^} ciU (so) is the responsible 
for the oscillations appearing in the spectral function p(t) which the OPE is unable to capture. 
Mathematically, the singularities of the Digamma function are absorbed into the cotangent of 
(|6.36(l . which has the form 

1 °° 1 

ir cot (ttz) = - + 2z V — - (6.44) 

z z — ' z z — n z 

n—l 
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On the contrary, T>^j nst (so) gets contributions mainly from the low q 2 regime, decaying expo- 
nentially to a constant 

2>£L.(*o) = CN+ o(e- 2 ^/ A2 ) (6.45) 

whose explicit expression is 

C [n] = ^ - B n+1 (§) } - F 2 K - A 2 )" (6-46) 

As it should, plugging in the results of (|6.41() - (|6.46(l in (|6.35|l . one recovers the results for the 
moments M t of (|6.28f) . 



6.3 l/Nc corrections 

In this section we will move to a more realistic model departing from the zero-width model 
and allowing the resonances to have some width. In the light of large-Ac QCD, this would 
correspond to including the l/Nc corrections in a consistent way. Lacking a solution as to how 
to implement this corrections, we have to incorporate them in an ad hoc way. Obviously, this 
has to be done in a somehow sophisticated way if we want U LR (q 2 ) to preserve its analytic 
properties. For instance, the cut in the complex q 2 plane for Re q 2 is not achieved through a 
naive Breit-Wigner function. A convenient choice was found in |103l I115| , which, translated to 
our model, would yield 



U LR (q z ) = 11- 



irN r 



F 2 



F 2 
p 



oo 



F 2 



OO 



Ml 



- 1 z + M 2 A (n) ^ z + Ml (n) 



(6.47) 



where the variable z reads 



A 



2 • \ C 

-q z — is x 



A 2 



C = i 



irN r 



(6.48) 



As in the previous section, the sums can be performed and expressed in terms of Digamma 
functions, 



n LR ( q 2 








Z + ra A 
A 2 







z + m\ 
A 2 



(6.49) 



Notice that this expression provides a smooth transition to the zero- width case through letting 
the parameter a go to zero (or equivalently, £ to 1). This parameter a is proportional to the 
width of all the resonances. Indeed, resonance propagators aquire an imaginary part of the 
form 



r» = w M p (n) 



1 + 



(6.50) 
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Figure 6.5: Comparison of the data from ALEPH \116^ and OPAL (overlaid on a single curve) 

with the model, ilmll(g 2 ). As with the experimental data, the kinematical factors for tau decays are 
included and the pion contribution has been subtracted. 



meaning that resonance widths increase with the excitation number n. Therefore, the width 
of every resonance is modulated by the parameter a and by the way that masses grow with 
the excitation number (cf. I|fi.8[l ). The resulting behaviour is precisely the one found in the 't 
Hooft model. Comparison with experimental data suggests a value for a 



a = 0.72 



(6.51) 



while the other free parameters are the ones inherited from the zero-width scenario and are 
given in (|6.29() . The shape of the spectral function of Ulr (q 2 ) as compared with data is shown 
in figure (6.5). 

The analytical structure of HLR(q 2 ) has not changed except for the replacement q 2 — ► z. 



This means that TIlr(q ) has an expression very close to the one given in 



iWg 2 ) 



rng£*(*) + o 



,-2^/A 2 |< 



\n£^(z) + A(q 2 ) + o (^i^i 

where the duality violating term A(g 2 ) now reads 
nF 2 ] 



Re(q 2 ) < 
Re(q 2 ) > 



namely 

(6.52) 



A(0 



A 2 C 



— < cot 



-q 

A 2 



2\ C 



A 2 



cot 



,2\ < 



A 2 



A 2 



(6.53) 



with £ as given in (|6.48|l . The OPE has an akin expansion to the one in l|6.20(l , but the variable 
z, once expanded in terms of q 2 , gives rise to logarithmic corrections in the OPE coefficients, 
to wit 

a ka , Q 2 . m ( a 2 \ " 



C 2k {Q z ) 



2\N C =3 



a 



i 



,n, 1o ^ + {n 2 ) 



(6.54) 



where C^ -00 are the zero- width OPE coefficients as given in (|6.20() . These logarithmic cor- 
rections mimic those found in real QCD. Therefore, it will prove convenient to split the OPE 
coefficients as follows 



c 2fe (Q 2 ) 



2\N C =3 



= a 2k + b 2k log + O , ^ 



(6.55) 
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Figure 6.6: From left to right: Mo, Mi (solid curves, first row), M2 and M3 (solid curves, second row) 
as a function of so together with the OPE prediction (dashed curves) as it stands in the square brackets 
of jttfy ) to PHfy . 



We must note, however, that &2fc coefficients, which yield the anomalous dimensions in QCD, 
do not arise there as a 1/Nq effect, unlike in our model. With this caveat in mind, we can plug 
in our set of parameters I|6.29I6.51[) to obtain the numerical values for the OPE coefficients 

0-2 = b-2 = = CL4 = 64 

a 6 = -2.8 x 10~ 3 GeV 6 , b 6 = -5.9 x 10~ 4 GeV 6 

a$ = +1.8 x 10~ 3 GeV 8 , b s = +5.1 x 10~ 4 GeV 8 (6.56) 

Contrary to what happens with the OPE coefficients, i.e., in the Euclidean, where we have 
seen that l/Nc corrections are rather mild, the spectral function in the Minkowski changes 
dramatically with the inclusion of widths. This follows from the fact that the Nc — * 00 and 
q 2 — * 00 limits do not commute, as already pointed out in |10H| . Indeed, in the Minkowski 
region, (q 2 =t>Q), the variable z takes the form 



t 



,2 



z= -\ 1 -^{^) +i N: +0 {N!)) (6 - 57) 

Taking the imaginary part of n(g 2 ) in the Minkowski, one finds 

ImII(i) = Imn OPE (i) + ImA(t) (6.58) 

where the term coming from the OPE is precisely the analytical continuation to the Minkowski 
of the OPE in the Euclidean, 
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Thus, the absence of t -1 and t~ 2 terms above is closely related to the absence of dimension-two 
and dimension-four operators in the OPE, a condition we have enforced with the two Weinberg 
sum rules. The last term in H6.58|) reads 

ImA(t) = 

AirF 2 _2 2L at . / 2t-m 2 A ~m\\ . / m\ - m\\ I ' _ / a 
— e ^ sm^ m ^_^j(i + (_ 

+ o(e~%?&} (6.60) 

and stems from the duality violating piece A(q 2 ). It is illustrative at this point to see how 
the picture has changed for the moments as one has gone from the Nc — > oo case of figure 
(6.2) to the finite Nc case depicted in figure (6.6). Expression l|6.60|l illustrates our previous 
statement of non-commutativity of the Nc — > oo and t — > oo limits. Indeed, taking Nc — > oo 
first leads to the delta function spectrum of the large- Nc limit, whereas reverting the order 
yields an exponential damping which eventually dies out (see figure (6.6)). 

This dramatic change in the imaginary part of A(g 2 ) also translates into its real part. 
Actually, the full function A(q 2 ) behaves as follows 

-2n( lQ | I sin(t((p— tt))| „ „ ((]<(£><- 

^%7^ e [ } ' ^ = k 2 |e-, {L-<V< 2 27r (6-61) 

Recall that the Nc — > oo limit (i.e., £ — * 1) in the previous equation kills the exponential 
fall-off for = 0, 2tt, but it is exponentially suppressed otherwise. This is in accordance with 
what we observed previously in figure (6.2): on the Minkowski axis oscillations are undamped. 
Therefore, in the finite Nc case, contrary to what happens in the Nc ~* oo case, the duality 
function T>^(sq) displays an exponential suppression everywhere 

£> W (s ) ~e-^% (6.62) 

which allows one to write 

2? W (so) = ^LiM+^lM ~ O(e-W) (6.63) 

V [ :J alL ( So ) = J°°dtt n hmA(t + ie) -CM+o(e-M) (6.64) 

^Lt.(so) = C [ " ] + 0(e- 2 ^) (6.65) 

Therefore, combining the previous equations and taking the large sq limit results in the follow- 
ing 

f°° 1 

C [n] = - dt t n -ImA(t + ie) (6.66) 
Jo * 

which also means that T>^(sq) can finally be cast in the simple form 6 

dt t n -ImA(t + ie) + O [e'^j (6.68) 



All along we have used that 
which can be shown to be satisfied from a relatively low sq upwards 



27ra sp o SQ 

e JVc > e a* (6.67) 
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Armed with this expression for the duality violating term, we can write the analog of H6.35J1 
for the finite Nc version of our model in the following fashion 



dt t 2 p(t) 



dt t 7 p(t) 



/■so 

/ dt t 3 p(t) 
Jo 



aw + he log 
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/ dt p(t) + 
Jo 
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bs 
2s 2 
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a 6 + 6 6 log - 


h — 
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F&8 




^6 So + 


ho 
so 



A 2 



be4 



b&4 



bis 
so 



+ 



= ^ [0] (s ) 


(6.69) 


= V [1] (S ) 


(6.70) 


= ^ [2l (s ) 


(6.71) 


= 2? [3] (so) 


(6.72) 


= ^ [7] (so) 


(6.73) 



where the OPE integrals (the terms in brackets above) get contributions from the hk param- 
eters, according to the following formula 



^ f dq 2 (g 2 

2m J\q2\= S0 



2\n T\OPE 



LB 



(q 2 ) = (-i) r 



O-ln+2 



->2n+2 



1 S ° 

l0g A^ 



OO 

E 



(-if 



J2k+2 (n 
" S 



-fc) 



k\ 



(6.74) 

The extraction of the a,2k OPE coefficients therefore gets polluted not only from the duality 
violation term but also from all &2fc terms. The existing analyses for the extraction of OPE 
coefficients in QCD do not include the influence of these hk terms. Certainly this seems 
reasonable for the hk which come with inverse powers of so> even though one should be aware 
of their growth with increasing dimension (eventually, they get bigger than the a,2k parameters) 



>2k 



N c 2k 



(6.75) 



What certainly deserves more attention are the hk terms coming with positive powers of sq 
out front. Despite their small absolute value, at so = m 2 . ~ 3.15 GeV 2 the so term alone can 
bring a huge number. In the particular case of the dimension-six condensate, Q6.71(l shows that 
there are no harmful hk terms coming with a positive power of sq. Moreover, in real QCD 
(only) &6 has been determined so far and shown to be negligible. This might very well explain 
the overall agreement in o,q of all the existing determinations 7 . Concerning the dimension- 
eight condensate, the discrepancies found in the different analyses on the market may signal 
at sizeable hk terms not taken into account. 

As already pointed out, the hk terms in our model do not have the right l/Nc scaling. 
This results in a value for the &6 parameter which is far bigger than the observed one. Probably 
the same happens with the remaining bik terms, so that our model would be overexaggerating 
the weight of the hk terms. Should this be true, they would have a milder impact in QCD 
than they have in the model. Anyway, at least for consistency, one should try to assess their 
actual impact. 

7 See, for instance, E||7H], ITP51 - ITT21 . ITT5I . 
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We now move to the analysis of the duality violating terms, which in our opinion bear 
the main source of difficulty in the determination of the OPE condensates. To show that 
point clear, we will compare a FESR and pFESR analysis without taking into account duality 
violations (which is what present determinations do) to later on show how to improve on that. 

6.3.1 Finite Energy Sum Rules 

The strategy to follow will be to identify the duality points Sq in the moments Mq and Mi, 
i.e., those points satisfying M (sq) = and Mi (sq) = 0. Present experimental data cover 
the first two duality points, which we will refer as and Sq 2 \ Contrary to what happens 
in the large- Nc version of our model, duality points need not be the same for the different 
moments, and actually they are not. However, and this is the main motivation of the finite 
energy sum rule approach 8 , they should lie close enough to allow for a reliable extraction of 
the OPE coefficients. In the light of lj6~oT)j) . finite energy sum rules take the following scenario 

M (sS) = = Mi(aS), M 2 (s%) ~ a 6 , M 3 (s* ) ~ a 8 (6.76) 

as a good enough approximation to the physical picture. The first duality point in our model 
sits at 4 1} = 1-472 GeV 2 , and one obtains the following predictions 

Al BSR = -4.9 • 10" 3 GeV 6 , A FESR = +9.3 • 10~ 3 GeV 8 (6.77) 

which are to be contrasted with the OPE expressions, which have to be truncated at some 
point. We choose 

to yield 

A 6 (1.47 GeV 2 ) = -2.4 • 10~ 3 GeV 6 , A 8 (1.47 GeV 2 ) = +2.6 • 10~ 3 GeV 8 (6.79) 

Contrary to what happens to Sq 1 ' , the second duality point does not simultaneously satisfy 
Mo(sq) = = Mi(sq). This is quite similar to what happens in real QCD. Following |108| . 
we choose Mi(sq ) = to define the duality point with which to make our predictions. This 

corresponds to Sq 2 ^ = 2.363 GeV 2 . This strategy seems to be well justified, since one expects 
duality points to shift minimally between neighbouring moments. The results we find are 

j^fesr _ 2 • 10~ 3 GeV 6 , A FESR = —1.6 • 10~ 3 GeV 8 (6.80) 

to be compared to the true OPE values 

A 6 (2.36 GeV 2 ) = -2.8 • 10~ 3 GeV 6 , A 8 (2.36 GeV 2 ) = +3.4 • 10~ 3 GeV 8 (6.81) 

In view of the preceding results, one can hardly argue in favour of a preferred duality point. 
In principle, one should expect that a larger s$ would reduce the amount of duality violation, 

8 For a review, see 1911 . 



ae + b 6 log -- + — , 

s* 
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(6.78) 
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Figure 6.7: Comparison of the pinched weights wi(t), ui2(t) and ui:j(t) as defined in the main text with 
their corresponding OPE contributions subtracted. Vertical axis in arbitrary units. 

thus improving the predictions. Indeed, the prediction for A§ is slightly better in the second 
duality point, but clearly this is not the case for A%. Most probably, the amount of duality 
violation is still big enough to ruin this picture. Remarkably, the scenario we find in going to 
the first duality point to the second one, i.e., Aq reducing to half its value and As changing sign 
mimics the one observed in real QCD. However, regardless of the duality point, uncertainties 
in the predictions amount to, at least, a 30%. In the large-iVc version of our model we saw 
that pinched weights indeed suppress the duality violating terms, but do not get rid of them 
entirely. We here reproduce an analysis in the same spirit as that of |114j . We fit the OPE 
contributions to the pinched- weights wi (t) and W2 (t) , choosing a window of equally-spaced 20 
points ranging between 1.5 GcV 2 < sq < 3.5 GeV 2 . The upper end is chosen so as to coincide 
with the would-be m T mass threshold 9 . The values obtained this way are 

Al mch = -3.8 x 10~ 3 GeV 6 , A p 8 mch = 6.5 x 10~ 3 GeV 8 (6.82) 

which also fail to reproduce the true values better than a 30%. 

6.3.2 Other approaches 

We can also use our model to test the minimal hadronic approximation we have been using in 
previous chapters as an interpolating function for different Green functions. As a starting point , 
we will restrict our spectrum to consist of a pion, a p resonance and an axial a\ resonance, all 

^Obviously we have not included the tau lepton in our model. In real QCD it falls around the p(1700) 
resonance, which in our model would correspond to rn? n ~ 3.8 GeV 2 . 
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of them infinitely narrow, namely 

= -F 2 S(t) + F 2 6(t My) - F 2 5{t ~ M\) 



1. 

7T 



-i m nf M ' 



It then follows that 
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F 2 
1 n 



(i 2 ) = ^ 



F 2 



M v - q 2 



M\ - q 2 



(6.83) 



(6.84) 



with free parameters Fq, Fy, Fa, My and Ma- We fix F$ to take its value from our model of 
previous sections, i.e., 

F = 85.8 MeV (6.85) 

and determine the remaining four parameters by demanding the two Weinberg sum rules to 
hold 

F 2 - F 2 A = F 2 , F v My-F 2 A M 2 A = {) (6.86) 

Additionally, we impose constraints on the slope at the origin and area below TlLR{q 2 ) hi the 
Euclidean. These are given by Lio and the electromagnetic pion mass difference, defined as 



Lio — 



d 

dj 2 
3q 



ml 



8ir 2 m n F^ J 



q ' 2 U LR (q 2 ) 
dQ 2 



q 2 ^Q 

Q 2 U LR (Q 2 ) 



(Q 2 



(6.87) 



The resulting matching equations for the previous two parameters as given by the model of 
(|6.47(l and the Minimal Hadronic Approximation read 



l^M 2 



L w - -7-P0' 



3a 



M 2 



M 2 M\ 



-(5.2 ±0.5) • 10" 



8TT 2 m^F 2 



f; 



M\M 



A 1V± V 



Ml- My 



/ M 2 \ 

log \ = (4.2 ± 0.4) • 10" 3 GcV (6. 

\ M vJ 



where we have taken a = 1/137 and m T = 135 MeV. Combining l|6.86|l and (|6.88() yields 
My = 0.70 ± 0.01 MeV 



F\ 



122 ± 6 MeV 



M A = 1.00 ±0.03 GeV 
F A = 84 ± 7 MeV 



(6.89) 



Plugging l|6.89|) and (|6.85|) into (|6.84|) . we obtain the following predictions for the OPE coef- 
ficients 



hi HA 



= -(3.6 ±0.3) • 10~ 3 GeV 6 



,MHA 



(5.4 ±0.7) • 10~ 3 GeV 8 



(6.90) 



which again are comparable to the previous estimates using FESR and pFESR. 

Just for illustrative purposes, figure (6.8) shows TlLR{q 2 ) of our model as compared to 
U L J R HA (q 2 ). We here reproduce the well-known fact that the MHA to Green functions, despite 
being extremely different from the experimentally known QCD Green function in the Minkowski 
region, is a very good approximation in the Euclidean. Remember that the Euclidean region 
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Figure 6.8: Comparison ofHi,R(q 2 ) 



as it stands in \6. H!ifl with the experimental data of ALEPH 



lllfif and OPAL \117^ . The kinematical factors for t decay are included and the pion pole subtracted 
away. As already emphasised, being a bad approximation to the Minkowski, the MHA is however a 
very good approximation to the model in the Euclidean region. 



is precisely where we have been using the MHA as an interpolator to compute nonlcptonic 
electroweak couplings in the previous chapter. Recall that electroweak couplings have the 
geometrical interpretation of being the area under the curve, and indeed the agreement is 
remarkable. 

Another commonly used method to extract the OPE coefficients relics on the Laplace 
transform of the Cauchy theorem. They are usually coined Laplace sum rules or sometimes 
Borel sum rules. Following |118| , we can write 

r 2 6 f s ° 2 f s ° 

«6 -ofo 777 + - = — / dte- tT p(t) + - dtte~ tT p{t) 
iZ T Jo T Jo 



12 



4 + aw 77 + ••■ = -^f / dt e ~ tT P(t) ~ % I dtt e ~ tT P(t) ( 6 - 91 ) 
2 ' Jo T Jo 

which yield the dimension six cig and dimension eight condensates a§ in terms of the (expo- 
nentially weighted) moments. The strategy in Laplace sum rules is to look for a stability region 
of the parameters ctg and ag in the Borel parameter (r) space. We follow |118j and fix so to 
be close to a duality point. As it can be seen in figure (6.9), no such stability region is found. 
We have also moved the so and nowhere stability was found. Laplace sum rules therefore do 
not work in our model, even though in real QCD they seem to do. 

6.3.3 Modelling duality violations 

The previous analyses suggest that indeed the inclusion of duality violations is necessary in 
order to improve our determinations beyond 30%. The challenge is then how to model duality 
violations in a generic way. Recalling (|6.60(l . we take as our ansatz 

-ImA(t) = k e -7 * sin {a + (it) (6.92) 

One possible strategy is to use (|6.92|) to fit the right-hand and left-hand sides of H6.69J) . (|6.70|) . 

where we set the &2fc coefficients to zero. This is obviously an approximation, but it seems 
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Figure 6.9: Plot of the coefficients a\ (left) and <2g (right) for a fixed so close to a duality point, 
together with the true values d6,8 from Eqs. 16. 561) . as a function of the Laplace variable, r. 



reasonable as all those coefficients come with inverse powers of sq. The fit results in the 
following values for the parameters 

K= 0.026, 7 = 0.591 GeV~ 2 , a = 3.323, /3 = 3.112 GeV~ 2 (6.93) 

which produce the first two plots of figure |6~Tn)l . The prediction for M 2 and M 3 is given m 
the last two plots, where one can see that the agreement between the true and predicted curves 
is remarkable, improving in the intermediate regions where Mi and M3 cancel. Once we have 
modelled the duality violations, we can use the predicted curves to find the duality points for 
each moment. For the sake of illustration, we show how this strategy works with Mi, M3 and 
M7. Predicted duality points sit at 

V [2] {s* ) = , s* = 2.350 GeV 2 

V [3] (s*)=0 , s* = 2.307 GcV 2 

V [7] (s* o ) = , Sq = 2.130 GeV 2 (6.94) 
Therefore, ijfiTTJl . (jeTY!^ and ijeTSj) give 

a 6 + b 6 log % + — - = [ ° dt t 2 p(t) = -0.00251 GeV 6 (6.95) 
A s o 2s o Jo 

a 8 + b s log i| - b 6 s* + ^ = - [ dt t 3 p(t) = 0.00329 GcV 8 (6.96) 

. 6lB log i| - ^ + h -^f- b u s* + ^ = - ['dt t 7 p(t) = 0.0179 GeV 16 (6.97) 

A 5 4 s n Ja 



aw 



whereas the true value on the left hand sides should be -0.00281 GeV 6 , 0.00344 GeV 8 and 
0.0161 GeV 16 , respectively, meaning a 10%, 4% and 11% error, clearly much precise determi- 
nations than those obtained by any method neglecting duality violations. 
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Figure 6.10: Upper row: Plot of the left-hand side (solid line) and right-hand side (dashed line) of 
Eqs. l6.69\6.7Q l, neglecting the b coefficients, as a function of so (in GeV 2 ). Lower row: Plot of the 
left-hand side (solid line) and right-hand side (dashed line) of Eqs. \6. 71\6. 7ty , as a function of so (in 
GeV 2 ). The right-hand sides are obtained with the fitted parameters of \6.9S\) . 



Chapter 7 

Conclusions 



In this work we have exploited some of the features of an approximation to QCD, called large- 
Nc QCD, in order to provide physical insight into the non-pcrturbativc aspects of the hadronic 
world. In chapter 3 we introduced a Lagrangian designed to go one step further than chiral 
Lagrangians in the sense of being capable of describing the interactions between the pseudo- 
Goldstone bosons (the pion octet fields) and the lowest lying resonance states. The Lagrangian 
admitted an interpretation in the framework of large- Nc QCD. This turned out to be a powerful 
tool which allowed a well-defined computation of quantum corrections, in powers of 1/Nc- We 
therefore tested the validity of the Lagrangian with a simple example, the prediction for Liq at 
the quantum level. The tree level determination of that parameter with the same Lagrangian 
was long ago shown to be in good agreement with experimental data. We however found that 
this agreement no longer persists at the quantum level and the Lagrangian is unable to yield 
the right evolution for Lio(^t), dictated by chiral perturbation theory. We suggest this as one 
of the easiest non-trivial exercises to put to the test every Lagrangian purported to go beyond 
chiral perturbation theory. 

In subsequent chapters we set the stage for the study of kaon physics parameters. In 
particular, we stressed the fact that in nonlcptonic weak interactions scale-independence in the 
hnal physical results is not easy to achieve. This has to do with the fact that matching long and 
short distances is highly non-trivial: Wilson coefficients encode the ultraviolet aspects of the 
calculation, i.e., the short distances, and can be computed with perturbative techniques; as for 
the long distances, one has first to find the chiral realization of the operators involved. It can 
be shown that the complicated part of the calculation thereafter amounts to a determination 
of the (associated) low energy chiral couplings, which can always be expressed as integrals over 
certain Green functions. Non-perturbative tools are then needed to further proceed. Yet they 
have to make sure that these two pieces, Wilson coefficients and hadronic matrix elements, are 
smoothly connected in order to make any spurious scale dependence disappear from the final 
physical answer. We showed how a well-defined approximation to large- Nc QCD, termed the 
Minimal Hadronic Approximation, solves in a very natural way the problem. 

The advantage of the Minimal Hadronic Approximation is that it provides an analytic ex- 
pression for those Green functions lying under the low energy chiral couplings, expressions 
which are valid at all energies. Those expressions can be thought as interpolating functions 
(rational approximants) between low energies and high energies, where calculational methods 
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exist (%PT and the OPE). In particular, we have shown explicitly how the correct imple- 
mentation of the OPE constraints is the feature that ensures the right matching and renders 
scale- independent results. 

We have applied this method to the determination of the parameters governing kaon phe- 
nomenology: we have provided a new value for Bk in the chiral limit and to next to leading 
order in the 1 /N c expansion, taking into account the effect of including dimension eight oper- 
ators in the OPE of the underlying Green function. This was motivated by a recent analysis 
concluding that these dimension eight operators were potentially sizeable, therefore suggesting 
its inclusion. Our result shows that Bk is essentially unaffected by such terms. 

Another issue was the evaluation of the corrections in the charm mass expansion to the 
kaon mass difference, Amj-, and ex- The charm mass expansion is known to have slower 
convergence due to its proximity to A x and therefore it is a source of potentially sizeable 
contributions which deserve attention. This calculation in inverse powers of the charm mass 
is tantamount to computing the leading short and long distance corrections to the AS* = 2 
Lagrangian. Such corrections give rise to two different scales: a long distance hadronic scale 
A x ~ 1 GeV due to the (non-perturbative) up-up exchange box diagram and a short distance 
parameter 5k ~ 0.35 GeV. It turns out that A x only has impact on Am^, while Ek receives 
a small 0.5% correction proportional to 8k- The net result is that, whereas Atuk has a 
sizeable correction of a roughly 10%, which pushes its value closer to the experimental one, 
£k remains basically unaffected. It is worth stressing here that our results were performed in 
the chiral limit, and thus there is room for potential improvement with the addition of chiral 
corrections. Furthermore, recall that we have only computed the short distance contribution 
to the kaon mass difference. There still remains an evaluation of the so-called long distance 



AS = 1 x AS = 1 



contributions, which lies beyond the scope of the present work. 



In the last chapter we have addressed the issue of quark-hadron duality violations, which 
arise as a consequence of the lack of validity of the OPE over the whole complex q 2 plane, 
most blatantly in the Minkowski axis. We have focussed our attention on the two-point Green 
function YlLn(q 2 ), whose OPE coefficients are related to kaon physics parameters. Armed with 
a model, we have tested the most relevant approaches to the extraction of the OPE coefficients 
and have shown how in principle one can improve on them. Our conclusion is that they all 
lie in the same ballpark of values and to really make an improvement, at least in our model, 
one has to take duality violations into account. We have put forward a proposal to model the 
duality violating pieces, which in our opinion goes beyond our particular toy model and can be 
considered as generic. As such, we suggest that it be incorporated in QCD analyses of ^lr{<1 2 ) 
in order to reassess the existing approaches used to determine the OPE coefficients. 



Appendix A 

SU(3) algebra 



A.l SU(3) algebra 



The eight generators of the algebra are usually written down in terms of the following matrices, 
known as the Gell-Mann matrices, 

1 \ / -i \ / 1 

Ai = | 1 , A 2 = \ i , As = -1 
000/ \ / \ 

1 \ /0 -i\ / o o o 

000 , A 5 = , A 6 = 1 
100/ \ i / \ 1 

000\ 1 ( 1 \ 

A 7 = | -i \ , A 8 = -= 1 (A.l) 
i / o -2 / 

It is relatively easy to check that they satisfy the following algebra 

[A , Ab] = 2i / Qbc A c (A. 2) 

{A a ,A c } = -^5 ab l + 2d abc \ c (A.3) 

JVC 

which leads to 

2 

A a A b = — — S ab l + d abc X c + if abc A c (A. 4) 

i\c 

from where it is immediate to verify that 

Tr(A Q A h ) = 2 ( 5 ab (A.5) 
Another useful relation is the following 

t) w (^)--K**"-W ( "» 



112 



A. 2 Fierz transformations 



113 



The values of the (completely antisymmetric) structure constants f a b c appearing in (jA.2IA.4p 
can be inferred from their expressions in terms of the Gell-Mann matrices 

fabc = -r. tr[A a ,A b ] A c (A.7) 

to yield 

/123 = 1 

/l47 = — /l56 = /246 = /257 = hib = ~ fi67 = ~ 



/458 — /678 — (A.' 



while the (completely symmetric) coefficients d a bc arc defined as 



dabc = |tr{A a , A fc } A c (A.9) 

whose values are listed below 

^118 = <^228 = ^338 = — ^888 = 

dl46 = ^157 = ~^247 = ^256 = ^344 = ^355 = — <^366 = ~ <^377 = 7^ 

C?448 = <^558 = ^668 = ^778 = (A. 10) 

As for the ST/ (2) group, their generators are usually written in terms of the Pauli matrices, to 
wit 

1 \ f -i \ ( 1 



i o J ' T2 = [ i o ) ' T3 - 

which satisfy the SU(2) algebra 

[Ti,Tj] = 2i£ij k T k (A. 11) 

A. 2 Fierz transformations 

Consider a product of Dirac bilinears of the form 

^r^i^rVa (A.12) 
where T is to be understood as the most general matrix in Dirac space, 

r = c s T s + c v IV + c T T T + caTa + cpTp (A.13) 

where 

i 

= 757^ i r P = 75 ( A -14) 
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Sometimes one may want to express l|A.12(l as 

i' 3 r fa i'4 r ip 2 

This amounts to a change of basis. Defining the Dirac bilinears as 



(A.15) 



O s (i,j;k,l) 
O v (i,j;k,l) 
T (i,j;k,l) 
A (i,j;k,l) 
P (i,j;k,l) 



one can write 



2J c M (T M )ij (T M ) M (D M (i,j;k,l) = 2J 



with the coefficients cm given by 
where 



cm — Aa/jv 



A 



MN = 77 



/ 


2 


8 


24 


8 


2 


\ 




2 


-4 





-4 


-2 






1 





-4 





1 






-2 
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-4 


2 
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-8 


24 
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/ 



(A.16) 

cat (T N ) k j N (i, l; k,j) 

(A.17) 
(A.18) 

(A.19) 



is the matrix associated to the change of basis. 



Appendix B 

Schwinger operator formalism 



This formalism |85j is a very convenient method to do calculations involving covariant deriva- 
tives because it preserves gauge invariancc at all stages of the calculation. This is unlike 
ordinary diagrammatic perturbation theory where the c* M and the g s G^ sitting in a covari- 
ant derivative appear at different orders in the expansion. The advantages provided by this 
technique are difficult to emphasize, since working in coordinate space in a covariant way is 
especially suited for extracting without much effort the operators and coefficients that build 
up the OPE. Unlike a Feynman diagram approach, there is no need to make any surgery: the 
formalism singles out the Wilson coefficient and its associated operator in a most natural way. 

B.l The background field method 

We start from the QCD Lagrangian 

£qcd = -\ G$ Gft + q (i p - m) q (B.l) 

where 

G$ = -d v G^ +g s f bc GfG^ (B.2) 
Gauge fields are split in the usual way in any background field method 

G,{x) = G,{x) + g^x) , 9fl (x) = ^ (B.3) 

where G^(x) is the classical field, which will play the role of a static background upon which 
the quantum g^(x) will propagate. The field strength tensor can then be rewritten as 

G$ = G$ + - D v gl + g s r hc g\gl (B.4) 

where G$ is the classical field strength 

G$ = D ^i a) - D V G^ + g s f abc G^G^ (B.5) 
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and the covariant derivative is redefined to contain only the (strong) background field 



0^ = 8^- ig 3 Gf_, 



(B.6) 



It is worth emphasising that the background field in the equation above is of gluonic nature 
because we are interested in the strong corrections. Other interactions can in principle be 
accounted for once the proper background fields are provided. The QCD Lagrangian with the 
splitting (|B.3() can be written 



Cqcd = -\(g$ + - D v gV> + g s f abc g<A + q{ip + g s 7 %)g (B.7) 

B.2 The Schwinger formalism 

One starts by defining quantum mechanical operators X^, P v satisfying the eigenvalue equation 

X„\x) = x„\x) , P„\x) = iD„\x) (B.8) 
from where it is immediate to compute the usual commutation relations 



Xft, Py 



P P 



~2 



One then has that 



(x\P\y)=il& 6{x-y) 



(B.9) 



(B.10) 



where the covariant derivative is understood as acting upon the delta function. With this 
definition the quark propagator can be expressed as 



Si(k)= I d 4 xe lk - x (x\^ 10)= / d 4 x(x\ 



+ ft —m. 



|0) 



(B.11) 



f - m t 

where the second equality follows from the relation 

e ik-xp v = (p v + ku yk-x (B 12) 

Analogously, the gluon propagator can be computed from (|B.7|) by conventional means to yield 



-I ab 



\y) 



(B.13) 



_P 2 g lu ,-2g.f abB G> a , J 

Following (|B.8|) . any quark field q(x) is to be understood as q(X), so that one has 

q(X)\y) = q(y)\y) (B.14) 
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B.3 Dimension six and dimension eight operators in the 
weak OPE 

With these definitions we can apply the formalism to AS = 2 transitions. As an example of how 
the method works, we will derive the leading and next to leading order operators in the OPE of 
the function Wlrlr(q 2 ) as defined in section 5.2. All along we shall regulate divergent integrals 
in dimensional regularization using dimensional reduction |119j . In this regularization scheme 
one keeps the algebra of Dirac matrices in four dimensions while the integrals over momenta 
are performed in d = 4 — e dimensions. Consider the four quark operator 



s L {x)^d L (x) s L (0)7"d L (0) 



(B.15) 



where the quark fields s(x) and d(x) are to be understood as the full ones, in contrast with the 
perturbative ones s^°\x), d(°\x). Their relation is the following 



sl(x) 
d L {x) 



sf{x)+ig a / ctxs^>{x)G^x)rS{x-x) 



j4~ -(0)/- 



4 0) O) + ig s I d 4 x S(x - x) G„(x)^ d { ° } (i) 



(B.16) 



Plugging the previous expressions in (|B.15J| leads to the following expansion 



Mi 



{x)^df{x) ^(0)7^d^(0) + 



(0)/ 



; (0) 



(x)Y 



d"x S(x - x) ig.G^x) ^ d { ° ] (x) 



(0)7* 



d?yS(0-y)ig a Gp(y)<fd$ ) (3j) 



d^x S(x — x) ig s G^(x) 7^ df (x) 



d A x sf^ (x) ig s G fJl (x) 7 M S (x — x) 



d^x (x) igsG^x) 7^ S(x — x) 



l,df{x) 



d 4 ys { °\y)ig s G^y)rS(0-y) 



7^l } (0) 



7m4 0) (0) 



^df(x)s^(0)r 

+ 0(ai) 



d 4 yS(0-y)ig s G,(y)rd^(y) 



(B.l 



The first term is the non-interacting one, whereas the following four lines account for the one- 
gluon exchange contribution, depicted in figure (5.3). Remember that G^ can be factorized 
as a background added to a dynamical gluon field. Plugging in the expressions for the quark 



B.lll) and gluon l|B.13(l propagators, and using (|B . 14|) for the quark fields, the second line of 
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B.17(l . i.e., the first diagram of figure (5.3), reads 



d 4 xe^ I d*x I d^HAsfl.j^sl^dfW- 



P 2 9al3 - 2g s f abc G b 



a/3. 



d 4 x I d 4 y (x\ s^' 7 M 



(P + qfg a p - 2g s f^G b a 



(o|4°V 



(B.18) 



where in the second equality i|B.12() has been used. Now, expanding the propagators at large 
external momentum 



1 



-y 



i 



i 



i „ i 



i d j 
-rfi- fl- 



ip + qfg al3 - 2g s ^G\ 



mi 



4 s — gap 



p2 

~^9m3 



2g s f abc G a0 



<P 



(B.19) 



and keeping the 0(P ) terms we end up with the dimension six contribution. Colour-fierzing 
can be done with the aid of 



A" 



A" 



-r 1 I — I =-[8 a S jk - 



1 

N~c 



8i 3 S k i 



while Dirac algebra can be simplified with the following identity (in 4 dimensions) 



.0123 



(B.20) 



(B.21) 



From the equation above it is straightforward to derive the useful relations 



YYl X ® lulalx = 2<T CT 7 a ® 7ct + 27 CT ® 7„ + 25" a -j a l5 ® 7^75 - 27^75 ® 7^75 

Y l l u l X ® ixlalu = 2^7" ® 7« + 27 ff ® lv ~ 2,5 l/ (T 7 Q 75 ® 7a 7s + 27 ff 7 5 ® 7„75 (B.22) 

where use has been made of the well-known expression 

e^ Xp e^ aP = -2{5\5 p p - S x ^ a ) (B.23) 

The procedure can be straightforwardly extended to the remaining diagrams of figure (5.3). 
The extraction of dimension eight operators (0(P 2 ) in l|B.17|) ) is formally identical to the one 
for dimension six, except for the fact that now one has to deal with covariant derivatives and 
a bit more complicated Dirac structures. Both points arc illustrated in the next section. 
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B.4 Determination of dimension eight operators in the 
I /ml expansion 

We concentrate on the box diagram of figure (5.1) with, e.g., charm-charm exchange. The 
procedure can be readily exported to the remaining diagrams. The contribution from this 
diagram can be obtained in coordinate space from the following expression: 



G 



F M^\ 2 U I d i xd A yd A zd i w (x|s(l + 7 5 )7 /i - 



1 



+ ff - mi 



-7^(1-75)%} 



\z)(w\s{l+~f 5 )j„- 



q 2 -M 2 v ' f + ^ _ m , 



■7^(1 ~ l5)d\z) (y\- 



\w) (B.24) 



Integration of the W particle amounts to the expansion 

1 , , -1 



S(x — z) + 



(B.25) 



This has to be done twice. One of the Dirac deltas can be used to reduce the number of 
integrals, whereas for the second one we use its representation in momentum space: 



5{x - y) 



rl 4 n 

» 1 Aq-(x-y) 



(B.26) 



which will eventually give rise to the loop integral. 

Now we can expand the quark propagators to any order in the soft momentum f in the 
following way: 



1 



+ ff - mi 



i 



i - m, 



1 



1 



1 



i - mi fa - m, fa - m, t 

1 ft 1 1 



— mi 



(B.27) 



Recall that this is the massive counterpart of the first line of (|B.19|I . Keeping terms of 0(P 2 ), 
since we are interested in dimension-eight operators, and after some diracology, we end up 
with: 



(x\Pr 



1 



(x\P L 



{P-q? 



P L \0) = (x\P L 
2{P-q)i 



Pl\Q) 



mf (q 2 - mf) 2 



P 2 



(q 2 — mf) 3 (q 2 — mf) 2 (q 2 — mf) 2 (q 2 — mf) 2 



(P-q) _ .^KPpq^Sls 



(q 2 - mf) 2 



(B.28) 
Pl\0) 



which is the expression to be employed for the quark propagator hereafter. 
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Now we use IjB.lOj) to convert operators into covariant derivatives. Integration by parts 
shifts the derivatives to the quark fields and delta functions can then be integrated in a trivial 
way. Thus, one can write, for instance: 

f d 4 x d 4 y s L {x)-f (x\ -J—-\Q)^d L (0)s L (y) lv (y\ - ^ |0) 7Al d L (0) = 



= s L {0)D^ 7u d L {0)s L (0) lfi d L (0) 



(2tt)^ (g2 - m |)3 



where the Dirac structure has already been simplified and the operator in the second line has 
been rewritten in terms of the dual gluon field strength tensor 

& v = \e^G aP , G^ = yG;„, £ 0123 = +1 (B.30) 



In order to get to the result in Eq. (|B.29|I one must employ the identity D 2 = 4f cr Alt/ G A " y + If> , 



2 

and the equations of motion for the quark fields in the chiral limit. This will leave us with two 
operators, one with the gluon field strength and one with its dual. Out of them, the former is 
chiral-supprcsscd. Indeed, using 

SLG av lv d L = y {s L (D a B - % D a ) d L - 0„ {sLl„D a d L )} (B.31) 

one finds that 

(K°(p)\s L YG av d L \0) = 1 d^K°(p)\sL^D a d L \0) ~ 0( P 2 Pa ) (B.32) 

and, thus, it can be neglected. 

Upon performing the integral over q in Eq. I|B.29(I one can easily extract the dependence. 
After adding up all the contributions, one obtains the result l|5.89|) - H5.90|l in the main text. 



Appendix C 



Bernoulli Numbers and Bernoulli 
Polynomials 



C.l Bernoulli Polynomials 

A sequence of polynomials A n (x), where n is a positive integer, is called an Appell sequence if 
they satisfy the following properties 

deg(A n ) = n (C.l) 

A' n (x) = nA n - 1 (x) (C.2) 

The last equation provides a method for recursively determine the full Appell sequence, given 
an initial Aq. Therefore, in general 



M*)=Y,[ ™)"rX n ~ r (C.3) 

where the combinatorial numbers are defined as usual 

n\_ r(n+l) 

r J r(r + l)T(n-r + l) v ' ' 

As with most sequences, one can define the whole sequence through its generating function, to 
wit 

G( V ) = V4(x)7 (C.5) 

n=0 

Given the above definitions, the Bernoulli polynomials B n (x) can be defined as the Appell 
sequence with 

B (x) = 1 (C.6) 

which satisfies the additional constraint 

• i 

B n (x)dx = 0, n = l,2,„. (C.7) 

o 
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C.l Bernoulli Polynomials 
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It is relatively easy to show that their generating function is 

oo n 

e z — 1 ' n\ 

Quite notably, Bernoulli polynomials have a reflection symmetry relation 

B n (l-x) = (-l) n B n (x) 

and a related one 



B n (x + 1) — B n (x) = nx 



n-l 



(C.8) 



(C.9) 



(CIO) 



In particular, using the above equation for certain values of x it is easy to conclude that odd 
Bernoulli polynomials satisfy 



B 2n+1 (x)=0 , (7i = l,2,...) 



when x = 0, -, 1 



(C.ll) 



Sometimes it is also useful to express Bernoulli polynomials in terms of their trigonometric 
expansion 



iw.)-(-ir»(».)ig^E=^ 



k=l 



B 2n+1 (x) = (-1)" +1 2 (2n + 1)! ^ £ 2ra+1 

v 7 fe=i 



which follow from the expansion in Fourier series 



S n (x) = 

The first Bernoulli polynomials read 

B (x 
B x {x 

B 2 (x 
B 3 (x 
B 4 (x 
B 5 (x 
B 6 {x 



(2m) 



An 



,2-irikx 



(C.12) 
(C.13) 

(C.14) 



fe=-c 



1 

1 

X ~2 



x -x + - 
6 

3 3 2 1 

* ~ 2* + 2 X 



c 5 a 5 o 1 

x s X + -x J X 

2 3 6 



2 2 42 



(C.15) 
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C.2 Bernoulli numbers 

Bernoulli numbers B r are denned as the coefficients a r of the Appell sequence for the Bernoulli 
polynomials. Therefore 

Bn{x)=jr( n \B r x n - r (C.16) 

It is easy to show that indeed Bernoulli numbers can be determined from Bernoulli polynomials 
in a straightforward way, namely 

B n = B n (0) or analogously B n = (-l) n S„(l) (C.17) 



by using the reflection property (|C.9|) . This justifies the use of the same symbol for both 
Bernoulli polynomials and Bernoulli numbers. However, sometimes it is also useful to deter- 
mine the Bernoulli numbers without having to resort to Bernoulli polynomials. The following 
formulae arc then quite useful. Bernoulli numbers can be expressed as an integral of the 
generating function, 

Bn = ^t*J^ (C.18) 
2ni J e z - 1 z^ 1 v ' 

as derivatives of the generating function 

B„ = lim -f- ( —^—] (C.19) 

or in terms of the Riemann zeta function 

S n = (-l) n+1 nC(l-n) (C.20) 

Recalling (fTTTTT^ . we can immediately conclude that only even Bernoulli numbers are non- 
trivial, i.e. 

B 2n+ i = , (n = l,2,...) (C.21) 
For even Bernoulli numbers there exist the following representation 



(-l )"~ 1 2 (2n)! 
Wf 3 



B 2n = V C(2n) (C.22) 



whose asymptotic behaviour can be determined readily if one uses the Stirling asymptotic 
formula 



(71 \ 2n 
— ) (C23) 
7re/ 



C.2 Bernoulli numbers 
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The first Bernoulli numbers read 



B = 1 

* - -J 

Bi = ~30 

B6 = 42 

fto = ^ (C24) 
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